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Preface 


The proceedings of the Second Conference on Sonic Boom Research 
held at NASA headquarters, Washington, D.C., on May 9-10, 1968, 
are reported in this NASA Special Publication. This conference was 
organized by the Research Division of the Office of Advanced Re- 
search and Technology. It was, in essence, a follow-on meeting to 
the first Sonic Boom Conference, the proceedings of which are reported 
in NASA SP-147. 

The purpose of this conference was to review the current status of 
the NASA University Program on Sonic Boom Research, to survey 
the current research program at the NASA centers, to determine those 
areas of sonic boom research that are the most pressing from the 
standpoint of commercial supersonic transport (SST) operation, and 
to determine which, if any, of the various avenues of research appear 
to be the most promising with regard to sonic boom overpressure 
reduction. 

The entire first day and the morning of the second day of the meet- 
ing were devoted to the presentation of the invited papers covering 
the status of the various NASA projects. A. Richard Seebass served 
as chairman of these sessions. The afternoon session of the second 
day was devoted to an informal open discussion of the state of the 
art of sonic boom theory. The subjects discussed in this session could 
be categorized into two principal topics which serve to answer the 
following two questions: How far have we proceeded in developing 
methods of analysis and understanding of sonic boom phenomena? 
How far have we proceeded in devising methods of reducing sonic 
boom for feasible aircraft designs? Wallace D. Hayes served as 
chairman of this discussion session. 

As was the situation last year, and to reiterate a statement made 
in last year’s report (NASA SP-147), “it is difficult to overemphasize 
the importance of the phenomenon we call the sonic boom.” In July 
1968 the President of the United States signed into law a bill requiring 
the Federal Aviation Administration (FAA) Administrator to . ... pre- 
scribe and amend such rules and regulations as he may find necessary 
to provide for the control and abatement of aircraft noise and sonic 
boom. . . Sponsors of the legislation said it was broad enough to 
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allow the FAA to ban supersonic flights entirely unless the boom 
problem is brought under control. Several other countries are con- 
sidering regulations restricting proposed SST's to subsonic flight over 
populated areas. It is the consensus of opinion that the SST's must 
be capable of supersonic operation over land to be a true economic 
success. Such aircraft should be able to fly supersonically anywhere, 
over land as well as water, without causing human discomfort or 
structural damage. 

It is not certain yet how much sonic boom overpressure will be 
acceptable for humans. Present thinking is that the overpressures 
may have to be reduced to values as low as 1 pound per square foot 
to be acceptable. Thus, it is apparent that every aerodynamic tech- 
nique of reducing the overpressure must be investigated. A large 
number of investigations and significant progress have already been 
made by NASA in cooperation with the Federal Aviation Administration. 
Much more scientific understanding of the problem is needed before a 
large reduction in sonic boom is accomplished. 

As a result of last year's conference and after discussing the sonic 
boom problem with the FAA and the U.S. Air Force, the following 
10 areas were defined for further research on the generation and 
propagation of sonic booms: (1) validity of the basic theory, (2) mid- 
field signature modification, (3) atmospheric effects, (4) effects of 
accelerations and maneuvers, (5) multipole effects and lateral spread- 
ing, (6) pressure distributions on various geometrical configurations, 

(7) lift contribution to the far-field downward momentum flux, 

(8) sonic boom of exotic configurations, (9) engine stream-tube area 
effects, and (10) techniques for predicting the sonic boom of hyper- 
sonic aircraft. 

Also, it was concluded at last year's meeting, April 12, 1967, that 
there was a clear need for the active participation of knowledgeable 
aerodynamicists and engineers in sonic boom research in order to 
elevate our level of understanding. Accordingly, the Research 
Division of the Office of Advanced Research and Technology, NASA 
Headquarters, established the following seven research programs in 
six universities to complement the in-house program at the NASA 
Ames and Langley Research Centers and to investigate the 10 basic 
research areas : 

New York University.... A. Ferri and L. Ting Low boom configurations, 

N-wave reflections 

Cornell University. E. L. Resler, Jr., No boom-lifting con- 

and F. K. Moore figurations, far-field 

theory 

The Aeronautical Re- M. T. Landahl Nonlinear effects 

search Institute of 
Sweden. 
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Colorado University- A. Busemann Boom reduction 

Princeton University W. D. Hayes _ __________ Second-order wave 

theory, geometric 
acoustics, propagation 
through a caustic 
(ray envelope) 

Cornell University_______ A. R. Seebass and Azimuthal redistribution 

A. R. George 

Columbia University M. B. Friedman Theory of the superboom 

The progress made during the past year by the NASA centers and 
contractors represents substantial advances in developing methods 
of analysis and understanding of sonic boom phenomena and in 
devising methods, particularly aerodynamic techniques, of reducing 
sonic boom for practical aircraft designs. For example, in regard to 
the analysis and understanding of sonic boom phenomena, the straight- 
forward technique for predicting the propagation of a sonic boom 
from a slender aircraft can now be conveniently employed with the 
assumption that nonlinearity does not affect the rays. Use of such 
an assumption permits one to calculate directly the nonlinear modifi- 
cations of the signature. Great advances have been made in the 
understanding of geometric acoustics with winds and variable sound 
speed. Also a linearized theory for flow fields around aircraft has 
been well known for many years. However, life is not all that simple 
in view of the fact that the linearized theory referred to above does 
not provide a dependable and consistent first-order description of the 
flow field far from the aircraft. In addition, nonlinear effects exist 
near the aircraft. 

For practical considerations and importance the following theoret- 
ical problems have been identified, and though better understood, 
remain unsolved : 

(1) Prediction of sonic boom signatures near a caustic 

(2) Nonlinear effects on ray tubes 

(3) Diffraction into shadow zones 

(4) Nonlinear effects near the aircraft 

(5) Effects of turbulence 

Our researchers believe that they can find practical solutions to 
these basic problems because the methods of solution are presently 
within focus. 

During the past year, the Aeronautical Research Associates of 
Princeton under contract to NASA developed a computer program 
that will compute the full overpressure signature anywhere within 
the field for an aircraft in unsteady supersonic flight in a horizontally 
stratified atmosphere with steady winds. Although the program 
cannot account for the vertical winds, these effects are usually small 
compared with those of the horizontal winds. 
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A major portion of the meeting was devoted to reporting the 
progress that was made during the past year regarding the develop- 
ment of aerodynamic techniques to reduce the sonic boom. Foremost 
in the minds of the researchers is the well-recognized fact that techno- 
logical advancements that produce improved overall efficiency for 
SST aircraft, such as improvements in engine thrust to weight ratio, 
lift-to-drag ratio, structural weight, and fuel consumption, will reduce 
sonic boom intensities. Above all, the reduction of the sonic boom is 
not to be realized at the expense of aircraft efficiency. 

Several papers presented at this meeting discussed promising 
aerodynamic techniques for reducing the sonic boom. Antonio 
Ferri showed that redistribution of lift along the length of present 
SST configurations results in near-field effects that provide maximum 
sonic boom overpressure of the order of 1 pound per square foot at 
no expense to the lift-to-drag ratio of the airplane at cruise. This 
is a significant reduction compared with the maximum boom over- 
pressure of 2 pounds per square foot for present configurations selected 
for the supersonic transport. For a given lift and volume of an 
aircraft there is a lower bound to the far-field overpressure. This 
results from the fact that, while the lift contribution to sonic boom 
can be reduced by the addition of suitable volume elements, there is 
no way of avoiding the boom attributed to lift alone without modify- 
ing the Bernoulli constant of the flow. Thus, Ferri has attempted 
to approach these lower bounds with aircraft configurations that 
have realistic skin friction and wave drag contributions and not 
unusual structural requirements. 

A. R. George’s paper has shown how quadrupole effects can also 
be used to approach this theoretical lower bound. Other interesting 
boom reduction techniques currently being studied are discussed by 
Resler, Seebass, Hayes, Carlson, and McLean in this volume. 

In an assessment of the status of our sonic boom research program 
at this meeting, it was concluded that there are no carefully docu- 
mented experimental results for the overpressure signature in the 
vicinity of a caustic that results from an aircraft acceleration maneuver 
or an ambient sound speed gradient. Because no theory presently 
exists for the overpressure signature at a caustic and because the 
pressures are greatly intensified there, it appears that much more 
experimental work is needed. It should be possible to simulate the 
maneuver effects in ballistic ranges. The ambient sound speed 
gradient effects can also be duplicated in shock tubes, supersonic 
wind tunnels, and certain simulators. The NASA Langley Research 
Center has recently obtained some qualitative information on the 
effect of sound speed gradient in a ballistic range. Further experi- 
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mental work on this problem is being planned by several groups in 
the NASA contract program to complement the in-house program. 

Finally an assessment of the proceedings of this meeting has 
indicated to me and to those researchers who participated in the 
meeting that significant progress has been made during the past 
year in our understanding and analysis of the generation and prop- 
agation of sonic boom. Further, it has been analytically demon- 
strated in the NASA program during the past year that the utilization 
of a particular aerodynamic technique can result in large reductions 
in sonic boom overpressures. This is especially significant in that 
it allows us to expect vast improvements in future SST conceptual 
designs from both improvements in aircraft operating efficiency as 
well as through novel aircraft designs. It is yet too early to predict 
whether any of these design techniques will lead the way to the 
development of a domestic SST that will be allowed to fly super- 
sonically over land as well as over water. It is to be expected that 
more definitive answers will be forthcoming as the various promising 
avenues of approach are explored in much greater detail. 

I. R. Schwartz 

Research Division 

Office of Advanced Research and Technology 
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Atmospheric Effects on the Sonic Boom 

I. Edward Garrick 
Langley Research Center , NASA 


INTRODUCTION 

The material presented at the conference was based to a large 
extent on a paper presented at the International Association of Me- 
teorology and Atmospheric Physics, XIV General Assembly of the 
International Union of Geodesy and Geophysics, Lucerne, Switzer- 
land, and published recently as NASA Technical Note D-4588 (ref. 
1) . For this reason, only a condensed discussion of this material is 
given here, and some newer information will be stressed. 

Acceptance by the general public of supersonic transports will 
depend in great measure on the characteristics of their sonic booms 
(ref. .2). These characteristics are affected by several factors which 
depend on design, operations, and environment or atmospheric 
conditions. Design factors such as weight, size, lift, and volume dis- 
tributions and operation factors such as altitude, Mach number, 
and flight path may be accounted for, or even predicted, by existing 
theoretical methods which have been closely verified under con- 
trolled conditions, such as in a supersonic wind tunnel. However, 
the effects of environment or atmospheric conditions, such as wind, 
temperature, and 'pressure gradients and turbulence, like the predic- 
tion of weather itself, remain probabilistic or statistical in nature. 

The report aims to summarize and discuss some of these atmospheric 
effects, including the findings of several flight studies. Much of the 
flight information has already been published; however, in the paper 
by D. J, Maglieri in this volume entitled “Sonic Boom Ground 
Pressure Measurements for Flights at Altitudes in Excess of 70 000 
Feet and at Mach Numbers up to 3.0,” some of the most recent data 
•are included. It is hoped herein to account for or explain, at least 
qualitatively, the kinds of statistical information on peak pressures 
that are found by field measurement. 
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ANALYTICAL CONSIDERATIONS 

Only a broad overview of the theoretical background is mentioned 
here. Prediction of sonic boom characteristics has been based on 
methods in which idealized atmospheres are assumed. For example, 
in his basic work on the flow fields of bodies of revolution, Whitham 
assumed a uniform atmosphere. Randall employed a correction factor 
to account for an isothermal hydrostatic (exponential) atmosphere. 

Effects of a layered or stratified atmosphere with arbitrary steady 
or slowly varying winds and temperatures were studied by Palmer, 
Friedman, Kane, and Sigalla. This work makes direct use of classical 
concepts of geometrical (ray) acoustics and permits comparison of 
results for a nonstandard atmosphere with those for the standard 
atmosphere. It also leads to means for predicting effects of maneuvers 
and explains in large measure the focusing and channeling of energy 
of acoustical signals. Many references are given in reference 3 by 
Kane and Palmer, a report which contains the most detailed study 
of meteorological effects on the sonic boom, based on an analysis of 
NASA data obtained prior to the summer of 1964, (See also ref. 4.) 

It is pertinent to mention briefly some recent work on propagation 
in a stratified atmosphere sponsored by the Langley Research Center 
under contract and conducted by Aeronautical Research Associates 
of Princeton, Inc, This work, to be published, is authored by Wallace 
D. Hayes, Rudolph C. Haefeli, and Helene Kulsrud. It presents a 
theory, based on geometric acoustics, and a calculation program for 
tracing the pressure signature from the vicinity of the airplane to 
the ground, as influenced by passage through a steady layered at- 
mosphere. The program will give the location, intensity, and time 
history for both level flight and for maneuvers. See also a paper en- 
titled “Similarity Rules for Nonlinear Acoustic Propagation Through 
a Caustic” by Professor Hayes in this volume. 

Turbulence and temperature fluctuations are known to have im- 
portant modulating effects on the distortion of the sonic boom pressure 
signatures. Theoretical considerations are needed to explain anomalous 
spiky signatures and to rationalize the statistical data. It may be 
noted that a large literature exists relevent to sinusoidal sound (or 
light) with substantial contributions by authors from the U.S.S.R. 
For an account of this work which starts with contributions of A. M. 
Obukov, see the books by Tatarski and Chernov (refs. 5 and 6). In 
particular, the method of Rytov, combined with the random variable 
properties of patches of turbulence, has led to log normal charac- 
teristics for the statistics of the scattered pressure amplitudes. Experi- 
mental measurements of Kallistratova for sound scattering and of 
Tatarski for light scattering tend to confirm trends indicated by 
theory. 
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Many others have contributed to the theoretical background of 
the interaction of turbulence with sound or shock waves, and only a 
few names are mentioned here: Ligh thill, Ribner, Palmer, Muller, 
and Matschat. Despite the encouraging analytical efforts, relevant 
statistical predictions for sonic booms at present seem to depend 
mainly on measurements from airplane flights. In presenting the 
measured data subsequently, it is convenient in general to normalize 
the measured data in relation to a nominal overpressure calculated 
from the simple theory of Whitham or Randall, or from that of a 
layered standard atmosphere. One such scheme programed for com- 
puting machines is given in reference 7. 

A recent significant contribution to the theory of effects of atmos- 
pheric turbulence on sonic booms by S. C. Crow was made known to 
the author through a preprint (ref. 8). It is believed appropriate to 
include, herein the following complete summary from this report. 

Recorded pressure signatures of supersonic airplanes often show intense, 
spiky perturbations superimposed on a basic N-shaped pattern. A scattering 
theory, incorporating both inertial and thermal interactions, is developed to 
explain the spikes. Scattering from a weak shock is studied first. An exact 
solution of the scattering equation, finite at the shock and everywhere behind 
it, is derived in the form of a surface integral over a paraboloid of dependence, 
whose focus is the observation point and whose directrix is the shock. The 
solution is found to degenerate at the shock into the result given by ray acoustics. 
Eddies in the Kolmogorov inertial subrange are found to be the main source of 
sonic-bang spikes, and Kolmogorov's similarity theory is used to show that, 
for almost all separations h between the shock and the observation point, the 
mean-square pressure perturbation equals (Ap ) 2 ( hjh ) 7/ ®, where A y is the 
pressure jump across the shock and h c is a critical distance predicted in terms of 
meteorological conditions. The mathematically predicted mean-square pressure 
perturbation at the shock itself is found to be finite but enormous. The exact 
solution of the scattering equation is generalized for arbitrary waveforms, and 
the approximate relation G(f) — G N (J) (1 + (/// c ) 7/6 ) is established between the 
energy spectrum G N ( f) of undeformed sonic bang and its mean perturbed spectrum 
G(f). The critical frequency f CJ inversely proportional to h C} is found to be as 
low as 63 cps for an idealized model of atmospheric turbulence. The conclusion 
is reached that scattering can augment considerably the psychological impact 
of a sonic bang. 


ATMOSPHERIC TURBULENCE 

In recent years, much information on the local structure of the 
atmosphere has been assembled (for example, ref. 9). It is known that 
characteristic lengths of turbulence in the lower atmosphere below 
1000 to 2000 feet (or about 305 to 610 meters) depend on height above 
terrain. It is also believed that much of the distortion of sonic boom 
pressure signatures attributed to turbulence may occur at the lower 
altitudes. An interesting detailed view of the vortex structure near 
the ground is afforded by figure 1 which is from ref. 10. This figure 
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Figure 1 . — Gust structure near the ground as indicated by isotachs (ref. 9) . 


shows equal velocity contours, or isotachs, obtained from measure- 
ments taken during 25 seconds of gusty weather. The horizontal 
section shows the isotachs obtained from measurements taken at 
eight stations along an array of 50-foot (15.2-meter) poles equally spaced 
over a linear distance of 420 feet (128 meters). The vertical section 
shows the isotachs obtained from measurements taken at five equally 
spaced stations up a 250-foot-high (76.2-meter) tower. In each sec- 
tion, the ordinate represents distance along the array or up the tower, 
and the abscissa is time. During the sampling period the average 
wind traversed a distance of over 1000 feet (304.8 meters). 

Turbulence in the upper atmosphere has been probed by instru- 
mented airplanes (ref. 11). Typical results of measurement for 
moderate turbulence are shown in figure 2 as the variation of power 
spectral density of vertical velocity with wave number (1/wave- 
length). The dashed extension to the measurements follows the 
theoretically indicated 5/3 power of the wave number. The shape 
of this curve tends to be relatively invariant while the amplitude is 
a function of the severity of the turbulence. The curve represents 
the Kolmogorov inertial subrange; for lower wave numbers than 
those shown, the dynamical processes originating the turbulence 
dominate and influence the spectral shape. For higher wave numbers 
than those shown, dissipative processes converting turbulence to heat 
begin to dominate. The curve shows the distribution of power in 
the various wavelengths and the area under the curve yields the 
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Figure 2. — A spectrum of atmospheric turbulence (ref. 10) . 

total power. This same type of figure can also apply approximately 
near the ground with relatively less power in the longer wavelengths. 
From its spectral distribution, it may be determined that a normal 
sonic boom pressure signature of 0.3- to 0.4-second duration covers 
acoustic wavelengths that range from about 1000 meters to fractions 
of a meter. Thus, the interaction of the pressure signature and the 
turbulence occurs over the entire spectrum shown and includes 
acoustic wavelengths less than and greater than the characteristic 
vortex sizes. The turbulence of smaller scale than the sound wave- 
lengths tends to produce isotropic scattering, whereas the scattering 
of larger scale turbulence is directional. 

FILTERING THE N-WAVE 

Insight into subjective effects of sonic booms and into frequency 
response requirements of measuring instruments is obtained by study 
of the amplitude and frequency content of the N-wave of various 
durations and rise times, as has been done. Further insight has 
been obtained by applying various filtering techniques to the analog 
signal of the N-wave. Two typical results are shown in figure 3: 
(a) filtering so that frequencies below 10 hertz are not measured, 
and (6) filtering out all frequencies above 30 hertz. The amplitude 
is not changed very much, but, whereas filtering out the high fre- 
quencies yields a rounded N-wave, filtering out the low frequencies 
gives a sharp rise and a distorted wave. Thus, the high frequencies 
315-766 0—68 2 



8 


SONIC BOOM RESEARCH 


TYPE OF INPUT 



(o) MEASURED RESPONSE FOR 

FREQUENCIES FROM 10 Hz TO 10000 Hz 




Figure 3. — Effects of filtering low or high frequencies from an N-wave. 

seem dominant in defining the rapid compressions of the N-wave, 
and the lower frequencies appear dominant in defining the expansion. 
However, while filtering out high frequencies suffices to obtain rounded 
signatures, it is not a necessary requirement, because these signatures 
may also be obtained by a shifting of the phases of the high-frequency 
components, as will be further discussed. 

MEASURED DATA 
Wave Shapes 

An indication of the measured variations in wave shapes at ground 
level is given in figure 4 for airplanes of three different sizes. The 
associated durations are approximately 0.1, 0.2, and 0.3 second. 
These signatures generally fall in several classes: 

Normal representing signatures closely resembling theoretically 
calculated N -waves 

Peaked representing signatures wherein the pressure peak is 
amplified relative to the basic N-wave 
Rounded representing signatures with longer rise times and 
lower peaks than the N-wave 
Combinations and variations of these signatures also occur. 

It is known that near the airplane the pressure signature is quite 
complex, regardless of atmospheric effects, and that with increased 
distance the nonlinear characteristics of the propagation process in 
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F-104 B-58 XB-70 



ROUNDED 

Figure 4, — -Variation in measured sonic boom pressure signatures at ground 
level for small, medium, and large airplanes in steady level flight (ref. 2). 



Figure 5. — Indication of signatures measured only a few thousand feet from the 
XB-70 airplane in flight and that at the ground (ref. 2). 

quiescent atmosphere tend to simplify the signature and to produce 
the N-wave signature. Figure 5 shows results of probe flights of the 
flow field of the XB-70 airplane flying at M=1.5 and at an altitude 
of 37 000 feet (11.28 kilometers) made by an F-104 airplane at 
various separation distances. The detailed geometry of the airplane 
is reflected in the near-field signatures both above and below the 
airplane. Of particular interest is the signature at ground level 
which indicates that for an altitude of 37 000 feet the far-field N-wave 
is not fully achieved for this large airplane, rather a hybrid-type, 
near-field signature exists. 
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Sonic Boom Pressures 

Measured sonic boom pressure signatures for two flights of the same 
airplane type are shown in figure 6. The pressure signatures on the 
right were obtained when surface winds were about 28 knots (14.4 
m/sec) and rather gusty, and those on the left were obtained for much 
lower wind velocities. These signatures were obtained at one station 
by microphones separated by only 100 feet (30.48 meters) in a cross- 
arrangement. However, as indicated in reference 11, similar varia- 
tions, indicating wave-shape distortion on the one hand, and lack 
of distortion on the other hand, were measured over a 150-square-mile 
(388.5-km 2 ) area during the same two flights. 

Statistics of Pressure Peaks 

A number of flight programs during which routine measurement of 
sonic boom signatures was made with carefully calibrated pressure 
instrumentation have been conducted and reported. The most ex- 
tensive of these was a 6-month flight program at Oklahoma City 
(ref. 13). However, only fragmentary observations were made of 
associated weather conditions. The investigation demonstrated that 
wide variations in the signature types were associated with dynamics 
of the atmosphere; moreover, the distributions of pressure peaks 
appeared to follow approximately a log normal probability curve — * 



(a) LOW WIND VELOCITIES. (b) HIGH WIND VELOCITIES (GUSTY). 


Figure 6.— Time histories of sonic boom overpressures showing wave-shape 
variations measured for two flights of a B-58 airplane on different days (ref. 11). 
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that is, if the pressure peaks are plotted as log pressure or in decibels, 
an approximately Gaussian or normal distribution occurs. Also 
shown by this study was the result that the impulse also followed such 
a distribution law, although it has smaller variance or standard 
deviation. The approximate theory for scattering of sound (mono- 
chromatic waves) by turbulence also yields a normal distribution 
for the logarithm of the amplitude (ref. 5). The theory indicates 
that given certain sufficient information about the structure of turbu- 
lence, a calculation of the variance of the distribution of peaks for 
sonic booms may be feasible. Palmer has already made some 
preliminary estimates in this regard. 

In conjunction with the XB-70 flights a number of accompanying 
flights of B-58 and F-104 airplanes were made. The Mach number 
and altitude ranges covered were from M— 1.5 to 2.5 and h— 11.3 to 
about 18 kilometers for the XB-70 airplane, M~ 1.5 to 1.65 and 
h= 9.7 to 12.2 kilometers for the B-58 airplane, and M— 1.3 to 1.4 
and h— 5.2 to 6.3 kilometers for the F-104 airplane. The nominal 
calculated overpressure for all these airplanes was about 1 millibar. 
The data in figure 7 for the three airplanes were obtained in late 
morning flights from November 1966 to January 1967, when generally 
small turbulence and convectivity occurred. Despite the large 
differences in airplane operating conditions and signature durations, 
a strikingly similar probability distribution pattern for the three 
airplanes appears (the log normal distribution is linear on this type 
of plot) . 

Tabular Summary 

Table I provides a convenient summary of some of the results of 
several flight programs (refs. 12 to 14). Presented are the ratios of 
the mean values measured to the nominal values calculated. The 
ratios of the ranges of ±1 <r and ±2 a- relative to the mean of the 
pressure ratio are also given. (These ratios may be regarded as corre- 
sponding to a chosen mean pressure ratio of unity.) The ±3<r values 
as extrapolated are also presented and represent 99.7 percent of the 
data. This latter result is only suggestive and, as mentioned in refer- 
ence 3, these extreme values may be conservative. The pressure meas- 
urements were made under the flight track, except where the lateral 
distance is noted. The largest variability in the ratio factors occurs at 
the farthest lateral distances and for the lowest supersonic Mach 
numbers where ray paths in the atmosphere are longest. However, 
the nominal calculated values are generally lower at the lateral posi- 
tions because of the greater distances. The numerical factor of the 
mean overpressure for the ±2 a range of ground-track data varies 
from a maximum of 1.78 (Oklahoma City data) to a minimum of 1.22 
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Table I. — Statistical Variations oj Sonic Boom Overpressures 


Test site 

Time period 

R ange 

Mean of 
A P°„.lo 

Ratios relative to mean for range 

Number 

data 

samples 

Mach no. 

Altitude, km 

±l<r 

±2cr 

±3ff 

Oklahoma City 1,,_ 

Feb. to Apr. 1964__ 

1. 3 to 2. 0 

6. 4 to 12. 5 

* 0. 825 

0. 76 to 1. 30 

0. 58 to 1. 70 

0. 45 to 2. 21 

652 

Oklahoma Oity 3 b 




* 1. 04 

. 73 to 1. 37 

. 53 to 1. 86 

. 39 to 2. 55 

654 

Oklahoma City 4 ° , 




* 96 

. 67 to 1. 51 

. 45 to 2. 23 

. 30 to 3. 31 

637 

Oklahoma City 1 

May to July 1964— 

1. 2 to 1. 6 

8. 5 to 14. 0 

* 925 

. 75 to 1. 34 

. 56 to 1. 78 

. 42 to 2. 37 

549 

Oklahoma City 3 b 




* 1. 20 

. 74 to 1. 35 

. 56 to 1. 92 

. 42 to 2. 37 

554 

Oklahoma City 4 c _ _ 




1. 38 

. 68 to 1. 49 

. 46 to 2. 39 

.31 to 3. 23 

548 

Edwards Air Force 

! Nov. 1966 to Jan. 

1. 3 

9. 4 

4 855 

. 78 to 1. 28 

. 61 to 1. 64 

. 47 to 2. 11 

1378 

Base. 

1967. 








Chicago, _ „ 

Jan. to Mar. 1965, 

1. 2 to 1. 66 

11. 6 to 14. 6 

4 92 

. 77 to 1. 29 

. 60 to 1. 67 

. 46 to 2. 15 

191 

Edwards Air Force 

Sept, to Oct. 1961 _ 

1. 5 to 2. 0 

9. 4 to 21. 3 

d 1. 035 

. 90 to 1. 11 

. 81 to 1. 23 

. 73 to 1. 36 

135 

Base. 









Combined (fig. 


1. 2 to 2. 0 

9. 4 to 21. 3 

4 942 

. 82 to 1. 24 

. 67 to 1. 53 

. 54 to 1. 87 

358 

14). 

! 








Edwards Air Force 

Nov. 1966 to Jan. 

1. 5 to 2. 5 

11. 3 to 18. 3 

d 1. 00 

. 91 to 1. 10 

. 82 to 1. 22 

. 74 to 1.36 

447 

Base. 

1967. 








Edwards Air Force 

Nov. 1966 to Jan. 

3. 0 

21.3 + 

.a. 

CD 

00 

. 86 to 1. 16 

. 74 to 1. 34 

. 64 to 1. 56 

578 

Base. 

1967. 









a Nominal overpressure calculated for isothermal atmosphere. 
b 5-mile offset from flight track, 
o 10-mile offset from flight track. 

d Nominal overpressure calculated for standard atmosphere. 
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No. No. DATA 



Figure 7. — Probability of exceeding a given value of the ratio of measured 
to calculated ground overpressures along the flight track for the XB-70, B-58, 
and F-104 airplanes. (Time intervals of flights of the three airplanes varied 
from about 2 to 5 minutes.) 

(Edwards winter data). The data in the table cover a wide range of 
Mach numbers; in particular, data at M= 3 obtained at Edwards 
Air Force Base are included. 

A histogram of combined data is shown in figure 8, representing 
over 4000 data samples taken on ground track for various aircraft. 
Because the calculation schemes for A po c&u differed slightly, this 
figure is basically only of qualitative interest. 

Measured Overpressure Signatures Along Range 

In one of the recent studies with the F-104 airplane, an 8000-foot 
horizontal (2.44-kilometer) range instrumented with many microphones 
was employed. The pressures were measured at about 40 stations 
along the range. Figure 9 shows some results for the overpressures, 
which indicate a remarkable wavelike pattern. Apparently long 
waves in the atmosphere of the order of 10® meters are influencing 
the pattern of the results. Experimental results suggest that this is 
a moving pattern and that the approximate Gaussian distribution 
pattern for the overpressures would exist at each station. A gradual 
change along the range from peaked to rounded signature is apparent. 
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PROBABILITY 

DISTRIBUTION 



Figure 8.— Histogram of combined data taken along ground track for various 

aircraft. 



Figure 9. — Overpressures as a function of distance on the ground track for 
F-104 airplane in steady flight at a Mach number of 1.3 and an altitude of 
30 500 feet, with two sample signatures. 


The measured positive impulse also shows a wavelike pattern, but 
considerably reduced in variance. 

The mechanism of the apparent embedding of the local pattern of 
turbulence within a systematic larger pattern is not understood. It 
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is of interest to examine the spectra of two signatures, one having 
a peaked and the other a rounded wave shape as indicated in figures 
9 and 10. The stations for these measured signatures were 600 feet 
(183 meters) apart. The calculated amplitude squared, or energy, 
spectra of the signatures are shown in figure 10. Phase plots were 
also obtained but are not reproduced herein. Only relatively small 
changes occur in the envelopes of the amplitudes, despite the large 
change in signatures. As indicated previously in the discussion of 
figure 3, differences between peaked and rounded signatures may be 
reflected in the presence or absence of high frequencies in the asso- 
ciated spectra. For the wave shapes and spectra of figure 10, the 
higher frequencies are apparent for both signatures; but calculations 
of the relative phases showed that the lower frequencies of the spectra 
appear well correlated and coherent, whereas for the higher frequencies 
the relative phases of the two waves tend to become random. The 
general theory of sound scattering also indicates reduced coherence 
of sound waves and increased fluctuations in phase with increase in 
frequency. After passage through uncorrelated patches of turbulence, 
the higher frequency components of the N-type wave tend to become 
random both in phase and amplitude. This behavior is also in 
accord with the central limit theorem in probability. Synthesis of 
waves with N-wave amplitude distribution and random phase distri- 
bution analytically should yield interesting statistics for the peak 
amplitudes. It is thus suggested that the turbulence scattering 
leads to a phase scrambling process taking place for the higher 
frequencies. 



Figure 10. — Energy spectra for the two different sonic boom pressure signatures 
in figure 9. (The relative amplitude is given by 10 logi 0 //(<o)/ 2 .) 
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CONCLUDING REMARKS 

The paper reviews pertinent information obtained in recent years 
relating to atmospheric effects on the sonic boom and, in particular, 
includes some results of various flight programs. These atmospheric 
effects are complex, and a statistical approach appears necessary. The 
statistics of peak pressures follow approximately a log normal distri- 
bution, a result that is indicated by existing theory for pure (sinusoidal) 
sound. A tabular summary of the flight data gives the standard 
deviations of pressure peaks relative to nominal calculated values of 
the mean. Information is included on observed variations of sonic 
boom signatures for different types and sizes of airplanes. The 
spectral content of some ideal and some measured pressure signatures 
is exhibited and discussed with reference to peakedness or roundedness 
of the wave. 
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Sonic Boom Ground Pressure Measuremenfs (or Flights at Altitudes 
in Excess of 70000 Feet and at Mach Numbers Up to 3.0 

Domenic J. Maglieri 
Langley Research Center , NASA 


INTRODUCTION 

During the past several years a substantial amount of sonic boom 
flight research data have been obtained for a wide range of operating 
conditions. (See Garrick’s paper “Atmospheric Effects on the 
Sonic Boom” in this volume and ref . I.) The material of this paper 
relates particularly to the effects of extreme altitude and represents 
the results of some recent sonic boom signature measurements asso- 
ciated with operation of the SR-71 airplane. The data were obtained 
during the Edwards Air Force Base Sonic Boom Evaluation Program 
conducted in 1966-1967 time period (ref. 2). About 2000 data points 
were recorded from 35 flights at altitudes in excess of 70000 feet and 
at Mach numbers to 3.0. 

Variations of the measured ground pressure signatures for on track 
and for lateral measuring station locations are illustrated and the 
statistical variations of the overpressure, positive impulse, wave 
duration, and shock wave rise time are presented. The above data are 
useful directly for evaluating the effects of altitude and Mach number 
on the sonic boom ground exposure patterns. The results of these 
studies will also permit validation of available theories regarding 
aircraft source generation and the propagation of shock waves from 
high altitudes. 

WAVEFORM VARIATIONS 

Three sample tracings of SR-71 sonic boom pressure time histories 
measured along the ground track are presented in figure 1. These 
signatures represent the range of wave shapes observed during these 
studies, and they are noted to vary from a normal N-wave, to a peaked 
wave, and to a rounded wave. This type variation is similar to those 
observed previously for the F-104, B-58, and XB-70 aircraft (ref. 1). 
The signatures of figure 1 are also used to define the quantities: 
overpressure, A p 0 , positive impulse, Z 0 , rise time, r, positive time 
duration, A to, and total time duration, AT. 
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Figure 1.— Tracings of measured SR-71 sonic boom pressure- time histories along 

ground track showing some waveform categories and signature definitions. 

GROUND TRACK MEASUREMENTS 

Figure 2 includes tracings of sonic boom pressure signatures 
measured along the ground track of the airplane at various altitudes 
from below 50 000 feet to in excess of 70 000 feet and Mach numbers to 
about 3.0. It can be seen that all the waveforms presented are of the 
N-wave type and are similar in nature to those obtained for the B-58 
aircraft of reference 3. In general, the peak-pressure values decrease 
and the time durations of the waves increase as the altitude of the 
aircraft is increased. Based on other related analytical work, it is 
estimated that the sonic boom signature, overpressure, and time period 
variations are attributable mainly to altitude rather than Mach 
number effect. 

The bow wave peak overpressures obtained from signatures of the 
type shown in figure 2 are presented in figure 3 as a function of air- 
plane altitude for steady level flight. The data of figure 3 were 
measured at distances within 3 nautical miles of the aircraft ground 
track. Each symbol represents an individual mission, the blocked-in 
symbols being the average of the measurements obtained on 20 to 64 
-microphones. The circle symbols represent data from flights made 
during the November 1966 to January 1967 time period for which 
atmospheric conditions were observed to be somewhat more quiescent 
than -those during the six flights of the June 1966 time period as 
represented by the diamond symbols. The plot of figure 3 repre- 
sents results of 28 flights for , a total of 928 data samples. 

At the present time the calculated sonic boom characteristics of 
the SR-71 aircraft are not available. However, predicted and meas- 
ured overpressures (ref. 4, fig. 16, and ref. 3, fig. 30) are available 



SONIC BOOM GROUND PRESSURE MEASUREMENTS 


21 


A^.lb/ft* 




70,000 



Figure 2.—' Tracings of SR-71 sonic boom signatures measured along ground track 
for aircraft at various altitudes and Mach numbers. 


*Po 

lb/ft 2 


• 20 to 64 MIKES (WINTER DATA) 
o 3 to 5 MIKES (WINTER DATA) 

♦ 20 td 64 MIKES (SUMMER DATA) 



Figure 3 “Measured peak overpressures along ground track as a function of 
altitude for airplanes in steady level flight. Data are for a range of Mach 
numbers to 3.0. 

for the B-58 aircraft, which is about the same size and weight as the 
SB-71, and the above B-58 data are shown as the stippled band in 
figure 3 for comparison. 

It can be seen from figure 3 that the SR-71 measurements are in 
general agreement with the measurements and calculations for the 
B-58 operating in the same altitude range. This result would suggest 
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that no unusual phenomena are observed because of the extreme 
altitude flight conditions of these tests. 

Since average values are plotted in figure 3, the data points do not 
represent the extreme variations of the measurements. Variations of 
the individual measurements will be illustrated in subsequent figures. 

LATERAL SPREAD MEASUREMENTS 

Figure 4 contains tracings of sonic boom signatures at various 
lateral distances out to about 27 miles from the ground track for 
flights at an altitude of about 70 000+ feet and a Mach number of 
about 3.0. As seen in figure 2, the pressure signatures have the gross 
features of N -waves except for extreme lateral distances. The rise 
times of the signatures of figure 4 become longer and the overpressures 
decrease as the lateral distance from the flight track increases. The 
wave period, however, does not seem to vary in a systematic manner 
with increasing distance as was the case for increasing altitude in 
figure 2. These findings are similar to those obtained for other 
aircraft (refs. 3 and 5). 

Overpressure data from signatures of the type shown in figure 4 
are plotted in figure 5 as a function of lateral distance from the flight 
track for two altitude ranges. The solid curves shown in figure 5 are 
the predicted overpressure variations for the B-58 aircraft (ref. 4). 
The dashed vertical lines are estimates of the lateral extent of the 
pressure patterns based on atmospheric refraction (ref. 6). The 

LATERAL DISTANCE, n.mi. 

■ - 0.09 


Ap 0 , lb/ft 2 





0.16 



[—0.10 see-) 


26.8 


Figure 4. — Tracings of SR-71 sonic boom ground pressure signatures as meas- 
ured at various lateral distances from the ground track for the aircraft at an 
altitude of 70 000+ feet and a Mach number of about 3.0. 
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(a) M = 2.0 @ 50 , 000 * (2 FLTS) 



LATERAL DISTANCE FROM TRACK, nmi 
•(b) M = 3.0 @ 70,000*+ (15 FLTS) 


Figure 5. — -Sonic boom overpressures for the SR-7I aircraft as a function of 
lateral distance for two different flight conditions. 

overpressures are seen generally to be a maximum along the ground 
track and to decrease with increasing lateral distance. The predic- 
tion of the lateral cutoff point appears to correlate well with the data 
for the Mach number and altitude ranges of these tests. 

STATISTICAL VARIATIONS 

Analyses of the SR-71 data have been made to determine the 
statistical variations of peak overpressures, A p 0 , positive impulse 
values, I Q , time durations, A t 0 and A T, and rise times. These data 
are presented in the form of histograms in figures 6 through 9 for all 
flights in the altitude range 70000+ feet and for measurements made 
within ±3 miles of the aircraft ground track. 

Figure 6 contains a histogram showing the variation of the measured 
ground pressures. The results shown are for 22 flights giving 704 
data samples. The overpressure interval for plotting was taken as 
0.05 psf. It can be seen that the largest number of events is associ- 
ated with overpressure values of the order of 0.9 psf. Values ranged 
from about 0.55 psf to about 2.0 psf, and one extreme value of about 
2.6 psf was encountered. 

A presentation similar to that of figure 6 is shown in figure 7 for 
the positive impulse values. An examination of the histogram of 
figure 7 indicate that the variability in the impulse is generally less 
than for the associated overpressures of figure 6, and this result is 
consistent with the findings of other similar studies. 

Figure 8 includes histograms showing the variability of the total 
time durations, AT, and the positive time durations, A£ 0 , of the 
SR-71 sonic boom signatures. It can be noted that variations of the 
order of ±25 percent exist, and these results are consistent with 

315-766 0—68 3 
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Figure 6. — Histogram showing variability of peak bow wave overpressure 
of SR-71 sonic boom signatures for aircraft operating at 70000+ feet and 
Mach number about 3.0. Data are for measurements within ± 3 nautical miles 
of ground track. 



Figure 7. — Histogram showing variability of positive impulse or SR-71 sonic 
boom signatures for aircraft operating at 70000+ feet and Mach number about 
3.0. Data are for measurements within ± 3 nautical miles of ground track. 
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at, sec Alb, sec 

Figure 8. — Histograms showing variability of the time durations of SR-71 
sonic boom signatures for aircraft operating at 70000+ feet and Mach number 
about 3.0. Data are for measurements within ±3 nautical miles of ground 
track. 

similar data shown for the B-58 in reference 2. Variations of the 
positive time durations and the overpressures can both cause varia- 
tions in the impulse; however, inspection of figures 6, 7, and 8 suggests 
that the overpressure variation is the dominant factor. 

The data of figure 9 relate to the rise times of the waves. Rise 
time data per unit overpressure are plotted in histogram form using a 
plotting interval of 0.001. For the types of waves illustrated sche- 
matically in figure 1, the rise time is defined as the time from the 
onset of bow shock overpressure to the time of the maximum over- 
pressure value. It can be seen that rise time ratios (and therefore 
absolute rise times since the overpressures are of the order unity) of 
from near zero to as large as about 0.065 are observed and the most 
frequently encountered value is about 0.01. Generally longer rise 
times are associated with rounded waveforms and much shorter rise 
times are associated with peaked waveforms. 

Since the rise time of the wave is a significant factor in the subjective 
response to sonic booms, an evaluation of this parameter was made 
for the range of flight conditions encountered. These data are 
presented in figure 10 as the variation of the bow wave rise times per 
unit overpressure as a function of altitude for measurements made 
within ±3 miles of the ground track. The data symbols represent 
the average values of the measurements and the vertical lines indicate 



26 


SONIC BOOM RESEARCH 



Figure 9. — Histogram showing variations of bow wave rise time to peak value 
of overpressure as measured from flights within ±3 nautical miles of ground 
track for the altitude range 70000+ feet. 



Figure 10. — Variation of boom wave rise time to peak value of overpressure 
as measured from SR-71 flights within ±3 nautical miles of ground track 
as a function of altitude. Data symbols represent average values for each 
flight along with maximum and minimum values observed. 
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the extreme spread of the data observed for each particular flight. 
It can be seen from the figure that there is considerable scatter, but, 
in general, the rise time per unit overpressure increases as the altitude 
of the aircraft increases. These results are similar to those obtained 
for the B-58 aircraft of reference 7. There is, therefore, the suggestion 
that the increased altitude results generally in longer rise times. 


p CONCLUDING REMARKS 

Sonic boom measurements have been obtained from 35 flights of the 
SR-71 vehicle at altitudes in excess of 70000 feet and Mach numbers 
to 3.0. No unusual phenomena were encountered for the extreme 
altitude and Mach number ranges of these tests, and the results fit 
generally into the established patterns of other available sonic boom 
flight data from F-104, B-58, and XB-70 aircraft. 
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Laboratory Sonic Boom Research and Prediction Techniques 

Harry W. Carlson 
Langley Research Center , NASA 


INTRODUCTION 

The present discussion is intended to serve as an addendum Jbo 
the report “Experimental and Analytic Research on Sonic Boom 
Generation at NASA” given at the NASA Conference on Sonic Boom 
Research, April 12, 1967, and published in NASA SP-147. Among 
the topics to be considered are refinement of the prediction methods 
which were discussed at some length in the previous report and the 
extension of wind-tunnel tests into the low hypersonic speed range. 
In addition, attention is directed to experimental data applicable to 
the multipole concepts for boom reduction under current consideration. 
Finally, some preliminary laboratory experiments dealing with shock 
propagation are described. 

SYMBOLS 

A effective cross-sectional area 
h airplane altitude or lateral distance from model 
l airplane or model length 

M Mach number 

A p incremental pressure attributed to flow field of airplane or 

model 

At time increment measured with respect to the arrival of the 
airplane bow shock 

Ax distance measured parallel to airplane or model longitudinal 
axis from bow-shock position to point on pressure signature 
a angle of attack 

d E elevon deflection 

Ah vertical separation distance between generating airplane and 
probe airplane 
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DISCUSSION 

Data from recently completed wind-tunnel tests of a 6-inch-long 
model of the XB-70 airplane are now being utilized in an effort to 
assess the applicability of present theoretical prediction methods to 
a complex configuration in the sensitive extreme near-field region of 
the flow. For the sample data shown in figure 1, it is seen that the 
basic theory which considers the area distribution for M=1.5 area 
rule cutting planes and the lift distribution for a flat wing surface 



Ax Ax 


Figure 1 .— XB-70 wind-tunnel data use in refinement of prediction techniques. 

does not completely predict the shape of the measured signature. A 
more complete theoretical treatment which considers the angle of 
attack attitude in derivation of the area and lift developments, and 
which accounts for interference effects by employment of a newly 
developed computing program for a warped lifting surface in con- 
junction with a canard or horizontal tail, is seen to result in an im- 
proved correlation. It should be pointed out, however, that for the 
large distances encountered in flight, where the signature is more 
fully developed, the more complete theory would result in an almost 
negligible change in the prediction. The impulse calculated for the 
more complete theory is only 2 percent greater than for the basic 
theory. The measured impulse is about 10 percent greater than the 
predicted value, which indicates a 5 percent difference in the extrap- 
olated value of the bow shock at large distances. In view of un- 
certainties concerning the nature of the model boundary layer, it 
is not known whether extrapolated experimental data or the theory 
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would offer the best estimate for the airplane. A completely laminar 
boundary was assumed for the model calculations. 

Analysis of signatures obtained in probe flights for the B-58 air- 
plane have brought to light an unexpectedly large effect of control 
surface deflection. As shown in figure 2, an improved correlation of 
experiment and theory results where the theory accounts for the 
estimated 7.5° deflection of the elevon (trailing edge up). Because 
of the increase in the impulse, this effect is expected to account for 


yi M = 1.6 

h □ 48,000 ft 
h = 1,340 ft 

Flight data & e «75° 



Figure 2. — B-58 probe flight data use in refinement of prediction techniques. 

about a 5-percent increase in ground overpressures. Elevon deflec- 
tion can have this degree of influence on the positive part of the sig- 
nature because more lift must be generated on the forward part of 
the airplane to compensate for the negative lifting force on the control 
surface. Thus it is suggested that control deflections be taken into 
account whenever high precision is required in theoretical predictions. 
However, it is believed that this particular case of airplane and flight 
condition represent a rather extreme situation. 

As a first step in the extension of the Langley wind-tunnel program 
to hypersonic speeds, tests were performed on a systematic series of 
body shapes at Mach numbers up to 4.63. An example of data 
from these tests is shown in figure 3. For the lower Mach number, 
the correlation of theory and experiment is quite good, at least for 
the positive part of the signature, and is representative of results ob- 
tained at Mach numbers below 3 for slender bodies. As might, have 
been expected, discrepancies between theory and experiment increase 
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h/J = 5 


o Experiment 
— Theory 


M = 2.96 


M = 4.63 




Figure 3.— Preliminary wind-tunnel tests at high supersonic speeds. 


with increasing Mach number and with increasing diameter-to-length 
ratio. Part of the discrepancy may be attributed to the increased 
model boundary-layer effect at the higher Mach numbers. The 
theory shown in the figure does not include the effect of boundary- 
layer thickness. However, it seems clear also that there is some 
failure of the Whitham method and that further theoretical studies 
are required. The body of revolution tests will be supplemented by 
an experimental investigation of simple lifting wings at speeds up to 
M— 4.63 planned for the near future. 

With the current interest in the use of multipoles as a means of 
achieving sonic boom reductions, it is thought appropriate to reexamine 
the results of wind-tunnel tests dealing with this subject. The pur- 
pose of the experiments, which were conducted about 10 years ago 
at the request of Dr. Busemann, was to examine the persistence of 
the directional properties of shocks created by quadrupoles. Measured 
signatures for a quadrupole device having a hatchet shape are shown 
in figure 4. As shown by the data, although there are quite large 
differences in pressure signature characteristics at one body length, 
the flow has become nearly axially symmetric at only eight body 
lengths. 

In another of the tests, a ducted quadrupole was studied. The 
shape was formed by distortion of a cylinder 2 inches in length and 
0.625 inch in diameter to form a 2:1 ellipse inlet section and a 2:1 
ellipse exit section with major axes 90° apart. Data for this con- 
figuration, shown in figure 5, indicate a significant difference at 




Figure 4.- — Wind-tunnel tests of a multipole device. Quadrupole, ikf «=2.0. 

eight body lengths between the signatures measured below and to the 
side of the configuration. Differences in impulse suggest that a 
10-percent difference in far-field overpressures should be expected. 
These results, then, indicate that significant directional properties 
are produced by quadrupoles only when the pressure-producing 
surfaces are displaced a significant distance from the body axis. 

Because these directional properties are achieved at the expense 
of large departures from conventional aircraft geometry and because 
the benefits may stem in part from the increase in effective body 
length defined by supersonic area rule concepts, the results are 
compared in figure 6 with an estimated signature for a parabolic 
body having the same effective length and the same maximum cross- 
sectional area. The area distribution for the quadrupole body 
(determined by M— 2.0 cutting planes) includes the effect of an 
estimated laminar boundary layer and has the stream-tube capture 
area removed. This analysis indicates that benefits of the quadrupole 
are slight when compared with the simpler expedient of lengthening 
the body. 

This survey of miscellaneous research items will conclude with a 
description of some recent laboratory studies of shock-wave propaga- 
tion conducted by Dr. Raymond Barger of the Langley Research 
Center. The peculiar behavior of a shock wave in the region near 
a focus was studied by means of the experimental setup shown in 
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Figure 5. — Wind-tunnel tests of a multipole device. Ducted quadrupole, M — 2.0. 

Experiment 
— Theory 
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Figure 6.— Ducted quadrupole compared with, a parabolic body. M— 2.0, 

hfl—B. 

the sketch of figure 7. An N-wave signature triggered by a spark 
was reflected from parabolic mirrors so as to converge the shock 
front at a sharply defined focus. Before the focus, as shown in the 
signature at the left, there is a normal N-wave with a leading positive 
pressure jump. At the focus the measurements show not only a 
large amplification of the shock intensity, but also a change in signa- 
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ture shape identified as a phase shift of 90°. Beyond the focus there 
is the interesting condition of a leading expansion portion of the 
signature and a signature which corresponds to an 180° phase shift. 
These results are in agreement, at least qualitatively/ with the pre- 
dictions given by the laws of geometric acoustics. 

A second laboratory experiment deals with focusing a shock front 
in the presence of a temperature gradient. In this experiment a 
small projectile was fired down an instrumented range in which was 
placed a chamber with a heated wall to produce a temperature 
gradient normal to the path of the projectile. The primary instru- 
mentation consisted of pressure transducers and apparatus for 
schlieren photography. Some of the photographic results are shown 
in figure 8. The path of the projectile traveling at a Mach number 
of about 1.1 is indicated at the top of each photograph. For the 
photograph at the left, there was no temperature gradient. At the 
center where a moderate temperature gradient was introduced, a 
noticeable curvature of the shock front is evident. The temperature 
is higher near the bottom of the picture, corresponding to the real 
atmosphere situation. At the right, the gradient is large enough so 
that the speed of sound near the center of the field is equal to the 
speed of the body. There, as would be expected, the shock front 
becomes normal to the flight path and loses its identity. 

With further development, these laboratory techniques could serve 
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Figure 8. — Laboratory study of shock field behavior in a temperature gradient 


a valuable purpose in the study of various aspects of propagation 
theories which have been developed and others which are now being 
developed. 


Configuration Design for Specified Pressure Signature Characteristics 

F. Edward McLean 
Langley Research Center , NASA 


INTRODUCTION 

Near-field design concepts, which have evolved from NASA studies 
(refs. 1-3), appear to offer a promising means for alleviating the sonic 
boom problem associated with the supersonic transport (SST) . 
These studies have shown that, for certain flight conditions, the 
effective shape of a large airplane can be tailored to provide ground 
pressure signatures which may be more desirable than the typical 
far-field N-wave. This paper presents design requirements for the 
plateau pressure signature and the signature with finite rise time, 
two signature shapes which exhibit interesting possibilities for sonic 
boom reduction. Consideration will be given to problems associated 
with the practical attainment of these pressure signatures within 
the normal operating regime of a supersonic transport. 

SYMBOLS 

A e effective area distribution of airplane or model equivalent 
body, including effect of lift as well as volume 
x distance measured along longitudinal axis from airplane or 
model nose 

A p incremental pressure attributed to flow field of airplane or 
model 

p reference pressure 

M Mach number 

h altitude 

W weight 

At incremental time from onset of pressure disturbance 

W g gross weight 

tr time required for maximum pressure rise within pressure 

disturbance (see fig. 5) 

Ax distance measured parallel to longitudinal axis of airplane 
from point on linear theory characteristics from nose of 
airplane to point on pressure signature 
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NEAR-FIELD SIGNATURES UNDER CONSIDERATION 

During the course of research on the sonic boom problem, NASA 
has investigated the pressure rise characteristics of a number of 
airplane configurations and equivalent body shapes (refs. 1-5). 
Considerations of near-field design concepts have been focused on 
the generating shapes and signature types illustrated in figure 1. 
In the left side of the figure is the plateau pressure disturbance. 
This signature shape, which is generated by an effective area devel- 
opment that varies as the three-halves power of distance along the 
axis, is characterized by a constant pressure region behind the bow 
shock. For airplane operating conditions where it is possible to 
generate this type of signature on the ground, the total pressure rise 
in the disturbance would be somewhat less than that of the correspond- 
ing far-field N-wave (ref. 1). The plateau pressure signature has 
found particular application in the reduction of sonic boom over- 
pressures during the critical transonic climb portion of a supersonic 
transport mission. 

A pressure signature with finite rise time is illustrated in the right 
side of figure 1. This signature is generated by an effective area 
distribution, which varies as the five-halves power of distance along 
the axis, to about the midpoint of the body length. Aft of the body 
midpoint, the area development is a much more complicated vari- 
ation with distance along the axis. The interesting feature of the pres- 
sure signature with finite rise time is that there is no instantaneous 

PLATEAU FINITE RI SE TIME 



Figure 1 . — Near-field pressure signatures under construction. 
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shock pressure jump, but rather a gradual pressure rise at both the bow 
and tail portions of the signature. If the buildup in pressure can 
be extended over an appreciable rise time, on the order of 10 to 15 
milliseconds, there are indications that the noise associated with the 
sonic boom disturbance can be virtually eliminated (refs. 3 and 4). 
However, because of the extreme design requirements for this type 
of signature, the practicality of its application is still under investi- 
gation. 

Note that the experimental pressure signatures of figure 1 were 
measured in the flow fields of bodies of revolution with the required 
area developments. The work of Hayes (ref. 6) permits relation- 
ships to be drawn between the pressure rise characteristics of equiv- 
alent bodies of revolution and lifting airplane configurations. These 
relationships have been widely used in sonic boom investigations. 


USE OF PLATEAU SIGNATURE IN DOMESTIC SST STUDY 


The sonic boom problem raises some serious questions as to the 
acceptability of domestic or overland flight of the supersonic transport. 
Because of the economic importance of such a domestic operation, 
however, the Langley Research Center of NASA has recently con- 
ducted a study of the feasibility of a relatively low-boom, domestic 
supersonic transport. In the development of the domestic SST 
study configuration, the most advanced aerodynamic technology 
and concepts were utilized. Near-field design procedures were used 
to take advantage of the sonic boom overpressure reductions which 
have been indicated for the plateau pressure signature. 

The manner in which near-field design concepts were applied to 



Figure 2. — Use of plateau signature in domestic SST study. 
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the domestic SST study configuration is illustrated in figure 2. Be- 
cause airplane shaping for a specified signature can be accomplished 
only at one flight condition, a near-field design Mach number (M— 
1.8) was chosen, intermediate between the critical transonic (M== 1.2) 
and cruise (M—2.7) Mach numbers. At the design Mach number 
(M=1.8), the volume and lift elements of the airplane were carefully 
tailored to provide a three-halves power effective area distribution 
for the assumed altitude of 60 000 feet and airplane weight of 338 000 
pounds. The plateau pressure signature which would be expected 
for these flight conditions is shown in the right side of figure 2. 

Near-field sonic boom overpressures which would be generated by 
the domestic SST study configuration are compared with correspond- 
ing far-field or N-wave overpressures in figure 3. For most of the 
flight Mach number regime, the previously mentioned application of 
near-field design concepts leads to overpressure levels considerably 
below far-field values. The sonic boom characteristics of the study 
configuration vary little with Mach number and are at a level of about 
1.4 psf. Note that the benefits from the plateau pressure design are 
felt on either side of the near-field design point and persist up to the 
cruise Mach number of 2.7. At cruise conditions, however, the sonic 
boom approaches far-field levels. 

Although only minimum overland flights of proposed international 
supersonic transports are currently anticipated, some consideration is 
being given to off loading fuel from these airplanes for the shorter 
domestic ranges of about 2500 nautical miles. It would be interesting, 
therefore, to compare the sonic boom characteristics of these off-loaded 


OVERPRESSURE ASSUMED CONDITIONS 
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Figure 3.- — Sonic boom characteristic of domestic SST study configuration 
compared with far-field theory. 
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Figure 4. — Sonic boom of domestic SST study configurations compared with 
sonic boom of proposed SST’s. 

international SST’s with the sonic boom characteristics of the NASA 
domestic SST study configuration. This comparison is made in 
figure 4. The band of values at the top of the figure represents the 
sonic boom characteristics of off-loaded versions of proposed inter- 
national SST’s as a function of mission range. For these airplanes, 
the overpressures approach 2.0 psf at the beginning and end of the 
mission and reach a minimum value of about 1.7 psf at the end of 
cruise. Note that the relatively small (W g ~ 300 000 pounds, l = 
190 feet), European-proposed international SST has slightly higher 
overpressure levels than the large (W g = 576 000 pounds, l— 306 feet), 
American-proposed SST. 

The curves at the bottom of figure 4 indicate the sonic boom over- 
pressure levels for the NASA domestic SST study configuration for two 
operating procedures. The solid curve corresponds to a one-step cruise 
operation at a Mach number of 2.7, and the dashed curve represents a 
two-step cruise procedure with cruise legs at M=1.6 and M— 2.7. 
For one-step cruise operation, the maximum overpressure of the 
domestic study configuration is 1.4 psf, and the minimum overpressure 
is approximately 1.2 psf. The two-step cruise procedure leads to 
approximately a 0.2-psf reduction in sonic boom overpressure at some 
loss in range capability. Note that the overpressure levels for the 
NASA domestic SST study configuration are some 0.6 to 0.8 psf less 
than the values shown for proposed international SST’s which have 
been adapted to domestic operation. This suggests the possible 
desirability of developing a supersonic transport configuration 
specifically for overland use. 
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CONSIDERATIONS OF PRESSURE SIGNATURES WITH FINITE RISE TIME 


One of the most undesirable characteristics of the sonic boom 
disturbance is the instantaneous pressure jump which is normally 
associated with the bow and tail wave of the disturbance. The elimi- 
nation or modification of these sharp pressure jumps would appear to 
be desirable. In reference 3, consideration was given to the design 
criteria necessary to spread the bow pressure rise over a finite time, 
and thus eliminate the instantaneous bow pressure jump. For normal 
operating conditions of the supersonic transport, this finite rise time 
consideration led to rather extreme design requirements of airplane 
length and slenderness. If the tail pressure jump is similarly modified 
to provide a finite rise time, more extreme design requirements would 
obviously result. 

Figure 5 indicates the rather extreme design requirements necessary 
to produce a ground pressure disturbance with finite rise times at 
both the bow and tail of the disturbance. The figure is drawn for 
typical cruise conditions of a supersonic transport, M— 2.7 and n— 
61 000 feet. To provide desirable rise times of 10 to 15 milliseconds 
at the ground, extreme airplane effective lengths are required. For 
anticipated supersonic transport weights, of the order of 600 000 
pounds, an effective length of approximately 1000 feet would be 
required for the assumed operating conditions. In addition, the air- 
plane effective area distribution would have to follow the development 
shown in figure 1 for a signature with finite rise time. 
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Figtjke 5. — Design requirements for signature with finite rise time. 
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The geometric problem associated with application of the extreme 
design requirements of figure 5 to a typical supersonic transport 
configuration is illustrated in figure 6. To produce a pressure signa- 
ture at the ground with a desired 10- to 15-millisecond rise time, the 
effective area development of the basic configuration would be modified 
to the area distribution indicated by the dashed curve. This modifica- 
tion would require a substantial increase in effective area and length 
extensions of 400 feet ahead of and 300 feet behind t,he basic configura- 
tion. Such a geometric modification appears to be impractical. 


M = 2.7, h= 61,000 ft 




BASIC CONFIGURATION 



Figure 6. — SST modification for finite rise time. 


Several schemes, such as the use of electrostatic effects or laser 
beams, have been proposed as possible ways of increasing the effective 
area and length of an airplane without resorting to geometrical 
modifications. Questions have risen as to how these schemes will 
provide the required control over area development or what the power 
costs would be. Recently, Carlson of NASA instituted a research 
program to provide answers to the latter question. Assuming that 
some unidentified force field or heat source is available, Carlson will 
determine the power required to create around the existing geometrical 
configuration a “phantom” body which will produce a pressure 
signature with finite rise time. 

The model for Carlson’s investigation is illustrated in figure 7. 
An entering stream tube of air, A 0 , will be acted on by force or heat 
to produce the phantom boundary illustrated at the top of the figure. 
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This boundary will alter the effective area distribution as indicated 
in the middle of the figure, with the cross-hatched region corresponding 
to the airplane and the solid curve representing the phantom body. 
Power required will vary with A 0 as indicated at the bottom of the 
figure. Preliminary results indicate that twice the maximum engine 
power available would be required to produce the modified signature 
illustrated in figure 6. 

The extreme geometric modifications or power required to produce 
pressure signatures with finite rise times would perhaps prohibit their 
attainment under the normal flight conditions of a supersonic trans- 
port, Note, however, that these extreme design requirements are 
based on theoretical considerations alone. Random atmospheric 
distortions, which are not considered in the theory, have produced 
unexpected rise times in measured signatures. With the assistance 
of these atmospheric effects, the design requirements for this desirable 
type of signature might be brought within practical limits. Further 
investigation appears to be warranted. 

x — CONCLUDING REMARKS 

Near-field design concepts have been discussed in relation to the 
\sonic boom problem of the supersonic transport. The practical 
application of the design requirements for a plateau pressure signature 
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to reduce the sonic boom overpressures of a domestic SST configura- 
tion has been demonstrated. Theoretical design requirements for 
pressure signatures with finite rise time have been presented for 
typical operating conditions of a supersonic transport. The problems 
associated with the practical attainment of these latter signatures 
have been considered. 
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Evaluation of Certain Minimum Boom Concepts 

Harry L. Runyan and Herbert R. Henderson 
Langley Research Center , NASA 


INTRODUCTION 

This paper is concerned with two study programs relating to sonic 
boom alleviation. In the first, an experimental program was con- 
ducted to study the growth of the pressure from a basic two-dimen- 
sional configuration to its three-dimensional flow field. In the second 
program, several planforms were examined theoretically, and an 
attempt was made to formulate a sonic boom efficiency factor, the 
use of which would provide for a consistent comparison of candidate 
configurations. 

EXPERIMENTAL PROGRAM 
Concept and Model 

The fundamental idea behind the theoretical work of Hayes 
(ref. 1) and Whitham (ref. 2) in the calculation of the sonic boom is 
the replacement of a three-dimensional aircraft with an equivalent 
body of revolution. This concept then presumes the accurate 
establishment of three-dimensional flow in the far field of the aircraft, 
and has been shown to be a very good approximation for normal air- 
plane shapes through comparison with both wind-tunnel studies and 
actual flight experience. 

The adequacy of the method, however, has not been thoroughly 
tested for shapes that do not appear essentially as bodies of revolu- 
tion. As a matter of fact, many concepts for shockless or low-boom 
airplanes are being put forward on the basis of two-dimensional 
reasoning. One such configuration, developed on the basis of two- 
dimensional flow, is shown in cross section on figure 1. The idea 
here was to create a lifting pressure on the lower surface for a small 
distance involving a shock wave, then expand the flow to the original 
free stream direction, creating a shock canceling expansion segment 
which would be very close to the shock wave and to the airplane. 
This shock cancellation procedure would be repeated downstream 
for several steps, thus leading to a number of lifting segments, followed 
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•PRESSURE DISTRIBUTION UNDER AIRFOIL 
(TWO-DIMENSIONAL) 


Figure 1 . — Conceptual model for shock cancellation technique based on two- 

dimensional flow. 


by shock canceling regions. The top of the aircraft would be flat 
and alined to the original stream flow direction in an effort to reduce 
any unwanted drag-producing shocks on the upper surface. Also, 
it was presumed that the engines could be designed so that they would 
exhaust at the same static pressure as the free-stream static pressure, 
thus eliminating closure shocks. 

For the model concept illustrated in figure 1, the pressure (on a 
two-dimensional basis) was calculated just below the model. These 
stepped pressures as shown on figure 1 were used as input in an un- 
published procedure using the method of characteristics for calcu- 
lating the pressure on the ground through a standard stratified atmos- 
phere. The calculated result for a 200 000-pound airplane at M— 3 
yielded an overpressure of one-fourth psf on the ground. This looked 
encouraging, so that a test model was subsequently constructed and 
tested in the Langley Unitary Plan Wind Tunnel at M= 2.7. The 
model of span 4 inches had four steps, each step having an angle of 
approximately 5° and being ^ inch long followed by a 3^-inch 
stream wise portion. The rear of the model was extended with a 
constant rectangular cross section which simulated the engine exhaust 
requirement of exhaust at free stream static pressure, and this constant 
cross section provided for the sting support. 

Results 

Pressure distributions below the model were obtained for several 
positions, and the flow field visualization was obtained by schlieren 
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photography. In figure 2 are shown the profile and plan view of the 
schlieren photograph for M— 2.7. On the top of the figure the plan 
view is shown, and it is interesting to note that no discernible shock 
waves are seen to originate from the corners of the leading edge. 



Figure 2. — Schlieren photograph of shock cancellation model test, Af— 2.7. 


The lower photograph shows rather clearly the shock and expansion 
pattern. The dark lines indicate shock waves, and the light areas 
indicate the expansion zones. The vertical white lines are caused 
by the schlieren window supports which are external to the flow. By 
following the shock waves from the model, one can see the coalescence 
of the shock waves, and at the lower right the four shock waves gener- 
ated at the model have been replaced by a rather strong leading edge 
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shock and two apparent weaker shocks. The remaining shocks in 
the photograph are attributed to the sting support. 

Two pressure distributions were measured below the model at 
distances of 5 and 20 inches, and the results are shown in figure 3. 



“5 0 5 10 15 20 

X, in. 

Figure 3. — -Theoretical and experimental pressure distribution for shock can- 
cellation model at 5- and 20-inch distance below model. 

{Because the model was 4 inches in length, these measurement stations 
correspond to 1.25 and 5 body lengths.) The upper curve gives 
the results for the 5-inch distance. At least for the near position, it 
was thought that a pressure distribution similar to that shown below 
the model on figure 1 (two-dimensional type) would have been ob- 
tained at 1.25 body length position. Instead, the formation of a 
typical N-wave is already apparent, with the shocks from the steps 
showing only slight perturbations on the flow. The theory of 
Carlson (ref. 3) was used as a basis for the calculations shown. 
Considering that the theory was never intended to be used so close 
to the body and the body was not axisymmetric, the agreement is 
thought to be rather good. For the 5-body length position, the 
. influence of the steps has been almost canceled, and relatively good 
agreement is found between theory and experiment, and a typical 
N-wave has been formed. For both cases, the influence of an adapter 
located between the model and the sting support was found to be 
rather strong, and calculated results with the effect of the adapter 
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included were made, and these should be compared to the experi- 
mental results. 

It was mentioned previously that the two-dimensional pressure 
distribution just below the model was used as an input to a computer 
program to calculate the pressure on the ground through a standard 
but stratified atmosphere, and the resulting ground A p was one- 
fourth psf. After the tests, the experimental pressure distribution 
at a distance of 20 inches below the model was used as input to the 
same computer program, and a ground pressure of 2.5 psf was calcu- 
lated. This again illustrates the danger of deriving shockless or 
minimum boom configurations based on two-dimensional flow fields. 

Discussion 

A concept developed on two-dimensional reasoning has been 
studied, both experimentally and theoretically, and the major con- 
clusion reached was the rapidity with which the three-dimensional 
flow was established and the relatively good agreement with a theory 
which has been developed for the so-called midfield or far-field cases. 
The results should serve as a reminder that minimum or zero boom 
concepts based on two-dimensional reasoning can be very misleading 
and that the development of three-dimensional flow is predominant 
and must be taken into account. 


SONIC BOOM EFFICIENCY FACTOR 

There are a large number of possible configurations for supersonic 
aircraft, some of which may provide low sonic booms. Factors, 
other than low boom, must be considered in evaluating a supersonic 
transport concept, and one of the more important is the lift-drag 
ratio ( LjD ) } because the range and efficiency of the aircraft are highly 
dependent on this ratio. Other factors such as structural weight, 
engine-airframe integration, stability and control, and low-speed 
landing characteristics are of course important, but these will not be 
included in an efficiency factor. In studying a number of config- 
urations, it seems worthwhile to form a ratio so that the various 
configurations could be systematically compared. Because it is 
desired to maximize LjD and minimize the sonic boom overpressure 
A p, the ratio 


would form a suitable ratio for a comparative analysis, where 



refers to the maximum overpressure on the ground. 
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Utilizing this concept, three configurations were studied by use of 
the theory of reference 3. The calculated Ap was corrected for 
altitude and standard atmosphere effects to the A p on the ground 
by use of the factor given in reference 3. The three configurations 
studied had rectangular, delta, and arrow planforms. 

Calculations were performed for a constant-weight airplane of 
400 000 pounds which might be suitable for transcontinental travel. 
Only the lift and drag attributed to a flat plate wing were calculated, 
and no skin friction drag was included in the drag estimate. For 
simplicity, it was also presumed that a fuselage could be designed 
which created no wave drag; e.g., by a Busemann biplane rotated 
90° about its longitudinal axis. This would conceivably provide 
zero wave drag and thus the fuselage contribution would be small. 
Results for the three planforms are shown in table I. The dimension 
l is the streamwise length of the wing in feet, and b is the span in 
feet. The overpressure on the ground is A p in pounds per square 
foot. The final column is the efficiency factor 

Wh 

It is desired to obtain a configuration that has the highest efficiency 
factor, while at the same time having a low A p. Results are given 
for three angles of attack. The planform having the highest efficiency 
factor of those studied is the 70° arrow wing of aspect ratio 20, at an 
angle of attack of 1°, which has a factor of 186 000 as well as having 
the lowest A p of 0.67 psf. This configuration is shown in figure 4, 
and it is obvious that structural considerations would compromise 
this configuration. A surprising result of this study is the low 
efficiency of the delta wing, when compared to the arrow or even the 
rectangular wing. The efficiency factor for the delta wing is the 
lowest for all cases, and the ground A p is the highest for all cases 
studied, when comparing the delta wing to the other two configura- 
tions. The rectangular wing of aspect ratio 1 rates next to the 70° 
arrow wing (AR— 20) in the cases studied. Structurally, a flat 
rectangular wing may be an efficient structure, and if the weight 
factor were brought into the picture, it could conceivably be a prac- 
tical configuration. It is felt that the efficiency factor concept should 
be broadened to include at least some structural weight parameter 
to form a more general efficiency factor. 



Table L — Sonic Boom Efficiency Factors for Several Configurations (M=S, Airplane Weight = IfiO 000 Pounds , and 

Altitude = 70 000 FeeJ) 


Wing configuration 

Aspect ratio 

C L 

Area 

Dimensions 

A p 


l 

b | 

«=1° 

R.p.e.tfin gnlar 

1. 0 

0. 021 

32 400 

180 

180 

1. 27 

95 400 

RAfttAllgulftf • • 

2. 5 

.023 

30 200 

110 

275 

1. 58 

76 800 

Rectangular 

5. 0 

. 024 

28 100 

75 

375 

1. 78 

68 100 

70° arrow - • _ _ _ _ 

20. 0 

.044 

15 700 

770 

560 

. 67 

186 000 

70° arrow-.-- 

2. 7 

. 210 

3 300 

130 

94 

1. 78 

69 800 

70° delta--- — 

1.5 

. 331 

2 050 

75 

54 

2. 18 

56 900 


a = 3° 


Re ct angular 

1. 0 

.063 

11 000 

105 

105 

1. 57 

25 800 

Rectangular 

2. 5 

.070 

9 000 

60 

150 

1. 85 

21 800 

Rectangular 

5. 0 

.072 

10 100 

45 

225 

2. 03 

19 900 

70° arrow - - 

20. 0 

. 131 

5 250 

445 

332 

1. 05 

39 600 

70° arrow 

2. 7 

.630 

1 100 

75 

54 

2. 05 

20 200 

70° delta 
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Figure 4. — Illustration of arrow wing configuration used in sonic boom efficiency 

factor study. 


CONCLUDING REMARKS 

The first part of this paper was concerned with an experiment a 
and theoretical study of a rectangular configuration, and a study 
was made of the development of three-dimensional flow. In spite 
of the fact that the configuration physically did not resemble a body 
of revolution, as the theory assumes, relatively good correlation was 
obtained between theory and experiment at a distance of approxi- 
mately only 1 body length. It emphasized the point that low-boom 
configurations based on two-dimensional reasoning may not perform 
according to expectancy, and that any configurations based on two- 
dimensional concepts must be studied in their real three-dimensional 
regime. 

The second part of the paper suggests the use of a sonic boom 
efficiency factor L/D(Ap/p ) 0 . Three planforms were examined by 
the use of this factor, and it was found that the arrow wing of very 
high aspect ratio was the best configuration, having low boom and 
high L/D, with the rectangular wing second and, surprisingly, the 
delta wing having the least efficiency from the sonic boom standpoint. 

It is felt that the structural aspects of these configurations should 
be studied, and a more complete efficiency factor formulated that 
would include structural weight. \ 
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Current Research in Sonic Boom 

Lynn W. H unton 

Ames Research Center , NASA 


INTRODUCTION 

The research program in sonic boom for this past year at the NASA 
Ames Research Center consists of three general study areas. The 
program covered the development and verification of new prediction 
methods (new computer methods and near-field experimental data 
method), a study of the effects of Mach number to 5.5, and a boom 
minimization study. The purpose of this paper is to review briefly 
the status of this research effort. 

PREDICTION METHODS 

There is a need for analytical methods that are capable of accurately 
predicting the near flow field pressures about arbitrary aircraft con- 
figurations for design optimization in connection with the sonic boom 
as well as the performance of the vehicle. Available theories are 
either too cumbersome mathematically for ease of handling on the 
computer (e.g., method of characteristics) or are limited in range of 
application (e.g., Whitham, ref. 1) because of some mathematical 
simplification such as the linear theory concept. As a result, prob- 
lems of analysis can occur in flow regions very near the aircraft in 
connection with the design and arrangement of configuration com- 
ponents. To help remedy this situation, two new flow field prediction 
methods are currently under development. One of these, utilizing 
the latest in. computer display technology, is described by Harvard 
Lomax in the paper entitled “Preliminary Investigation of Flow Field 
Analysis on Digital Computers With Graphic Display” in this 
volume. The second method, which will be developed under contract, 
consists of a modified version of an existing supersonic wing-body load 
calculation program wherein the configuration is represented with a 
spatial distribution of singularities over the entire wing and body 
surfaces. Because the classical method of Whitham, in contrast, 
uses an equivalent linear distribution of singularities located only 
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along the body axis, it is expected that somewhat improved solutions 
in the near field will be obtainable from this new calculation procedure. 

The final prediction technique relates to a method that w T as de- 
veloped at Ames for extrapolating experimental near-field pressure 
signature data obtained from wind-tunnel models to predict the sonic 
boom characteristics of airplanes flying at any altitude (ref. 2) . This 
new procedure overcomes some of the limitations in previous experi- 
mental and theoretical techniques by permitting studies of relatively 
large models in the wind tunnel. Some verification of the procedure 
has been made with the XB-70 and the X-15 airplanes for Mach 
numbers of 1.8 and 5.5, respectively, as shown in figure 1. Using 
scale models of these airplanes, near-field pressures were measured 
in the wind tunnel in a plane one body length below each model. 
This near-field pressure signature was then used to predict the dashed 
curve for the XB-70 at a distance ratio of 4.5 to compare with avail- 
able wind-tunnel data measured at that same distance. It can be 
seen that the comparison is good despite the definite near-field char- 
acter of the curve. On the basis of the correlation shown for two 
distance ratios it is apparent for this case of a lifting airplane con- 
figuration that a reliable F-function can be partially defined from 
pressure measurements made as close as one body length from the 
configuration. As described in reference 2, this incomplete F-f unc- 
tion is then believed valid for the calculation of the pressure signature 
at any distance ratio larger than that used in defining the F-function. 
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Figure 1.— Evaluation of experimental procedure. 
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Consequently, extrapolations of the overpressures are shown carried 
out to flight distance ratios for the actual airplanes flying at altitudes 
of 50 000 and 90 000 feet, respectively, to compare with available 
flight data obtained from NASA Flight Research Center. Again, 
the comparison shows very good correlation except behind the tail 
shock where there usually is a problem of simulating the actual air- 
plane with a sting mounted model. It should also be mentioned that 
the experimental flight data as presented represent an average of the 
pressure readings from several pressure cells. 

EFFECTS OF MACH NUMBER 

The current interest in exploring the potential of the hypersonic 
air-breathing vehicle for purposes of global transportation raises 
several questions regarding the sonic boom problem. First, will the 
good correlation between experiment and the Whitham theory found 
at moderate supersonic Mach numbers persist at hypersonic Mach 
numbers? Second, how do the levels of overpressure compare between 
a supersonic and a hypersonic cruise mission? Third, at the higher 
Mach numbers, is the overpressure any less sensitive to changes in 
configuration geometry? These are some of the areas that have been 
investigated, as described in reference 3, and for this review examples 
of each will be considered briefly. 

In figure 2 is shown a comparison of overpressure signatures 
from experiment and from the Whitham theory for a 7.5° half- 
angle cone cylinder for Mach numbers from 2 to 5.5. Experimental 



Figure 2. — Comparison of experiment with theory. 
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results were obtained at a distance of one body length and these data 
have been extrapolated to a distance ratio of 100 for the comparison. 
This distance ratio was chosen to ensure that the near-field restric- 
tions on the application of the Whitham theory would not be a 
factor in the comparison. For this simple case of a body of revolution, 
for which the theory should be particularly suited, it can be seen that 
the correlation of experiment and theory for the bow shock is fairly 
good only to a Mach number of 3. Beyond Mach 3 the theory leads 
to a definite underprediction of the strength of the bow shock. This 
result is not too surprising in view of the assumptions used in the 
development of the Whitham theory which place a definite Mach 
number limitation on the theory (see ref. I), Also to be noted is a 
significant increase in the signature length and an unbalancing of 
the positive and negative impulse areas. While these trends were 
predicted by the linear theory, the magnitude of these effects fell 
somewhat short of the measured characteristics. 

In figure 3 consideration is given to a comparison of the over- 
pressure levels for a supersonic and a hypersonic transport (HST) 
mission using geometrically similar blended wing-body configurations 
for the two cases. The SST is a smaller and lighter airplane with a 
gross weight of 500 000 pounds compared to 600 000 pounds for the 
HST. On the left of the figure is shown the assumed flight path 
schedule of altitude versus Mach number for a supersonic cruise at 



Figure 3. — Maximum overpressure for mission profiles. 



CURRENT RESEARCH IN SONIC BOOM 


61 


Mach 2.7 and a hypersonic cruise at Mach 6. Starting at the same 
Mach number and altitude, each mission includes a climb and accelera- 
tion leg, a Breguet cruise, and a descent leg. On the right is shown 
the maximum overpressure as a function of Mach number that has 
been derived by extrapolation of wind-tunnel data obtained for the 
model at a distance ratio of one body length. At the start of the 
mission with both airplanes at the same Mach number and altitude, 
the overpressures are separated initially by differences attributed to 
size and gross weight of the airplanes. As can be seen, the hypersonic 
mission imposes rather severe overpressures during the climb, whereas 
in cruise it reduces the overpressure to about half that for the super- 
sonic mission. If, with further technological developments in propul- 
sion, the climb schedule for the hypersonic mission can be altered as 
shown by the dotted line, then the overpressure could be improved 
considerably. 

Effects of configuration geometry on the maximum peak over- 
pressures have been investigated at Mach numbers of 2, 3, 4, and 5.5 
for three configurations representative of different concepts for air- 
fare .thing hydrogen-fueled hypersonic transport aircraft. The test 
models, illustrated in figure 4, were complete with empennage and 
simulated engine inlets and were 4 inches in length. Based on near- 
field measurements obtained at a distance of one body length, com- 
parisons of the overpressure characteristics as shown in figure 5 were 



Figure 4.— Test models. 
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Figure 5 .— Configuration effect: altitude, 50000 feet; weight, 600000 pounds; 
constant volume derived from wind tunnel h[l~ 1. 

derived for the airplanes having constant volume and flying at a 
constant altitude of 50000 feet and a weight of 600000 pounds. 
Because of the volume constraint, the lengths of the assumed air- 
planes were not constant as shown by the silhouettes in the figure. 
It may be seen that the all-body configuration gave the highest over- 
pressured level which worsened with increase in Mach number. 
While no attempt was made to optimize any of the configurations 
for boom purposes, it is believed that the canard on the all-body 
configuration in combination with the shorter length of this configura- 
tion may have been largely responsible for the measured differences in 
overpressure that are shown. From these limited results it is clear 
that the overpressure characteristics at the higher Mach numbers are 
sensitive to configuration differences at the moderate altitude 
considered. 

BOOM MINIMIZATION 

Based on current knowledge of the problem, the boom minimization 
program divides logically into two principal phases as shown below: 
Wing configuration study 
Planform 
Dihedral 
Camber 

Two-dimensional 
Three-dimensional 
Wing-body integration 
Conventional bodies 
Shock canceling bodies 

Because it is generally conceded that the thickness effects on sonic 
boom, in principle, can be handled by shock interference methods, 
the initial effort here is being directed at the lifting aspects of the 
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problem for which there are no known solutions. The approach 
taken, as the first step, is to examine wing configurations that are to 
some extent unconventional in order to search out possible new design 
concepts to boom minimization, that may have been overlooked in 
earlier studies. Currently being explored are wing planform, di- 
hedral, and camber. The investigation of camber includes studies 
of airfoil thickness distribution which is being examined in two-dimen- 
sional as well as three-dimensional tests. 

The second phase of this program will involve an integration of 
the more promising lifting arrangements with bodies having optimum 
thickness. Included here will be conventional bodies with area 
ruling as required to yield favorable midfield signatures as well as 
shock canceling bodies that utilize perhaps a three-dimensional 
application of the shock interference concept of Busemann. 

Shown in figure 6 is a display of the initial set of wing configurations 




being examined. Each has been selected to feature some particular 
arrangement of the wing shock system, such as intersecting, curved, 
or vertically displaced shocks, to explore our understanding of shock 
system combinations. Consequently, this selection of wings is 
intended only to serve as a baseline for future work. In addition 
to these symmetrical configurations, the arrow wing is also being 
tested with a cambered airfoil section and the delta wing is being 
studied with wing dihedral of ±10°. The geometric parameters 
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held constant in this series included: the body with a parabolic arc 
nose and a length of 7 inches, the total wetted area, the wing span, 
the exposed aspect ratio, and the wing thickness consisting of a double- 
wedge section with a maximum thickness ratio of 5 percent. The 
one exception to these stated conditions was the circular wing which, 
to satisfy the wetted area requirement, had a smaller span and hence 
a lower aspect ratio. 

A preliminary summary of the overpressure characteristics measured 
for these various wings in the presence of the fixed body is given in 
figure 7. Maximum peak overpressures resulting from the wing as a 
function of lift coefficient are presented for 12 wings for a Mach 
number of 1.4 and a distance of four body lengths. The assigned 
numbers, in addition to identifying the curves, can be seen to indicate 
the general progression of overpressure from the highest to the lowest 
level at the approximate lift coefficient for transition flight of 0.15. 
Also shown is a table on the right which is intended to summarize in 
somewhat gross terms the measure of success with which the charac- 
teristics of these several wings could be predicted by theory. 

The principal conclusions to be drawn from this data summary are 
as follows: 

(1) A relatively large spread in overpressures was obtained for this 
series of wings ranging from the unswept trapezoidal down to the 
highly sweptback arrow and the limits of this spread were predicted 
by theory. 



PREDICTION 

GOOD 


FAIR 

GOOD 


FAIR 

GOOD 

POOR 


Figure 7. — Summary of wing peak overpressures : M = 1 .4 ; h[l = 4. 
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(2) Camber for the arrow wing was advantageous at the higher lift 
coefficients and warrants further detailed study. This effect was not 
predicted by theory. 

(3) The variation in lift effectiveness on the overpressure as 
evidenced by the differences in the slopes of the curves is quite large 
and preliminary estimates of these effects were generally successful. 
As already mentioned, some problems arose in the case of the cambered 
arrow wing. 

(4) Wing dihedral on the delta configuration showed a surprisingly 
large effect on the overpressure which has been predicted with only 
fair success. 

It is of interest to examine the dihedral effect in a little more detail. 
Using curve 5 for the plane delta as a reference, it can he seen that 
curve 4, for 10° anhedral, indicates a higher overpressure, whereas 
curve 10, for 10° dihedral, shows a significant reduction in over- 
pressure. Also note, however, that these separated curves are nearly 
parallel which means that the effect of dihedral measured here is 
attributed primarily to a volume effect, presumably from the vertical 
displacement of the wing volume, and not to any favorable lift effect 
as had been questioned. 

To help put these dihedral wing results into a little better per- 
spective, figure 8 has been prepared which shows some additional 
characteristics measured at a Mach number of 1.7 and a lift coefficient 
of 0.2. First note on the left of the figure the small silhouette of the 
test configuration which is scaled to show the relative size and distance 
of the model from the probe station located at the ground level as 
tested in the wind tunnel that is equivalent to a distance of 4.5 body 
lengths. Also shown on the left are the peak overpressures measured 

h 11 = 4.5 CENTER-LINE VALUES 




Figure 8.— Dihedral wing characteristics: M— 1.7; C/,= 0,2. 
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as a function of lateral to vertical distance ratio for the three values 
of wing dihedral. For this limited amount of data obtained to date 
the overpressures are orderly and only partial agreement with theory 
(ref. 4) has been found. 

The curves on the right of figure 8 are intended to show the atten- 
uation of the centerline maximum overpressures for the three dihedral 
wings out to large distance ratios that are commensurate with flight. 
These results have been derived from the measured centerline sig- 
natures on the assumption that these near-field characteristics would 
remain orderly out to any distance. This assumption needs further 
experimental investigation. Although it might appear from these 
curves that the effect of dihedral virtually disappears at the higher 
distance ratios, a close scrutiny of the results will show that the 
relative effectiveness of dihedral at this distance ratio of 200 is still 
some 60 percent of the relative effectiveness measured in the wind 
tunnel at the near-field distance ratio of 4.5. These results, however, 
are preliminary and a realistic assessment of the significance of wing 
dihedral cannot be made until a more thorough study is completed. 
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The problem of calculating the sonic boom amounts to the problem 
of tracing the discontinuities (shocks) through a piecewise continuous 
function space (the atmosphere as time proceeds). For axisymmetric 
bodies this problem can be formulated mathematically as that of 
finding solutions to the axisymmetric wave equation. The linearized 
form of this equation, modified so that the local characteristics have 
the proper slope and so that local shocks are introduced to satisfy 
uniqueness properties, has been shown (by Whitham) to. give quite 
satisfactory results when compared with experiment. Further in- 
vestigations have shown that airplane configurations of the type now 
being flown at supersonic speeds are sufficiently close to being axi- 
symmetric (or, more precisely, are predominantly described by the 
first terms in an axial Fourier expansion) for the Whitham theory to 
be applicable with good results. A fundamental aspect of such studies 
is the question of whether or not any airplane configurations exist 
which have properties (in particular, produce sonic booms) that can- 
not be disclosed by the predominantly axisymmetric Whitham theory 
as it is now applied. This question has led, in one direction, to the 
study of multipoles, that is, higher order terms in the Fourier expan- 
sion; and in another direction, to the study of second-order effects, 
that is, higher order terms in the governing partial differential 
equations. 

It seems reasonable that the basic premises of the Whitham theory 
are sufficient to serve for a mathematical model that will describe 
with sufficient accuracy the boom signature of any (practical) airplane. 
This means, simply, that the form of the linearized wave equation 
could be used, provided the characteristic surfaces were properly 
traced through space according to the local stream direction and Mach 
number. It is, of course, the coalescence of characteristics of the “ same 
family” that forms the shocks, and the shocks further coalesce to 
form the final N-like wave which is the sonic boom. If a significant 
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unsymmetric twisting of the characteristics can be brought about by 
reasonable airplane shapes, there is a possibility that the effect would 
be significant in studies concerning the minimization of the boom 
pressure impulse. 

Unfortunately, the local adjustment of characteristics in three- 
dimensional flow is not as simple as for the axisymmetric case. The 
characteristics lines, appearing in the Whitham theory, generalize to 
conoids developed about curving streamlines. Further, if the three- 
dimensional properties are significant, the direction cosine between 
the axis alined with the free stream and the local velocity vector will 
depend, to a first order, upon both of the cross-perturbation velocities. 
This effect, together with the locally corrected Mach number, which 
depends principally on the third perturbation velocity, determine the 
regions of coalescence and the shock structure. 

Because of the complications mentioned above, an investigation was 
initiated for developing methods for the direct numerical calculation of 
the flow about airplanes moving at supersonic speeds. Here the word 
“direct” means that the basic partial differential equations are reduced 
at once to difference equations which are programed for a digital 
computer. This concept is not new; nevertheless, at this time the 
numerical methods that are available for such calculations are not in 
a satisfactory state. The principal difficulty lies in the fact that the 
reductions of the partial differential equations to difference schemes 
implies a continuity in the dependent variables that simply does not 
exist for the initial and boundary values that are pertinent. To help 
offset this difficulty a digital computer system was assembled which 
contained a cathode-ray tube connected directly to the high-speed 
core of the computer. Methods are now being constructed which 
will make use of the visual display in an attempt to isolate the shocks 
from the continuous portions of the flow and calculate both by their 
individually appropriate techniques. These developments are in a 
preliminary stage and no results can be presented that apply to actual 
sonic booms. The remainder of the paper is devoted to a brief 
discussion of the numerical analysis that has been developed so far 
under the consideration that immediate graphic representations are 
available. 

One important aspect of adding graphics to the sort of problems 
being discussed is to throw open the possibility that old (precomputer) 
methods may compete favorably in the new environment. In this 
sense, the method of characteristics in three-dimensional supersonic 
flows can be considered to be much less used in present-day computing 
than the fixed-mesh methods made popular by Lax in the midfifties. 
Under the new conditions, it seemed appropriate to reexamine the 
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two possibilities.' This reexamination led to the results discussed below. 
Consider the two-dimensional wave equation 

bu ~ 2 bw 
bx~ L bz 

bw bu 
bx~ bz 

If the method of characteristics is introduced, these reduce to the 
two ordinary differential equations 


d u± <7dw=0 (2) 

which apply along the two characteristics. If the notation 

Ui +1 =u(x+Ax, z) 

and 

Ui+i=u(x, z+Az) 

is adopted, equation (2) can be written 

u” +1 —u” +1 +C(wt +1 —W{+ 1 )=0 
U n i+ l -U n t - x —C{w n i + 1 -w n t . x) =0 

Solve for u” +1 and w" +1 and there results 


«i + 1 — §('!i?+l+«"-i)+^(w?+i— ^ W?-x) 
M)? + 1 =|(w?+i+W?_x)+^,(w?+x— «?_,) 


where the mesh must be chosen such that 


CAx 

A '% 


(5) 


Under the conditions imposed by equation (5), equations (4) are 
identical to those which result from the differencing scheme com- 
monly referred to as the original Lax method. This and similar 
studies showed that when the Lax and Lax-Wendroff differencing 
schemes work best, they do so because they are very close to local 
applications of the method of characteristics. 
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Turning next to a consideration of the stability of a numerical 
method, consider the simplest wave equation 


dtt 

da; - ° ds 


( 6 ) 


One form of numerical approximation is to represent the u dependence 
upon z by evaluating u along a discrete number of lines equispaced 
in z and continuous in x. In such a case the vector w is defined so that 


w T =(ui, u 2 • • ■ u M ) (7a) 

and the partial derivative is replaced by 

^=[B]w+0(Az*) (7b) 

where [B\ is a matrix heavily clustered about the principal diagonal. 
If equations (7) are used in equation (6), there results 





which is a set of M coupled, ordinary, differential equations. The 
solution of such equations depends upon the eigenvalues of the 
matrix —C[B]. If [5] is the matrix tridiag 

r_i±i, ±, i=fi 

L 2Ar Ax 2AxJ 

these eigenvalues are 

cos ■ ■ -M (9) 

The point of this analysis can now be made. The order of accuracy 
in representing the derivative in equation (7b) is no worse than k~ 2 
for any value of /3. However, the stability of equation (8) is pro- 
foundly affected by the choice of B., and one can see from equation (9) 
that equation (8) is stable or unstable depending on whether fi is 
greater than or less than 1. It is customary, in the study of ordi- 
nary differential equations, to regard the eigenvalues of equation (8) 
as fundamental numbers and to associate them with the physical 
interpretation of the phenomena being studied. Here we see that the 
eigenvalues are completely artificial. They may induce damping 
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(“artificial viscosity”), bring about neutral stability (the central 
differencing scheme, for which /5 =0) or they may cause divergence. 
At any rate, they no longer have a fundamental connection with the 
physical problem. 

The arbitrary nature of the eigenvalues of equation (8) depends, 
of course, on the approximation in equation (7b); an approximation 
which, in turn, depends upon the continuity implied in equation 
(7a). In flow fields with shocks this continuity does not exist every- 
where. If it is assumed to exist where it does not, and if the resulting 
errors are forced to decay by choosing (in effect) a value of P which 
induces an artificial viscosity, the results cannot, always, be justified. 
The entire flow field will be forced to decay, and this is not, in general, 
acceptable. 

In conclusion, the ability to compute flow fields about airplanes 
traveling at supersonic speeds with the aid of an immediate visual 
display of the calculations as they proceed opens the possibility of 
devising new, or reviving parts of old, numerical techniques. In 
exploring these possibilities, preliminary studies indicate the following : 

(1) Differencing schemes for hyperbolic partial differential equations 
governing discontinuous functions are good when they approach 
spacings that are equivalent to the method of characteristics, and this 
property is not necessarily connected to their stability. 

(2) Artificially stabilizing methods based on continuous functions 
not only force errors to decay, but also can force the entire solution 
to (artificially) decay. Such methods do not appear to be appro- 
priate for problems involving the sonic boom. 
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Part I — Analysis of Configurations 
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New York University 


INTRODUCTION 

Present configurations selected for supersonic transport produce 
booms having maximum initial overpressure on the order of 2 lb/ft 2 . 
A substantial part of the activities which has been carried out under 
NASA Contract N GR-33-016-1 19 has been directed toward the 
definition of future configurations that can reduce such values of 
maximum overpressure, without compromising too strongly other 
characteristics of the airplane. Several parameters affect the sonic 
boom, such as the weight of the airplane, altitude of flight, dimensions, 
etc. In this analysis the only parameter considered as variable is the 
configuration of the airplane. Therefore, the following assumptions 
have been made in the analysis : 

(1) The length of the airplane is kept constant and equal to 300 
feet. 

(2) The weight, altitude, and Mach number of flight are kept 
constant and equal to 465 000 pounds and 60 000 feet and M— 2.70. 

(3) The variations of configurations investigated must not affect 
the drag too much at the required lift at the flight Mach number. 

The analysis has been limited only to conditions corresponding to 
the maximum Mach number of the airplane. Therefore, the con- 
figurations obtained are possible configurations for cruise conditions 
and usually have acceptable transonic qualities. No attempt has 
been made to analyze subsonic cruise and low speed performance of 
such configurations. 

In order to use realistic parameters in the analysis performed, a 
configuration proposed by industry for the supersonic airplane design 
has been used as a basis for comparison. 

In the analysis, many of the design and aerodynamic data used 
have been selected on the basis of the values of equivalent parameters 
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selected in such configuration. However, the distribution of lift along 
the span has been modified. The analysis presented here shows that 
by conveniently selecting the shape of the equivalent area distribu- 
tion much lower values can be obtained for the sonic boom than 
the lower bound obtained before for similar conditions of flight. Con- 
sider first the sonic boom given only by the lift. The expression 
corresponding to configurations used in reference 1 is given by: 

Ap _0.60K r p l/ * IffW m 

p ■ v 2 {i) 

This value corresponds, for example, to a IS ting surface where the 
lift per unit span is constant and the integrated lift is proportional to 
the length from the leading edge, and there|ore L(r)==i3Lr. A lower 
value of Ap/p given by: . 

A p 0.53 ITVlS 

P ;'(h/Z) m \2 P w 

is possible with a configuration having large lift in the front part of 
the wingspan (ref. 2). Such values correspond to a maximum over- 
pressure of the order of 1.20 to 1.31 lb/ft 2 for the conditions specified 
above, and a value of K r equal to 1.8. The value of Ap/p given by 
equations (1) and (2) . are determined by assuming far-field effects. 
This assumption is correct for the configurations analyzed in reference 
1 and for the configurations selected for the supersonic transport 
designs. However, lower values can be obtained for the conditions of 
flight selected and for the size of the airplane selected by the U.S. 
design. For these conditions, substantial reduction on maximum 
overpressure can be obtained if configurations are used that utilize 
near-field effects. This can be attained by introducing an appro- 
priate distribution of lift that emphasizes such near-field effects. In 
order to prove this possibility, three steps are considered in the 
analysis presented here. First, the present selected configurations 
are reviewed in order to determine the reason for the difference 
between the minimum values given by equation (2) and the actual 
values of Ap/p of practical configurations selected. As a second 
step, the lift distribution and therefore the equivalent cross-sectional 
area distribution for a given lift and length is changed. The third 
step, to reduce the value of max Ap/p , analyses the implications on 
the actual airplane design of the differences between such equivalent 
area distributions and the area distribution used in the actual con- 
figurations. 
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SONIC BOOM ANALYSIS 

In figure 1, the equivalent cross-sectional area of the supersonic 
transport selected for the competition is presented. Curve 1 corre- 
sponds to the complete airplane; curve 2, to the lifting surface alone. 
Curve 3 corresponds to a distribution as required by the values of 
equation (1) for the same lift. 

In figure 2, the sonic boom signatures obtained for the three distri- 



Figure l. — Equivalent cross-sectional area of a present supersonic transport 

configuration. 


lb /ft e 



Figure 2. — Sonic boom signatures for conditions of figure 1. 
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butions are shown. All the values of the sonic boom presented here 
have been obtained with a numerical program developed by H. W. 
Carlson and associates at the Langley Research Center of N ASA and 
made available to New York University. The sonic boom of the 
complete airplane gives a Ap max — 1.92 lb/ft 2 , which agrees with the 
value calculated by others for the same configuration, while the wing 
alone is responsible for Ap of 1.59 lb/ft 2 . The reflection coefficient 
value K ri assumed here and in the following results, corresponds to 1.8. 
These results are far from the optimum considering the lift is con- 
centrated only in a small fraction of the length (about 42 percent), 
while the length of the airplane is 300 feet. To indicate the sensi- 
tivity of sonic boom to the lift distribution and equivalent cross- 
sectional area distribution along the length, selected shapes of several 
equivalent area distribution have been considered as shown in figure 3. 

Curve 2 of figure 3 is close to the distribution required for minimum 
boom. Curve 1 corresponds to a constant distribution of lift distri- 
bution for unit length of the airplane. If C L is the local lift coefficient, 
referred to the span C at the given station considered, then 

U=C L C=2£°AOpiy 

where the quantity ACp is the local pressure rise in the lower surface 
of the wing at a given station y along the span of the position x. 
Such distribution could be obtained by a sweptback wing with sonic 



Figure 3.— Equivalent area distributions corresponding to a given value of 

total lift. 
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edges. Curve 3 corresponds to a wing that is initially triangular and 
has appropriate camber. Curve 4 corresponds to present airplane 
configurations. The sonic booms obtained are shown in figure 4. 
The results are of interest because they indicate that a large departure 
of S(x) distribution from the optimum value affects the maximum 
pressure of the sonic boom only slightly, provided that the lifting 
surface has an actual length equal to 300 feet, and that the lift is 
produced in the front part of the airplane. It must be noted that 
for the lift distributions of curves 1, 2, and 3, the sonic boom signature 
has still near-field effects while the boom signature given by the area 
distribution of curve 4 has no near-field effects in spite of the fact 
that the total lift, the altitude of flight, and length of the airplane 
are the same in all cases. 


1.92 



Figure 4. — Sonic boom signatures for the equivalent area distribution of figure 3. 

The maximum overpressure given by the far-field solution for the 
curves 1, 2, 3, and 4 corresponds to 1.35 lb/ft 2 for curve 1, 1.28 lb/ft 2 
for curve 2, 1.37 lb/ft 2 for curve 3, and 1.92 lb/ft 2 for curve 4. 

These results are important because they indicate that the lower 
values given by equations (1) and (2) for far-field solutions are not 
indicative of the present problem, provided that large amount of 
lift is concentrated in the front part of the airplane. In addition, 
the results show that conventional configurations do not utilize the 
effect of near-field solutions in spite of the fact that the dimensions 
of the airplane are sufficiently large to make near-field effects usable. 
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MULTIPLANE CONFIGURATIONS 

These results indicate that in order to reduce sonic boom, the 
airplane must have additional lifting surface in the front part of the 
airplane and therefore imply the use of highly sweptback wings. 
To have some indication of the minimum amount of lift required 
in the front part of the airplane to obtain near-field effect, a pre- 
liminary investigation has been performed of a tandem lifting system 
having different proportions of lift on the front wing and rear wing. 
The configurations considered are shown in figures 5, 7, and 8. The 
conditions of curve 5 correspond to a lift in the front wing on the 
order of one-sixth of the total. The shape of the equivalent cross- 
sectional area for figure 5 is shown in figure 6. Curve 1 in the figure 
corresponds to the equivalent area of the SST configuration selected 
as a comparison. In figure 7, the lift of the front wing has been 
increased to one-fourth of the total lift and in figure 8 to one-third of 
the total lift. The root cord of the front wing of figure 6 is 45 feet long. 
The wing starts at the nose. In figure 7, the root cord of the front 
wing is 85 feet; in figure 8, the root cord is 95 feet. Figures 5, 7, 
and 8 indicate the sonic boom signature for the configurations shown. 
The results of this analysis give some interesting indications. A 
small amount of lift or equivalent nonsectional area in front is in- 
effective for reducing the sonic boom. In addition, any increase of 
length is ineffective unless the lift is better distributed. The differ- 
ence in signature between two configurations, one 154 feet long and 
the other 300 feet long, is very small as shown in figure 5. However, 
if the lift of the front wing is increased, the value of the maximum 
overpressure decreases as shown in the results of figures 7 and 8. 


(2) (3) (4) 



Figure 5. — Signature of tandem configuration shown. 
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The following interesting conclusions can be reached from the 
results of figures 7 and 8. For lengths of the order of 300 feet, near- 
field signature is obtained if one-fourth to one-third of the lift is used 
in the front part of the airplane. Then, the value of the first shock 
can be kept to values of overpressure of the order of 1 lb/ft 2 . In 
addition, this result indicates a possible direction for producing 
gradual pressure rise without increasing the length of the airplane. 

A comparison of the results in figures 7 and 8 indicate that the A p 
of the first step of the overpressure depends mainly on the lift of the 
front wing when near-field effects are obtained. This value can be 
decreased by increasing slightly the length of the airplanes. From the 
results presented, it can be inferred that a larger airplane of the order 



Figure 6.— Equivalent cross-sectional area of configurations shown in figure 5.. 



Figure 7. — Sonic boom signature of the configuration shown. 
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Figure 8.— Sonie boom signatures of the configurations shown. 

of 350 to 400 feet long can have sonic boom signatures where the 
first step can be below 1 lb/ft 2 . 

PRACTICAL CONFIGURATIONS 

The distribution of lift along the length of the airplane is an effec- 
tive parameter in order to decrease the maximum overpressure 
attributed to the sonic boom. However, the effect of such distribu- 
tion on the drag attributed to the lift must be investigated. Usually 
in order to maximize the L/D of the airplane, two steps are considered. 
First, the drag for zero lift is minimized, then the drag attributed to 
lift is minimized. In order to reduce the sonic boom, interference 
effects must be utilized. The introduction of lift in the front part of 
the airplane makes the equivalent area distribution similar to the 
cross-sectional area distribution of a blunt body. This effect could 
be erroneously interpreted as equivalent to a corresponding increase 
of the drag attributed to lift. That this is not so can be easily seen 
by using the theorem of equality of drag attributed to lift in direct 
and reversed motion. The theorem states that within the approxi- 
mation of the linear theory, the drag of a given distribution of lift is 
unchanged by a reversal of the direction of motion. The reversal of 
the direction of motion for a given lift distribution produces large 
variation of the equivalent area distribution in spite of the fact that 
it does not change the drag. This effect is illustrated in figure 9. A 
triangular wing is considered in the upper left side of figure 9. When 
the direction of flight without changing the distribution of lift is 
reversed, a wing having the same spanwise and cordwise distribution 
of lift is obtained, figure 9 left. Then the drag attributed to lift does 
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Figube 9. — Equivalent area distribution of a given lift distribution in direct 
and reverse motion. The drag attributed to a given distribution of lift is 
unchanged by a reversal of the direction of motion. 

not change by reversing the flow, and the value of 0 L C of the two 
wings remains the same in the two cases because the lift distribution 
must remain the same. The equivalent cross-sectional area for the 
two cases is quite different in spite of the fact that the lift and drag 
are the same because the direction of integration changes sign. This 
consideration indicates that the drag attributed to lift does not 
introduce any limitation to the requirement of large lift in the front 
part of the airplane. 

The other consideration that can be made is related to the effect 
of volume. Again, the contribution of volume can be incorporated 
in the effect of lift. The fuselage, if considered as an isolated body, 
can be represented by a line source sink distribution. If the fuselage 
is placed underneath a flat plate, with sonic edges, then the source 
distribution induces lift on the plate proportional to the source’s 
strength, while the back part of the fuselage is placed above a plate, 
then the sink distribution induces lift proportional to the strength 
of the sink. This can be seen immediately by using Evvard’s theory. 
The sources and sinks required for the lift can be used to re- 
distribute the equivalent area of the airplane. In figure 10 the re- 
quired source sink distribution for the transport is indicated. If the 
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fuselage is placed under a wing and the leading edges of the wing are 
sonic, then the inviscid LjD of the system for a given C L is not sub- 
stantially changed. This is shown in figure 11. In one case a 5° 
cone is placed symmetrically across the wing or in front of a wing. 
In the second case, twice the same cross-sectional area is produced 
on the same length, by increasing the cone angle, and only half of 
the cone is used and placed underneath the wing. The pressure 
field of the cone extends under the wing and produces lift. The 
curves show that L/D in the absence of viscous drag is about the 
same in the two cones, in spite of the fact that the drag for zero lift 
is higher for the second case. The presence of the fuselage under the 
wing increases the lift of the front part of the lifting surfaces. The 
fuselage is closed in the rear on the top of the wing as indicated in 
figure 12. Assume that the wing area considered for the configuration 
is the same as for the configuration selected for the supersonic airplane. 
In the SST configuration used, the wing area is of the order of 8000 
ft 2 , while the fuselage wetted area is of the order of 8300 ft 2 . In the 
configuration of figure 12, the wetted area of the fuselage is reduced 
substantially, because the fuselage is incorporated with the wing and 
the additional wetted area attributed to the fuselage is reduced 




Figure 10.— Schematic utilization of volume to produce lift. 


Schematic of U tilization area, ft* 

Lift to be produced .... 780 

Sources , „ _ _ . - ...... 780 

Sinks. ........ . _ _ ___ . . _ . ... _ - _ - - , 780 

Sources available attributed to thickness. . 220 

Sinks available attributed to thickness. 220 

Total attributed to lift if interference is utilized: 

Sources ... 560 

Sinks..., _______ 560 
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Figure 11. — LfD attributed interference. 
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COMBINATION DOES NOT CHANGE DRAG 
LIFT HIGHER IN FRONT 

WING 



TOTAL DRAG 

DRAG, OF WING C+ FUSELAGE 


Figure 12. — Schematic design of fuselage integrated with lifting surfaces. 

substantially. Therefore, if the wing area is kept the same, a reduc- 
tion in drag can be obtained. An analysis indicates the wing area 
could be increased by 3000 ft 2 to a value of 11000 ft 2 without increas- 
ing the drag at a given lift. Figure 13 gives the polar diagram for 
the SST configuration. The wave drag of the airplane is a small 
part of the drag and is of the same order of the skin friction drag 
of the fuselage; therefore, the reduction of skin friction drag of the 
fuselage can balance some increase in wave drag. In addition, any 
decrease of wetted area corresponds to a decrease in weight of the 
airplane. If the surface area is increased and the total wetted area 
of the airplane is kept constant, then the shock drag of the wing is 
decreased because the percent thickness of the wing decreases, and 
the drag attributed to lift is also decreased. For a given weight of 
the airplane, C L S is constant (S is the area of the wing). The total 
drag of the wing is given by three components. 
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SCz> Wing — $Ci>wave~f“$ C& friction ~{~S C& lift 
and for constant lift and span, it is 

S Ci> wing^-K/ 8 -f* Cx> friction. S 

For actual conditions, an increase in S decreases the total value of 
$C D W i ng 



POSSIBLE CONFIGURATIONS 

In order to indicate possible configurations, the analysis of three 
airplane designs is presented here. These configurations have been 
selected on the basis of being not too different from present con- 
figurations. Many other configurations having more radical changes 
could be considered. The configurations have approximately the 
same L/D at M= 3. The configurations are shown in figures 14, 
16, and 17. 

Figure 14 is a triangular wing airplane with higher sweptback 
wing than present configurations. The fuselage and wing thickness 
are used to increase the lift in the front part of the airplane. The 
sonic boom of such an airplane is shown in figure 14. The maximum 
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Figure 14. — Possible configuration for supersonic airplane. 


overpressure is of the order of 1 lb/ft 2 . The lift distribution of the 
airplane is indicated in figure 15. The quantity C L is an average 
value integrated along the span of the wing (not along the cord). 
Figure 16 indicates the same configurations having larger wing area 
and the same sonic boom. 



Figure 15. — Distribution of lift along the length. 
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Figure 16. — Schematic design of airplane having large wing area. 


Figure 17 indicates another possible configuration, and the corre- 
sponding sonic boom. The same type of C L C distribution is used 
in this configuration. 



ib /ft z 
1.03 



Figure 17. — Possible configuration for low boom supersonic airplane. 


/ x SUMMARY 

Configurations have been investigated for supersonic transport 
which meet present SST design and requirement. It is shown that 
optimization on the waves of far-field analysis is not satisfactory. 
For a 300-foot airplane length, near-field effects can be obtained 
when sufficient lift is generated near the nose of the airplane. Because 
of near-field effects, sonic boom maximum overpressure on the order 
of 1 lb/ft 2 can be obtained with configurations that do not penalize 
the L/D of the airplane at cruise. Large improvements can be expected 
by limited increase of airplane dimensions. 


315-766 0—68 7 
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CONCLUDING REMARKS 

From these preliminary considerations, it appears that airplanes 
having lifting surfaces extended to the front of fuselage, and utilizing 
interference effects can be effective in reducing maximum sonic 
boom overpressure, because it utilizes near-field effects. | 
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INTRODUCTION 

When an N-wave hits a structure, two aspects of the reaction are 
of interest. One is the intensification factors and the locations where 
the transient pressure. variation becomes a maximum. The other is 
the assessment of the. transient load or pressure distribution on the 
structure. The .transient maximum pressure rises usually occur at a 
concave corner facing the incident N-wave. The intensification factor 
at the corner in certain instances can be determined directly by 
making use of the reflection principle, the conical solution (ref. 1), or 
the averaging principle (ref. 2) without the knowledge of the solution 
to the entire flow field even its neighborhood. However, for the 
determination of the transient load on the structure, it is necessary 
to obtain the history of the pressure distribution on the structure. 

A procedure for the determination of the transient pressure distri- 
bution on a two-dimensional structure in the shape of a rectangular 
block by an incident N-wave is developed in the next section by the 
superposition of the solution of a unit plane pulse incident on the 
structure. When the plane pulse hits the first convex corner, the 
pressure distribution is given by a conical solution (ref. 1) . When the 
conical flow field propagates to the next convex corner, the diffracted 
flow field, which is no longer conical, is obtained by the method 
developed for the diffraction of a general disturbance by a right 
convex corner (ref. 3). This method has to be employed whenever a 
disturbance attributed to diffraction at a convex corner propagates to 
its adjacent convex corner. A numerical program is developed to 
compute the pressure distribution on the two-dimensional structure 
incident by a plane pulse or by an N-wave. 

The problem of an N-wave incident on a three-dimensional structure 
is very complicated. Even the basic conical solution of a plane pulse 
incident on the corner of a cube is not yet available. In the last 
section, the transient intensification factors at several corners of a 
three-dimensional structure are computed when the averaging prin- 
ciple (ref. 2) is applicable. The results are used to demonstrate the 
dependence of the intensification factors on the direction of the 
incidence. 


89 



90 


SONIC BOOM RESEARCH 


TWO-DIMENSIONAL PROBLEM 

When the structure is represented by a rectangular block in the x-y 
plane and the incident waves can be represented by superposition of 
parallel plane pulses inclined at an angle a to the horizon and with its 
normal in the x-y plane, the pressure disturbance and the velocity 
potential are governed by the linearized theory, the simple two- 
dimensional wave equation 

~Q 2 <Plt <pxx <Pvy = 0 

where C is the speed of sound at the ambient condition. The shape 
of the N-wave, i.e., the function G(Ct)j is assumed to be given from 
the sonic boom calculations where G is the overpressure nondimension- 
alized with respect to the strength of the front shock and Ct is the 
distance behind the front shock. If represent the overpressure 
on the surface of the structure at an instant f>0 attributed to a unit 
plane pulse incident at the front top corner o at the instant t — 0 
(fig. 1), the overpressure on the structure attributed to the N-wave 
is given by the integral 

am 

=1bm-T t )[G{CT t )} 

t— 1 

n Pt—Ti+O 

+oe G'ict-cnj&ndt' (i) 

t — l •/ O 

where s is the arc length along the front, top, or the rear sides; n is 



Figure 1. — Pattern of diffracted waves. 
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number of shock waves in the N-wave; and [ G(Ct *)] is the jump of 
G across the ith shock at r*. According to the definitions, it should 
be noted that [G{Ct1)} — 1 and f(s,f) ~ 0 for << 0. The problem now 
reduces to the determination of the function the overpressure 

attributed to the incident of a unit plane pulse. 

Immediately after the plane pulse incident at the corner o (fig. 1), 
i.e., £>0, the diffracted wave is given by a conical solution, pco(x/t,ylt) 
(ref. 1), originated at 0. The plane pulse and its reflected wave travel 
across the top surface and arrive at the rear corner A at the instant 
ti—L sin a/O. For t^>t h the diffracted wave at the corner -A is given by 
another conical solution, p c i [(x—L)/(t—ti), y/(t— *i)], originated at 
point A. 

At the instant t 2 —L/C, the conical solution p co arrives at the 
comer A . For C>t 2 , the diffracted wave at corner A is no longer a 
conical solution. The solution for diffraction of a general incident 
wave by a convex right corner was presented in ref. 3. To facilitate 
the numerical evaluation of the solution, the incident conical solution 
Pco along the extended top surface can be approximated as 

Pco(%/t, 0)—a 0 x/(tc)+bo[l—x/(tc)f 12 

The conical solution and the coefficients a Q , b 0 depend on the incident 
angle a. The diffracted waves can be written as a 0 Si ( x , y, t) +b 0 S u 
( x , y, t)j where Si and S n are the parts attributed to ^l—x/(icj and 
[1 — x/(tc)f /2 , respectively, and are independent of the incident angle 
a. S x and S u can therefore be computed by subroutine. At the 
instant U—L (1+sin a)/C , the conical solution p C \ originated as 
corner A travels back across the top surface and reaches corner 0 . 
For the conical solution along the extended top surface can be 
approximated by 

Pci[ (x-L)/(t-td, 0]^T^x'l{a') + b 1 [l-x'/(Ct')r> 

where x' — (L—x) and t'—t—U and and b x are constants depending 
on a. The diffraction at corner 0 for t^>t 3 is then given by biS u 

with appropriate shift of time scale, and shift and rotation of the 
x and y variables. 

When t A — ( L sin ot+ITj/C, the conical solution p ci travels across the 
rear surface and reaches the ground. For 0>(L sin a+2H)/C , the 
reflected wave travels upward and passes the corner A . The addi- 
tional diffracted waves can again be represented by the linear combina- 
tion of the fundamental solutions S t and S n with appropriate changes 
in the variables x , y, and t. 

It should be pointed out that for 2>0, the plane pulse and its 
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reflection will travel down the front surface, be reflected by the ground 
when t^>H cos a/O, and be diffracted by the corner 0 when £>2 H 
cos a/G with the addition of a conical solution. The conical solution 
Pco will travel down the front surface, be reflected by the ground for 
0>H/G, and be diffracted by the corner 0 again when G>2H/C. It 
is clear that the intensification factor at the lower front corner B is 4 
for the duration H cos a/C<t<H/C. 

With the conical solutions and the fundamental nonconical solutions 
computed by subroutines, the additional solution attributed to 
diffractions and reflections at corners 0, A , B , and D can be computed 
and added at the appropriate instant, and the numerical program for 
overpressure f(s,t) on the surface of the structure is thereby con- 
structed. 

Figure 2 shows the pressure variation j(s,t) on the surface of a 
square block (L—H) attributed to a unit plane pulse incident at 45° 
(a— 45°). From the discontinuities in the curves, the propagation of 
the waves and the additions of reflected or diffracted waves can be 
observed. As t increases from zero to ZL/C, the pressure variations 
on the front, top, and rear surfaces increase from zero and finally 
approach the uniform value of 2. 

With the aid of the solution the pressure variation on the 

same block attributed to an N-wave incident at 45° to the same 
structure can be obtained by equation (1). Figure 3 shows the 
numerical results for an N-wave with a front and rear shock of the 
same strength; the thickness of the N-wave is L , i.e., n= 2 , [G(0)]= 
[(7(1) j—1, and G'——2 for (?<CV<Z. The pressure variation 
gradually returns to the ambient value for 0>ZL/C. 

With the same procedure, the pressure variation on a rectangular 
block (. L—2H) attributed to the incident of the same N-wave is 
computed and shown in figure 4. 

Figure 5 replots the pressure variation at the front bottom corner B 
from the preceding two figures to demonstrate the influence of the 
geometry of the structures. The front shock arrives at t=H cos a/O 
and is followed by a linear expansion wave. Further expansion occurs 
after t=H/G , when the diffracted waves from the top corner 0 arrive. 
The rear shock arrives at t= (H cos a+L)/G because Cr 2 —L . After- 
ward, the pressure returns to the ambient value rapidly. The basic 
differences in these two curves are attributed to the shift of time 
scale for the different ratio of HjL. Additional small differences are 
attributed to the different time of arrival of the diffracted waves 
from the top corner 0. 












Figure 5. — Pressure variation at the front bottom corner. 
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THREE-DIMENSIONAL PROBLEM 

Because the near-field solution for the diffraction of a plane wave 
by a three-dimensional body is not yet available, only some partial 
information at corners will be obtained by the averaging principle 
(ref. 27). The principle states that the intensification factor at the 
vertex of a cone with solid angle fi is equal to the ratio of the solid 
angle of the flow field at the vertex without and with the cone, i.e., 
471-/(47!-— S2). In the proof, the conical properties of the boundary 
are employed, but the solution of the flow field is not restricted to 
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Figure 6— Three-dimensional diffraction of incident wave. 


conical flow. Near the vertex or corner of a body it is locally a conical 
surface, and the averaging principle can be applied to compute the 
overpressure until the arrival of the disturbances from the nonconical 
portion of the body. This principle is now applied to a V-shaped 
building with a concave angle X (see fig. 6). For a given airplane 
flying at given Mach number and altitude, the N-wave signature and 
the incident angle to the ground are the same. The diffracted waves 
and the pressure distribution on the building depend also on the direc- 
tion of the flight path. Three distinct types of flight path are con- 
sidered, and the intensification factor at the corners and vertices are 
computed by means of the averaging principle. 

Flight path I approaches the concave side of the V-shaped building 
and makes a finite angle with the two legs of the V. Flight path II 
is parallel to one leg of the V- Flight path III is opposite to that of 
I and is approaching the convex side of the V. The ratio of overpres- 
sure to the value of incident N-wave for various corners and vertices 
is tabulated also in figure 6. 

For paths I and III the reflections of waves from the ground and 
the two vertical surfaces are assumed to be regular reflections (ref. 4) . 
For path II, there is no reflection from one of the vertical surfaces, 
i.e., a degenerated case of Mach reflection. The intensification factor 
at the lower concave corner of the V is 4?r/X for flight path I and is 
2 w/X for flight path II; they are 8 and 4, respectively, when X=7r/2. 
In between flight path of types I and II it will be necessary to admit 
Mach reflections from one of the vertical walls. On the other hand, the 
intensification factor at the convex lower corner for flight path III 
is 47r/(27r — X), i.e., % for X=7r/2. The lower concave corner is much 
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lower than 4ir/\ times the value of incident wave. These examples 
clearly illustrate the importance of flight path direction. It should be 
pointed out again that the intensification factors, i.e., the ratio of 
overpressure to that of incident wave, given in the table are valid 
prior to the arrival of diffracted or reflected waves from the adjacent 
edges and corners. For a complete pressure distribution on the 
surfaces of the buildings it is necessary to develop the near-field solu- 
tions. The first step will be the construction of a conical solution 
for the diffraction of a plane pulse by a convex corner of a cube. 


r 


CONCLUSIONS 

A numerical program for the pressure distribution of an N-wave 
incident on a two-dimensional structure in the shape of a rectangular 
block is completed. For a three-dimensional structure, the averaging 
principle is employed to compute the intensification factors for 
corners and vertices for certain durations and to demonstrate the 
importance of the flight path direction. j 
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Reduction of Sonic Boom Attributed to Lift 

E. L. Resler, Jr. 

Cornell University 


INTRODUCTION 

The expected sonic boom associated with presently planned super- 
sonic aircraft, particularly the commercial supersonic transports, is 
not likely to be received graciously by the general public. To alleviate 
this problem, ways to modify the aircraft's configuration to reduce 
or eliminate the sonic boom would be highly desirable. Physical 
arrangements to eliminate completely the boom attributed to lift 
are discussed in this paper. The particular configurations discussed 
derive their lift by imparting downward momentum to the fluid in 
the vicinity of the aircraft and canceling the far-field effects by other 
means, namely, a reduced area stream tube that could be accom- 
plished by special engine designs. These principles are most easily 
demonstrated in the two-dimensional case and then extended to 
three-dimensional cases. 

TWO-DIMENSIONAL CONFIGURATIONS 

The earliest and most famous example of an aerodynamic shape 
with thickness but no wave drag, and thus no boom associated with 
it, is the Busemann biplane and its counterpart with cylindrical 
symmetry. These configurations exhibit no boom because all the 
waves are confined to the internal flow and canceled by surface 
contours. The linearized version of a Busemann biplane is shown in 
figure 1. To continue the discussion, it is convenient to make use 
of the theory of H. Lomax (ref. 1). This reference relates equivalent 
cylinders in two dimensions or bodies of revolution in three dimensions 
to the cross-sectional area and normal forces on any aerodynamic 
configuration. In the far field the pressure disturbance attributed 
to the configuration is the same as would be caused by the equivalent 
bodies. The equivalent bodies may be different at different azimuthal 
angles depending on the configuration. Two-dimensional flows are 
particularly simple in that there are only two directions or angles of 
consequence, below (— 7r/2) and above (+tt/2). 
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Using the theory of Lomax to analyze the two-dimensional Buse- 
mann biplane makes clear the reason for no wave drag in terms of 
equivalent cylinders. The forces acting on the biplane give an equiv- 
alent cylinder that identically cancels the equivalent cylinder because 
of thickness. There is another way to visualize the equivalent cylinders 
for the Busemann biplane which will be used later to discuss lifting 
bodies. Consider the stream tube that passes through the biplane 
and consider it a solid cylinder that extends from — <» to + 03 • 
This cylinder plus the solid parts of the Busemann biplane when 
superimposed (following the Lomax scheme) result in an equivalent 
cylinder of constant area that extends from - » to + 00 and again 
such a cylinder obviously does not disturb the flow at all and has 
neither wave drag nor boom (fig. 2). 

Of course, there is another class of bodies that exhibit no boom but 
do have wave drag. Such bodies can be constructed from the equiva- 
lent bodies of the Busemann biplane. Such a body, which from the 
ground still looks like an infinite cylinder, is shown in figure 3. Any 




Figure 2.— Equivalent cylinders according to area rule for Busemann biplane. 
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Figure 3. — Body with thickness, no sonic boom, but wave drag and equivalent 
cylinders from ground (— ti-/2). 

two-dimensional shape with a streamline that is straight and parallel 
to the ground on the under side will give no two-dimensional boom 
but will have a wave drag. It is apparent that there are also three- 
dimensional counterparts of this class of body, namely, bodies with 
thickness, that are not hollow and would cause no pressure disturb- 
ance, and thus no boom directly beneath them, but would have a 
wave drag. These configurations would exhibit a moment but no 
net lift. Because skin friction is very important for supersonic 
flight conditions, the wave drag of these bodies may be a reasonable 
price to pay for boom reduction without having to cope with the large 
friction drag of a Busemann-type configuration which accomplishes 
the same boom reduction. 

Suppose we now extend these considerations to construct a lifting 
body that would exhibit no boom. Again the two-dimensional case 
is most easily discussed. Begin with a flat plate at angle of attack 
which ordinarily would have waves extending above and below it. 
It is the downward directed waves that result in the two-dimensional 
boom. Thus the front compression wave must be intercepted with 
another plate as shown in figure 4. In this case, to maintain a 
horizontal streamline, the air in the stream tube between the plates 
must be processed to prevent “phiming” of the stream tube as it 
leaves the duct, which would cause a boom. The processing of the 
air in the cross-hatched region of figure 4 is such that the stream tube 
leaves the configuration at the correct ambient pressure (therefore 
preventing pluming) but with a reduced stream-tube area. The 
equivalent cylinders for this two-dimensional lifting configuration are 
shown in figure 5. These are found by replacing the stream tube 
where the air is processed by a solid cyclinder before applying the 





Figure 4. — Two-dimensional lifting configuration with no sonic boom. 

area rules, because these rules were derived for flows where the 
Bernoulli constant is the same on all streamlines and the processing 
will change the Bernoulli constant in that particular stream tube. 
The equivalent cylinder because of lift just fills in the reduction in 
stream-tube area because of the processing and from below (■— w/2) 
the equivalent cylinder is, as in the case for the biplane, one that 
extends from — to + 00 and has no boom. Above the configura- 
tion, however (+ir/2), the cylinder is different from the biplane and 
is shown in the figure. Thus in two dimensions the configuration 
discussed has lift but there is no boom resulting from it. 

THREE-DIMENSIONAL CONFIGURATIONS 

The three-dimensional bodies meeting the same conditions are 
now easily visualized in the following manner. Assume a body of 
revolution extending from — <» to + <» of constant area. This body 




Figure 5, — Equivalent cylinders according to area rule for configuration in 

figure 4. 
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of revolution obviously has no boom associated with it. From t his 
body of revolution remove the stream tube that contains the fluid 
that will be processed so that it will leave the configuration with a 
reduced stream-tube area. The difference in volume between the 
initial assumed body of revolution and the stream tube of reduced 
area can be used for a “lift body” so that from the ground at least 
there is no disturbance. This lift can be achieved with an ordinary 
wing located, of course, properly in accord with the area rules relative 
to the affected stream tube. This whole process is indicated sche- 
matically in figure 6. Thus boom attributed to lift can definitely be 
reduced by programing engine parameters appropriately if the engine 
is used to accomplish the required stream-tube area reduction. 

In principle there is a class of configurations which have lift and 
no boom. Because this is accomplished by reducing a stream-tube 
area, in general there is a high velocity jet leaving the configuration. 
It should be realized that eventually dissipative processes will convert 
the kinetic energy of the jet into heat, and, as a result, the jet will 
eventually plume out to a larger area than at its origin at — °° . How- 
ever, because it is the dissipative processes that control this rate of 
energy conversion, it is expected that this pluming will take a very 
long distance in terms of aircraft lengths. Therefore while the pluming 
resulting from an expanding stream tube will cause a pressure dis- 
turbance on the ground beneath the aircraft, the rate of expansion 
is so slow that the equivalent body length will be very long and the 
ground disturbance will hardly be perceptible. 

AERODYNAMIC PROPERTIES OF BOOMLESS CONFIGURATIONS 

It is interesting to note that the configuration discussed above 
does not give the characteristic N-wave in the far field at any azi- 



Figuhe 6 . — Construction of three-dimensional configuration with lift and no 
sonic boom from equivalent body of revolution according to area rxile. 

315—706 0—68 8 
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muthal angle. In the far field at angles other than — -ir/2 (ground), 
the reduced stream tube more than compensates for the “lift volume.” 
Thus the far-field pressure signature is an expansion followed by a 
shock and then the tail wave. For the case where the lift boom is 
completely canceled on the ground ( — tt/2) , the shock wave is strongest 
above the configuration (+*Y2). 

It is also of interest to note that there is a drag penalty that is 
paid for the boom cancellation on the ground. Referring again to 
the two-dimensional case (fig. 4), it is obvious that the drag of this 
configuration is twice that of an ordinary flat plate at angle of attack 
giving the same lift. (Or alternatively the drag is the same as for 
the flat plate without the channel but one-half of the lift has been 
canceled by the bottom plate pressure rise, or the downward moving 
wave and its associated downward momentum has been intercepted.) 

In the three-dimensional case, as suggested by Hayes at this con- 
ference, the factor is 3. That is, the drag for a given lift using this 
scheme is three times what it would be for an ordinary wing with 
the same lift. This can be computed as follows. Using Lighthill’s 
result (ref. 2) that the wave drag of a body of revolution is just equal 
to an integral of the square of the F-i unction and Lomax's relation 
(ref. 1) between the lift and the equivalent body of revolution, namely 


S' e '(y,e)=-£l'(y,e) 

where 6 is the azimuthal angle, j8=- v / M 2 — 1, g is the dynamic pres- 
sure, and V{y,0) is the streamwise (^-direction) gradient of the lift 
force in the proper obliquely cut sections and S" the second deriva- 
tive of the area of the obliquely cut sections. The F-function F{y) 
for no boom on the ground is then: 

(Lighthill’s notation) 


ei/\ _! p / _£ m sin 0 dt _ £ l(t) d t \ 
{ - y, ~2irJ 0 \ 2 q (y-t) 1 ' 2 2q(y-ty' 2 f 

Thus for 6==—r/2, F{y)— 0. Note also that 

F(y) oc(sin 0+1) 
so 

F\y) oc(sin 0+1) 2 


( 1 ) 

( 2 ) 

(3) 


and integrating from 0=0 to d=2r gives a drag proportional to 3x. 
For the lifting element only which means keeping only the first term 

fir 

l sin 2 6 dd=7r. Thus the ratio of the drag of the boomless lifting 
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configuration to that of the wing above which gives boom on the ground 
is in the ratio of Stt/tt—S regardless of the lift distribution with y. 
The same result is evident from George's work (ref. 3) using his 
formula (16) and in his notation noting that for no boom B t =Ao 
(eq. 13) so the wave drag is just 

D W =2D 0 ^+D^ and 2? 0 <A) =A (B) (4) 

In this formulation A (B) is just the drag attributed to a pressure 
dipole or the wing alone. Again the ratio is just 3, because this is 
the same calculation. 

Fortunately this penalty is paid in wave drag alone. The drag of 
supersonic aircraft is chiefly viscous drag and the penalty paid in 
wave drag alone may not be large compared with the total drag. 

ENGINE DESIGN 

In another related paper (4) the possibility of using an engine to 
accomplish the stream-tube reduction is discussed. In particular, a 
ducted fan was used to effect the area reduction. By setting the 
thrust of this unit equal to the drag of the aircraft, it is found that the 
area reduction A A possible by this means is just 


A r lAAn/Ai 

In equation (5), A f is the inlet area of the unit and 


Ar ~ (L/D)yp { M 2 

where W= weight, £/D=lift-drag ratio, 7 = ratio of specific heats, 
Pi= ambient pressure at cruise altitude, and M = flight Mach number. 
Note that even if A t -^ » a A is at most A R , so a limited amount 
of reduction is possible with this scheme. It is also shown in this 
reference that in order to produce the power aboard the aircraft 
with a turbine, the turbine exhaust stream- tube area necessarily is 
larger than the inlet stream-tube area. This effect reduces the net 
possible area reduction. If LA' is the area increase attributed to 
producing the fan power, one finds: 


.LA 1 

A r 2 1 A A! T 1 

7^1 v W~Al~0 


( 6 ) 


In this equation n is the turbine cycle efficiency, T t is the turbine 
blade temperature and 0 is the ambient air temperature at cruise 
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altitude. To be able to alleviate the boom it is necessary that 
AA^>AA' which implies that the cycle efficiency tj be 

,>( 7-1 for AA>AA' (7) 

The efficiency of the cycle used here is related .to the temperature 
rise through the compressor ahead of the turbine, namely 7 ?= 1 -d/T c , 
where T c is the temperature of the air out of the compressor. Thus 
equation (7) dictates a certain compression ratio for the compressor. 
This analysis also identifies the important parameters involved in 
attempting to reduce sonic boom with an engine. Of course there 
are other schemes, some of which are better than the one described 
here. 

SUMMARY 

It has been demonstrated here that there are a class of supersonic 
configurations that have lift but no sonic boom. Using these prin- 
ciples, the sonic boom attributed to lift can be alleviated by inte- 
grating the engine design with the design of the rest of the aircraft. 
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The Approach to Far-Field Sonic Boom 

F. K. Moore and L. F. Henderson 
Cornell University 


INTRODUCTION 

In this study we have considered the approach to far field as 
described by the Whitham theory (ref. 1). Assuming that the 
observed signature would be nearly an N-wave, we have asked how 
the Whitham theory describes the first departure from that result as 
one moves toward the body. The hope is that, if attention is restricted 
to the approach to far field, signature improvements could be related 
in a simpler way to the aircraft configuration than would the descrip- 
tion of changes throughout the entire flow field. If some such 
simplicity of influence could be found, perhaps effects of high Mach 
number, density stratification, and the like could also be simply 
expressed. It should be emphasized that in this problem the Whitham 
theory has been used without modification, The question is one of 
efficient presentation of results. 

SYMBOLS 

p static pressure 

h miss distance 

l length of body 

y ratio of specific heat 

M Mach number 

P t/AP—1 

F pressure on equivalent body 
I SIF(t) dr 

r dimensionless distance along body axis 
k function of M (ref. 2) 

e small parameter, equation (3) 

e transformed version of e 

E small parameter, equation (15) 

x streamwise coordinate 
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Subscript 

o evaluation where F=0 (smooth body) or just ahead of disconti- 
nuity of F passing through zero. x 

EXPANSION FOR APPROACH TO FAR FIELD 

For purposes of discussion we may limit attention to the strength 
of the shock compression as given by the Whitham theory: 


y 


A p /§\ 8/4 

p \lj 



1/2 


0 1 ' 4 


bow 

shock 




( 1 ) 



bow 

shock 


2 1(f) 

kF\r) 


where F{r) is the customary integral of the curvature of an equivalent 
area function. In equations (1), r can be regarded as the parameter 
enabling connection to be made between pressure jump and “miss 
distance.” The ultimate N -wave is achieved when r approaches 
r 0 , the axial location where the ^-function passes through zero and 
the characteristic staying parallel to the free-stream characteristic 
originates. 

To seek the description of the approach to the far field it would 
seem reasonable to expand the foregoing solution about r 0 , but 
eliminating explicit dependence on r in the result. The following 
expansion results from this procedure. 


V7T + 2/. « 

K_Fo 

P F'i ZF' 0 ZFC 
F 0 2/ 0 + 4 7? 

1 FI 

Iq Iq 

K_1K 

2 II 


_ 

s F 0 2 Iq j 



* 2 +... ( 2 ) 


where 



Subscripts 0 in the foregoing equations refer to quantities evaluated 
at to. 
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Clearly, e vanishes as r— »r 0 , in view of equations (1), and I — >/„ 
as it ought. Thus, e is the small quantity governing approach to 
the far held, whether the equivalent body is smooth near r 0 (so that 
F o=0) or has a shoulder (so that F 0 has a positive value just ahead 
of the shoulder, as in fig. 1). 



If one has a smooth body for which F o —0, the asymptotic series 
just described contains no linear term in e. For example, the first 
two terms of equation (2) would be 


VT" 1+ »« 


( 4 ) 



no 
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This result is consistent with the equation cited by Lighthill (ref. 2) 
for pressure gradient in the far field, when the spreading of the N- 
wave is taken into account. If the body is not smooth, the linear 
term is present, and the N-wave is achieved at the finite distance 
defined by «— >0. 

REDUCTION OF SHOCK STRENGTH— SMOOTH BODY 

For the moment, concentrating on the smooth body, we may ask 
if it may fairly be concluded that the strength of the shock given by 
equation (4) in the “almost-far” field may be reduced by making 
F 0 ' small (it is negative anyway). That is, we consider giving the 
F-funetion a pronounced S-shape near r 0 . A little study shows 
that this approach is not fruitful. The reason is that, apart from 
convergence problems, equation (4) is physically misleading in the 
following sense: If F 0 ' were zero then F would tend to remain near 
zero in the vicinity of t, and the inclination of nearby characteristics 
would tend to remain the same as that emanating from r 0 itself. This 
group of characteristics would therefore not overtake the shock, 
because they would not be spreading. Thus, although the signature 
would tend to have an S-shape, the inflection region would be confined 
to a band of characteristics which, going out from the body, would 
be a progressively smaller proportion of the signature width. 

REDUCTION OF SHOCK STRENGTH— DISCONTINUOUS BODY 


When the equivalent body is not smooth (F 0 5 ^ 0), then the series 
of equation (2) would appear to have a more useful application. Of 
course, when F passes discontinuously through zero, one expects an 
N-wave at a specific distance from the body, namely, when the 
numerator of € goes to zero. Figure 1 is intended to illustrate this 
process for a cone cylinder. As explained in LighthilPs article, one 
may locate the beginning of the N-wave at the confluence of charac- 
teristics emanating from the endpoints of a line which cuts equal 
areas ahead of and within the F-curve, as shown shaded on the sketch. 
The slope of the line which cuts equal areas is y 2 Fo 2 /I 0} by simple 
geometry. Next, we ask if boom reduction can be expected when the 
denominator e is small and negative: 


FJ 


1FI 

2 Iq 


In effect, one considers tilting the F-curve to reach r 0 at nearly the 
same slope as the line cutting equal areas. Figure 2 illustrates this 
sort of change, in which, by the way, the area / 0 is kept the same. 
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Although again, we ean hardly expect convergence when the 
denominator of « goes to zero, we note that because there is a gra- 
dient of F, and F determines the spreading of the Mach waves, they 
can in fact undergo mutual interference. Figure 2 sketches this 
effect. The expansion emanating from r 0 has a higher F, and therefore 
will spread more rapidly than will the compression emanating from 
just ahead at a lower value of F. Thus the compression will prevent 
the expansion from the shoulder from reaching the shock as soon as 
it otherwise would have done, and the N-wave formation is delayed. 



Figure 2. — F-curve and wave system for perturbed cone cylinder. 
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In consequence, if 7 0 is conserved, the strength of the bow shock at a 
given location will be diminished. The overall energy in the N-wave 
would not be changed, though at a given location the shape would 
presumably be more S-shaped. It seems, then, that placing a com- 
pression ahead of the expansion corner may diminish the bow-shock 
strength for this shape. 

IMPROVEMENT OF CONVERGENCE 


Clearly, equation (2) does not converge fast enough to be useful 
in general circumstances. By considering certain specific cases, we 
can find an expansion parameter which gives better convergence 
than €. We must distinguish, of course, between smooth and discon- 
tinuous F-functions. 

(1) Three-quarter-power body and cylinder . — -For this configuration, 
the F-function is a square wave, F~F 0 , between r— 0 and r 0 . The 
exact formula for yj in this case is 


virv/oVsW 


(5) 


and, from equation (3), the parameter e is 


2lof VTo 1 
F 0 F, 


( 6 ) 


Now, we ask, what parameter e, expressed in terms of e, would make 
the first two terms of equation (2) exact: 


VTo l+ W 


(7) 


Combining equations (5) and (7), and using equation (6) gives 



( 8 ) 


or 



Thus, the square-wave example suggests equation (8) as a Euler 
transformation to improve convergence. 
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(2) Cone cylinder . — The ^-function, is a parabola, F—F 0 
Applying the foregoing procedure in this case gives 



( 9 ) 


Actually, equation (9) supports use of the simpler expression, equation 
(8), because the two agree to first order in e, agree at e= — <», and 
agree within 12 percent at e——I 0 /F 0 . 

(3) Sine-wave F-f unction, F=A sin vt/tq . — We choose this as a 
simple “smooth” case, and ask, what e would make the two leading 
terms of equation (2) exact in that case: 


we find, 


VL 




( 10 ) 


( 11 ) 


Equation (11) might serve as a suitable variable for smooth bodies. 

Now, if equation (8) is adopted, for discontinuous-slope bodies, 
the following formula results, using only two terms of equation (2) : 


VI F'o-Fl/2I 0 



Of course this formula has the feature previously discussed, that 
for the F-curve of figure 2, the numerator vanishes and the shock 
strength vanishes, except that when the N-wave is formed, numerator 
and denominator become equal, and I/I 0 is presumably one. 

Returning to the complete expansion, equation (2), we may rederive 
the series in powers of a variable 

e 

chosen in place of equation (8) for simplicity, and because r 0 is of unit 
order. This Euler transformation is evaluated on figure 3 for a 
particular cone cylinder for which the Whitham theory is in good 
agreement with experiment (ref. 3). It is seen that convergence of 
the transformed series is very good quite close to the body, even 
closer than «= — 1. Equation (12) is not displayed on figure 3, but 
its accuracy is comparable to the three-term Euler series, for the 
body chosen. 






Figure 3. — Comparison of expansion formulas with Whitham theory for cone cylinder. 
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CONCLUDING REMARKS 

It is clear that simple formulas can be derived by suitable choice of 
variable to give the shock strength quite accurately in the region of 
approach to far field and, as hoped, these formulas (e.g., eq. (12)) 
clearly suggest how changes in equivalent body shape can modify the 
signature in that region; the addition of a compression ahead of the 
expansion comer would improve matters, not everywhere in the flow, 
but as the N-wave is approached. The expansion method suggests 
that no modification is possible for smooth bodies which will favorably 
affect the nearly far-field shock strength. 

The approach described here should be extended to consider the 
signature shape behind the bow shock, hoping to find simple formulas 
for rise time, for example. The series upon which this extension 
would be based is 



Geometrically, E is the distance that the characteristic from r stands 
ahead of the front characteristic from r 0 . 


SUMMARY 


A series is derived for small departures from an N-wave for bodies 
of smooth and discontinuous slope. For the latter, the leading terms 
suggest that adding a compression just ahead of the shoulder reduces 
shock strength in the approach to far field. A Euler transformation 
is found to improve convergence of the series. 
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Nonlinear Effects on Sonic Boom Intensify 

M. T. Landahl, I. L. Ryhming, and L. Hilding 
The Aeronautical Research Institute oj Sweden 


INTRODUCTION 

The purpose of the investigation is to study those higher order 
effects on sonic boom strength that are normally neglected in the 
currently used first-order Whitham theory. In this exploratory phase, 
the method of matched asymptotic expansions is employed to assess 
the importance of second-order effects on the flow field in some simple 
cases. The associated change in shock strength can then be calculated 
from the disturbance flow field by use of the ordinary shock relations. 

TWO-DIMENSIONAL FLOW 

A complete account of the second-order complete flow field pattern 
for the two-dimensional case has been presented by Lighthill (ref. 1) . 
As an introduction to the more complicated axisymmetric case, the 
two-dimensional far field will be analyzed by continuing to second 
order the matched asymptotic expansion method developed to first 
order by Van Dyke (ref. 2). Accordingly, the perturbation velocity 
component u (£, ij; e) in the far field, where %=x—By, n=By, and e is 
a thickness parameter, is expanded as follows: 

«(€, r, e)=eUi^,H)+e 2 U 2 (i,IT)+ - - - H=e V ( 1 ) 

By inserting the expansion into the partial differential equation for 
the velocity potential and equating terms of like powers in e, one 
obtains 


U 1b +KUiUi ( = 0 (2) 

U, B + K(U l U 2 ) k =\ VlBB +M 2 (yUJJ lg +<p lB U 1( -l KUV^ (3) 
where K=±M\y+l)B~ 2 , and B 2 =M 2 - 1. 
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These equations can be solved most readily by employing the 
Monge theory (ref, 3). One finds that the solutions may be expressed 
parametrically as follows 


and 

where 


Ui—6(z) 
S=KHG(z) + z 


f7 2 =[iisT— (T+ 1) M 2 ] GP ( z)+Cr l 
t=B(l+KHG'(z)) 


( 4 ) 

( 5 ) 

( 6 ) 
( 7 ) 


The lines s(f, H)= constant describe the characteristics. Matching 
to the inner first-order solution gives 


1 


G(z)=~T'(z) 


( 8 ) 


where y=±tT gives the airfoil shape. For determination of the 
constant C, matching to the second-order inner solution is required. 

AXISYMMETRIC FLOW 

By employing a gage function analysis, it is found that the appro- 
priate expansion for u in the far field in this case reads 

u(S ) v;e)=e*U 1 (!t,H)+* ,> U a &H)+ . . . H=S V (9) 
The equations for U x and U% become 

Ui g +KUiUi ( +tjfjgU i=0 (10) 


< PlHB +M^U 1 U la +< P , a U H -\ KUiUr^ 


+ 


2 H 


<Pu 


(u) 


The solution for Ui may be written 


U i H ll2 =G{z) 

k=2KG{zW< 2 +z 


( 12 ) 


where G(z) is a fimction to be found by matching with the inner 
solution. According to Van Dyke (ref. 4), the inner solution to 
second order is given by 
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where 


0=«Vl + e 4 0>* 

1 r x ~ BT /(s l) da;i 

2irJ 0 V(*-*i) 2 -HV 

« 2 M 2 

<p?r=M}<px<pi x + Krtpi <pi x rv>f r 


(13) 

(14) 

(15) 


and e 2 /(x) is the source strength. The last term in the expression 
for must be included, because for r=0(t), <p lf is of order e _l , and 
this term is hence of the same order as the others in the neighborhood 
of the body surface. For r= 0(1), however, the term maybe neglected. 

To the particular solution for 02 one could also have added a 
homogenous solution of the same form as 01 . However, it is actually 
simpler to incorporate this term in 01 and instead assume / to have 
a correction term of order « 2 (or, of order e 2 In e for a slender body) . 
Thus, it is assumed that 

i=i (i6) 


For the matching, the outer limit of the near-field solution is taken 
by replacing ri=Br by H/e* in (13) and expanding in powers of «. 
This yields 


«=€W _1/2 i r (e)-2iiCF(f)F'(|)]+6 8 [--H- 3 / 2 F< 1 >a) 


+H- 


where 

and 


H4M 2 -H:)0F’(f) 2 +2F«>(S)F'(e))] 


+ 


F(0- 


J_ dt 

kJ2jo 


Fa ' w — 


(17) 

(18) 
(19) 


If a homogenous term had been included in the solution for 02 , the 
outer limit would have included a term of the form 


e 6 H- 1/2 F 2 (£> (20) 

which would have required, for matching, an additional term of order 
e 6 in the outer solution. Comparing the two expressions (12) and 
(17) for u, one finds that 

G{z)=F{z) (21) 


315-766 0—68 9 
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and the solution for Ui is thus given by 


U l =F(,z)H~ l i i 


* =2KF( Z )fl 1 «+s 
To lowest order in e for a slender body 

/=e-*S'=S' 


( 22 ) 


(23) 


where S(x) = irt 2 R 2 (x) is the cross-sectional area, and the solution (22) 
is found to be equivalent to that of Whitham (ref. 5). By inserting 
the expression for TJ\ into the equation for Z7 2 and solving as above, 
one obtains 

U 2 H 1/2 (1+2KF'H 1 ' 2 ) =K ^K-M 2 (t+|)] F 2 F' log H 

+ ( AW-2K ) F 2 +2F^F'^ H~ i ' 2 -F™H-'+ const (24) 


where F(z) and F lt) (z) are the same functions as before. To deter- 
mine the free constant, one needs to match to a third-order (of order 
e 6 ) inner solution not presently available. The other terms match 
completely with the second-order inner solution. 

To determine the relation between / and the body shape, we expand 
<Pi for small r in the manner of slender-body theory. Hereby we 
find after some manipulations that 

*H^+5K w »-£>+--- (25) 


Inserting this expression into the formula for <p 2 and carrying out the 
differentiation with respect to r, one obtains 


* r = 2^ (a!)+e4 


. m 2 a .. . 

, K 


(2v) 2 r JJ ^ 2 \2irJ r 2 4 ir J f 


(26) 


The tangency condition at the body surface requires that 


(Id ~<Pz)<p r —eR' at r=tR(x) 


(27) 
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By introducing the expansion (16) for /, and iterating by employ- 
ing the lowest order approximation /i in the second-order terms, one 
obtains 

f 2 =B 2 [~SS' "-(s^/J-M 2 [«3"+x* (S') 3 ]-g (S% (28) 

in which/— /i should be used in the calculation of <pi x . 

One finds that / 2 has the finite value of —M 2 <piS" at the rear end 
of the body. Hence, the modified source strength will have a dis- 
continuity similar to what would appear in the first-order source 
strength for a body with a corner. Because corners have not been 
allowed for in the above simple analysis, the calculation of the second- 
order correction to the flow field for regions influenced by the aft end 
of the body cannot yet be carried out. 

From /i and / 2 one may then determine the second-order Whitham 
function, F Wj 

F w = e 2 F= e 2 Fi + € 4 F 2 (29) 


by inserting j—ji + e 2 f 2 into equation (18). 

As an example, the first- and second-order source strengths and 
Whitham F-functions are computed for a parabolic body of revolu- 
tion of thickness ratio e— 0.07. As seen in figures 1 and 2, the second- 
order effects are fairly small near the forward portion of the body 
but quite noticeable toward the rear. Hence one would conclude 
that the higher order effects would have most influence on the rear 
shock. 

THREE-DIMENSIONAL FLOW 

Three-dimensional nonaxisymmetric flow is most conveniently 
treated in cylindrical coordinates x, r, and 0, The small parameter 
e in this case is assumed to be the thickness ratio of a body of revolu- 
tion having the same cross-sectional area distribution as the three- 
dimensional body. For the far field one then assumes 

u^U l (^H } d)+ e 8 U 2 (£,£T,0) (30) 

with £ and H defined as before. Upon introduction into the differ- 
ential equation for the velocity potential, one finds that the ilnear 
term r~ 2 <p ee will only appear in the equation for U 2 and that conse- 
quently, the equation for Ui will be the same as for the axisymmetric 
case. Hence, the solution for U x may be expressed as 

U l H il2 =F(z; e) 

£=2 KF(z) 0)H 1/2 +s 


(31) 

(32) 
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Matching to the inner first-order solution shows that F is given to 
lowest order by (18) as before, but with j{x; B) calculated from the 
equivalent body of revolution for each 6 as determined by the stand- 
ard supersonic area rule procedure. Thus the correctness of the 
commonly employed sonic boom calculation method for lifting wing- 
body combinations is verified to first order. The second-order 
solution is found to be 

U 2 H l ' 2 ( 1 +2 KF'H 11 *) = ^ KF 2 F Z [if-M 2 (r+|)] 

+2K 2 F z (^+^)j-log H 
-\^f 2 +f^fa z ^(~-m^2k(ff m +f^ 
+2KF z jF es 

-(i j> dz const (33) 


where subscripts 6 and z indicate partial derivatives. The solution 



Figure 1 . — First- and second-order source strength distribution for a parabolic 
body of revolution of thickness ratio c=0.07. Radius function, R— 2x(l — x), 
M— 2.5, 7=1.4. 
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Figure 2. — First- and second-order Whitham function for parabolic body of 

revolution in figure 1. 

differs from the one for axisymmetric flow by the terms containing 
derivatives with respect to 0, which all are generated by the term 
r~ 2 <pee. For the nonaxisyinmetric case no second-order near-field 
solution is readily available, and the matching cannot therefore be 
completed. Hence, the second-order correction to the Whitham 
function cannot be calculated at present. 

CONCLUSIONS 

It was found in the analysis that for three-dimensional flow the 
first-order equation describes the far-field wave structure very ac- 
curately with a relative error of only 0(c 4 ) (the second-order term is of 
order e 8 as compared to e 4 for the first-order term) . As a consequence 
the dominating higher order effect is that attributed to the nonlinear 
relationship between body slope and disturbance velocities at the 
outer boundary of the near field, which produces a second-order cor- 
rection to the Whitham F-function. A calculation of source strength 
change, and hence the correction to the Whitham function, for a 
slender axisymmetric body indicates that the second-order contribu- 
tion may be quite noticeable, particularly for regions toward the rear 
of the body. Also, the correction term increases in magnitude with 
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the Mach number approximately as M 2 , so that it becomes more 
important at high supersonic speeds. 

For three-dimensional configurations of general shape, no simple 
second-order solution for the near field is available, and one would 
then have to resort to some suitable approximate method. One such 
method, which has been employed recently (ref. 6) with considerable 
success for bodies of revolution with and without angle of attack, is 
the method of parametric differentiation (ref. 7) in its “local approxi- 
mation” and work to apply this method has been initiated. 

The possibility also exists to determine the near field from experi- 
ments. The theory indicates that by employing the measured 
disturbance flow field along a cylindrical surface at the edge of the 
near field, one should be able to predict very accurately the far-field 
flow pattern and shock strength. 

In the analysis for the axisymmetric body the slender-body approxi- 
mation was employed. To be valid the slender-body theory requires 
that the body is smooth without slope discontinuities. For a body 
with sharp corners there exists no simple relationship between body 
slope and source strength, and, in addition, the expansion for the far 
field will need to be modified. Another failure of the slender-body 
theory was that it produced a finite correction on the source strength 
and hence a jump in the source strength at the rear-pointed end of 
the body, thus precluding the calculation of the Whitham F-function 
for characteristics that meet the body axis downstream of the body. 
Consequently, the second-order effects on the rear shock cannot be 
analyzed until the slender-body approximation is removed. 
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The Feasibility of Larse Sonic Boom Reductions 

Adolf Busemann 
University of Colorado 


INTRODUCTION 

Small sonic boom reductions can be achieved by a large variety of 
changes in the airplane design until the optimization reaches the limits 
of absolutely rigid constraints; the possibility of finding new tricks 
to deal with the laws of nature may be remote in fields where our 
knowledge is well established and where the connection of all effects 
with the laws of nature is rather transparent. Supersonic flight and 
even supersonic aerodynamics do not have a long history as yet, but 
linearized theories in both two-dimensional and axially symmetric 
flow have existed for about 40 years and have lost their initial 
opaqueness with the years. Nevertheless, the central problem of this 
meeting, the sonic boom itself, came as a surprise with the beginning 
of supersonic flying some 20 years ago, though projectile waves and 
their investigations as Mach waves and shock waves count almost 
a century. Such historical facts may encourage our hopes for a 
compensatory surprise, or they may discourage us on account of the 
clarity of the supersonic theory. A linearized theory wins in time an 
apparent clarity, but can develop a widespread superstition as a 
result of its oversimplification. A nonexistence proof is still a dif- 
ferent matter from a generally accepted transparency. 

INVERSE RELATIONS 

The sonic boom, as the footprint on the ground of a supersonic 
airplane in steady flight, seems to create more trouble than the benefits 
of high-speed transport can outweigh. Under such conditions the 
permissible footprint is getting the status of the independent variable 
or “master/ 5 while the shape of the airplane becomes the dependent 
function or “slave.” Such reversal of a mathematical relation is the 
daily occurrence in engineering problems where “on second thought 55 
everything is interdependent; but one should never forget that the 
seemingly trivial inversions of the simplest arithmetical operations 
created great mathematical revolutions which brought our number 
system or number space to perfection. 
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Because the pressure distribution in the far field, according to 
linearized supersonic aerodynamics, consists of a half differential with 
respect to the distribution of the cross-sectional area combined with 
the lifting forces along the airplane length, the inverse relation is 
simply a half integral (actually it is a 3/2 differential in versed to a 3/2 
integral). As the name implies, the half integral is nothing but a 
half differential of an ordinary integral, which makes the mutual 
relations more or less equal. The trouble with noninteger differential 
or integral relations, however, is that the master function (in our case 
the airplane) may only be different from zero within a finite piece of 
the abscissa, while the slave function (the pressure in the far field) 
always gets a tail end stretching out to infinity. A complete reversal 
of such relations would attach the tail to the wrong partner. Impos- 
ing a tolerable footprint shape on an airplane of finite length must, 
therefore, lead to a compromise: prior to a given abscissa (say x=0) 
the footprint is the master and the airplane shape is the slave, but 
beyond that critical point the airplane is terminated and cannot play 
the slave anymore; thus the footprint falls back into its regular posi- 
tion as the slave with the tail. As long as infinite tales of decreasing 
strength are automatically unobjectionable, if started right, the 
proper precautions can be taken prior to the critical abscissa as a 
part of the prescribed footprint. 

The result of such considerations combined with the fact that 
linearized relations permit the superposition of particular solutions 
leads to the following compromise plan: 

Starting at a normalized abscissa prior to the terminal point (say 
at x= — 1 prior to x=0) a delta-function and integer powers of x+1 
are prepared for superpositions to a permissible footprint. For each 
particular function the half integral function F { {x) prior to the ter- 
minal point is established. Behind the terminal point Fi(x) continues 
as a Taylor (Maclaurin) series with integer powers of x limited to 1/2 
power below the power used in the starting portion of F { . The tail 
created by this convention is less objectionable than the particular 
function in the pressure component pt(x) prior to the switching point. 
Out of the integer powers of x the rear portion of the superimposed 
total function F{x) can be made to simulate the constant lift visible 
by the vortex wake behind the body, the zero cross section of the 
airplane body revised by a constant wake of cooled or heated air, 
and by positive or negative propulsion, if desired. 

CENTERED LIFT 

As soon as the F-function is determined, the area distribution of 
the body A(x) and the accumulated lift L(x) have to share the limited 
values of F along the abscissa. Because a resulting unbalanced 
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moment must be avoided, the center of mass and the center of lift 
have to coincide. This moment constraint can be seen clearly for the 
limit of a very soft body shell loaded with a constant specific weight 
throughout the body volume. In such a case the local lifting force 
l(x)~dL(z)fdx has to be proportional to the local cross-sectional 
area A(x) over the whole body length. The two relations, that 
A{x) plus a constant (depending upon altitude and Mach number) 
times the lift L{x) share the function F(x), while A(x) is proportional 
to L'(x), makes the dividing of F into its parts A and L unique. If 
L is plotted first and A is the difference left between F and Z, the 
slope of Z must be made proportional to F—L. A hatchet with a 
given length of its shaft guided with the end of the shaft along the 
incurve will cut automatically the Z-curve (displaced by the length of 
the shaft in the abscissa). This so-called “pursuit curve” is the first 
step to divide F between A and Z. Because the airplane shell is not 
really so soft as to require such a well-adjusted lift distribution, trading 
between local A’s and Z's can be made to improve the design of body 
and sweptwing. But the mass center and the lift center still have to 
coincide. Such constraint makes the footprint a genuine dictator for 
the complete airplane design. 


UNCONVENTIONAL SHAPES 

The considerations of the previous two sections of this paper are 
given in such detail because they can illustrate silent assumptions 
that have to be brought to light. Engineering compulsions and 
traditions very often introduce silent assumptions which we try to 
discover in “brainstorm” sessions. In case of severe restrictions, 
as they are imposed by sonic boom limitations, the conventional 
shape may not be able to survive. If, however, the details are exposed 
and the reasons for failing become obvious, the atmosphere for an 
unconventional exit is better prepared and the climate for a “break- 
through” may be created. The severe restrictions on the allocation of 
the given total F(x) for area and lift are a combination of the necessity 
of making lift center and mass center coincide, but under the silent 
assumption that wing and body are in the same plane. If we were 
allowed to raise the wing high above the body (fig. 1), the mass 
center and the lift center would still be in the same vertical line, but 
the wing abscissa and the body abscissa that create a combined 
footprint are inclined under the Mach angle. Disregarding the 
silent assumption that wing and body are coplanar makes it possible 
to interrupt the fighting between the body tip and the wing, which 
are additive in F(x ) , and allow the wing to cooperate with the receding 
tail end of the body of opposite sign at the latter part of F(x) by 
making fast exchanges without any pressure creation in the far field. 
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I center of lift 



The proposal to remove the coplanar wing and to add the elevated 
wing may be called the addition of a “quadrupole” by inspecting the 
change in the lifting vortices. But, not worrying about the name, 
the application speaks for itself as an unconventional shape to remove 
severe constraints. 


CONCLUSION 

The present abstract shows in detail only a small part of the 
feasibility study for large sonic boom reductions which is underway. 
The author is well aware that one day the full account has to be given 
of all the roads that were found blocked. But since optimism is a 
prerequisite for success, the few positive findings may still get priority. 
To mention just the general field explored, the linearized theory was 
investigated with respect to the stratified atmosphere in the gravity 
field of the earth including the temperature differences for the cutoff 
at the sides and the changes of Mach number along the way. The 
next higher order terms, neglected in the linearized approach, were 
explicitly determined and investigated as to meaning. The Mach 
number range was not kept as close to Mach number 1 as in the 
Brooklyn paper 1 of 1955, in which the reflected shock wave was 
allowed to create a long Mach stem with the incoming shock. Be- 
cause the pressure limits are in the order of one-tenth of a percent of 
the atmospheric pressure on the ground and because the altitudes 
reduce the pressures less than two orders of magnitude, the higher 
order terms are usually unimportant. Local superpositions of pres- 
sures high above the ground to speed up the decay without reaching 
the ground are hard to invent, but are not completely explored yet. 
Lower altitude flight with positive control over the finite slope in all 
pressure rises replacing the common N -shape of the footprint has 
advantages, if the ear responds much less, and if the airplane can 
sustain droplets and ice crystals or is allowed to change altitude on 
bad days. 

1 Busemann,\A<: Relation Between Minimizing Drag and Noise at Supersonic 
Speed. Paper presented at High Speed Aeronautics Conference (Polytechnic 
Institute of Brooklyn), Jan. 1965, pp, 133-144. 


Review of Second-Order Wave Structure 

David A. Catjghey and Wallace D. Hayes 
Princeton University 


The study of the wave systems associated with a body in supersonic 
flight is both of practical importance as it relates to the intensity of 
the sonic boom produced by a supersonic aircraft, and of theoretical 
interest as it illustrates techniques useful in solving problems of 
weakly nonlinear wave propagation. As has long been recognized, 
the classical linearized theories of supersonic aerodynamics do not 
present an accurate description of these wave systems at large dis- 
tances from the body. This is because cumulative nonlinear effects 
must be taken into account if we wish an accurate description of the 
wave after it has traveled a great distance. The nature of these 
effects has been considered by various writers (including Landau, 
Friedrichs, Whitham, and Lighthill), and Hayes has given the simi- 
larity variable governing the first-order wave structure (ref. 1). 

The present study is concerned with determining the second-order 
corrections to this wave structure. Thus, we attempt to find the 
second term in an asymptotic series representation of the solution, 
valid to all distances from the body. This is accomplished by the 
method of matched asymptotic expansions, wherein the solution 
valid at large distances from the body is matched to a local solution 
valid near the body which satisfies the actual boundary condition. 
In the large-distance, or wave-structure solution, account is taken of 
these cumulative nonlinear effects and shocks are inserted according 
to the appropriate conservation laws. Although the flow field is 
irrotational to the order of accuracy considered everywhere in the 
flow field (except behind the trailing shock), and a velocity potential 
is used, entropy jumps at the shocks are important and must be 
considered if the resulting shock shapes are to be correct to second 
order. 

The relationship of this theory to the other perturbation theories 
of high-speed aerodynamics is shown in figure 1. The first-order 
wave-structure theory is obtained by the inclusion of those nonlinear 
effects of second order which have a cumulative effect of first order. 

129 



130 


SONIC BOOM RESEARCH 


+( cumulative 
nonlinear 
effects) 


A 


First-order 
Wave Structure 

+ (higher order effects) 

Second-order 
Wave Structure 

<rM 4 TiT 3 ) 


-2 2-2 
1,M ,r*T3 ,u,v 


(Hayes, Whitham) 


(Caughey & Hayes) 


First-order 

i 

j 

Second-order 

Linearized 

+(higher order effects) 

Local 

(tB) 

i 

-2 2 -2 
i,m ,nre. 


(Ackeret) 


(Van Dyke) 


+(nonlinear N 

effects) 


First-order 

Transonic 

4 -3 

TM t.B 


(von Karraan) 


-Hhigher order effects) 


Second-order 

Transonic 

i*P*v 


(Hayes) 


Figure 1 .< — Relationships between perturbation theories of high speed 

aerodynamics. 

Similarly, third-order effects which exert a cumulative influence 
must be retained in the second-order wave-structure theory. The 
nature of the effects which exert a cumulative influence on the super- 
sonic wave structure is similar to that which must be considered 
locally in the theories of transonic flow. This is suggested by the 
likeness of the similarity parameters in the respective theories. 

The similitude for the second-order theory is not of the classical 
nature— i.e., it is not particularly useful for relating different flows 
while keeping a limited number of parameters fixed. Because the 
number of parameters obtained in the second-order theory is con- 
siderable, such a statement would be overly restrictive and of little 
or no practical value. Rather, a similitude like that of Van Dyke 
(ref. 2) for the local second-order theory is expressed: the important 
parameters of the problem are thought of as similarity coefficients, 
and are interpreted as giving the form of the second-order solution. 
Thus, the second-order potential is expressed as 

1 rM 4 


where T, and v are thermodynamic properties (defined in ref. 3) 
of the (general) fluid, evaluated in the free stream, M is the free- 
stream Mach number, and j3==^M 2 — T. The <Ai are functions 
invariant under the first-order similitude. Thus, the similitude given 
the form of the second-order solution and its explicit dependence 
upon these parameters. 

The precise manner in which cumulative higher order effects are 
retained in the theory is by contracting the coordinate which measures 
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distance along the wave system (which is very nearly along the 
direction of the free-strearn characteristics). The new contracted 
variable is found to be proportional to the so-called “age” variable 
in the analogous one-dimensional unsteady problems. These wave- 
structure variables are defined as 

1 = 1 — v 
YM*t — 

for planar flows and 

£=i— i 

2rM 4 r i=z 

n=- y" Vi? 

for flows about finite bodies, where (1,5?) are the Prandtl-Glauert 
coordinates in the plane of the free-stream direction and the azimuth 
angle of interest, with £ in the direction of the free-stream flow. 

For planar flows, the wave-structure solution is uniformly valid to 
the body surface, and the boundary conditions may be applied directly 
without the intervention of a local solution. Complete solutions 
have been obtained in this case, and confirm the results of Friedrichs 
(ref. 4) as presented and corrected by Lighthill (ref. 5), when 
specialized to their case of a calorically perfect gas (ref. 6). In 
addition to allowing treatment of a general fluid more easily, the 
present method is easily extended to the case of flows about finite 
bodies, where the cylindrical nature of the problem precludes easy 
application of Friedrichs' simple wave concept. 

For flows about finite bodies, the wave-structure solution must be 
matched with a local solution valid near the body surface. Here an 
intermediate solution (termed of 1 1/2 order) must be inserted to 
match with the outer asymptotic behavior of the local second-order 
solution satisfying the homogeneous equation. In addition, the outer 
asymptotic behavior of the third-order local solution is required to 
determine the second-order wave structure solution. Complete 
particular integrals have been found for the second-order solution, 
however, and the wave-structure solution may be easily constructed 
once the appropriate inner boundary conditions are determined from 
the local solutions. 

In the region behind the trailing shock, entropy and pressure 
disturbances are of the same order, and a velocity potential no longer 
exists. This region is important because the broad third-order 
pressure wave that is reflected from the shocks below the body even- 
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tually catches up with the rear shock above the body and alters its 
position. This effect must be considered if the theory is to be uni- 
formly accurate to second order. The appropriate scaling for this 
region is 


A rM 1 2 3 4 5 6 r 


A 

n 


=( 


r M 4 rV 

—)v 


in the case of planar flows, and 


II 

<>*AP 

vtm'tW 

< P* ) 5 

* i 


iH 

l P )' 


for flows about finite bodies. The rotational method of characteristics 
gives the strength of this wave, which is convected along lines parallel 
to the free-stream characteristics. In the region defined by the 
above scaling (which we term the Lighthill region, ref. 5), vorticity 
is again negligible, and the velocity potential exists to the required 
order. 
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Multipoles, Waveforms, and Atmospheric Effects 

A. R. George and A. R. Seebass 
Cornell University 


MULTIPOLES 

The azimuthal distribution of the pressure field about an aircraft in 
supersonic flight has been investigated by analyzing the various multi- 
pole contributions to this pressure field. In brief, it has been shown 
that for supersonic flows multipole effects are important in the far 
field, and that there are reasonably efficient ways of using multipole 
contributions to reduce sonic boom overpressures. 

The primary disturbances attributed to the flow about a supersonic 
aircraft propagate outward near a Mach cone originating at the air- 
craft. The pressure and momentum of the disturbance flow are con- 
centrated near this cone and in the wake. Thus, in contrast to sub- 
sonic flow where a disturbance effect is felt over a large ground area, 
a supersonic aircraft’s weight and volume effects are concentrated at 
the ground primarily between the intersections of the fore and aft 
Mach cones and the ground. 

There is no fundamental reason the motion of a closed volume at 
supersonic speeds need cause far-field flow disturbances and boom, as 
can be seen from the examples of the Busemann biplane or the analo- 
gous ring wing with central body. However, there appears to be no 
applicable way of eliminating the far-field disturbances attributed 
to lift without involving a net source or sink of energy or mass. How- 
ever, the boom attributed to volume and lift can be significantly re- 
duced by other methods. 

This section reports on a detailed analysis of the possibility of redis- 
tributing the variation of an aircraft’s pressure field around the Mach 
cone to reduce the overpressure directly below practical aircraft con- 
figurations (ref. 1). The utility of this idea is based on two factors 
affecting the propagation of disturbances in different azimuthal planes. 
First, disturbances in other than the vertical 0 plane will travel a 
longer distance before intercepting the ground and will thus have de- 
cayed somewhat more. If the dependence of portions of the flow field 
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is arranged properly, the maximum strength of the disturbances felt 
on the ground can be reduced. The second factor is the lateral cutoff 
of the “boom corridor,” which is attributed to the “total reflection” 
of waves traveling along rays inclined at more than the critical angle 
to the sound speed (temperature) gradient of the atmosphere. 

Previous analyses did not consider multipole effects explicitly 
because undue importance was ascribed to the fact that the disturb- 
ance from any slender configuration is essentially axisymmetric away 
from the axis. This need only be true for slender bodies. Nonslender 
configurations generally excite significant multipole contributions. 
An analysis of multipole solutions and the closure conditions that they 
must satisfy to be physically realizable leads to the conclusion that 
for any given azimuth angle a zero net dipole or quadrupole distribu- 
tion may be used to cancel the far-field disturbance of a closed non- 
lifting body. The boom attributed to a positive lift which corresponds 
to a positive net dipole distribution can only be completely canceled 
by canceling the lift or by a negative net equivalent source distribu- 
tion. (An example of such an approach is given by E. L. Resler, Jr., 
in “Reduction of Sonic Boom Attributed to Lift” in this volume.) 
However, the effects of lift can be reduced using multipole effects. 
An additional zero net lift distribution can be used to reduce boom 
but is generally not practical because of the moments introduced. On 
the other hand, additional quadrupole distributions can be practical. 
They can be looked upon as equivalent to zero net lift distributions 
where the lifting elements have been moved along Mach planes to give 
zero net moment. This also reduces the physical length of the aircraft 
compared to the length of the equivalent singularity distribution as 
discussed in reference 1. 

As was shown and discussed in more detail in reference 1, certain 
conditions must be met in order to excite significant multipole effects. 
The region where the boundary conditions are applied must satisfy 

Ax«Pr ( 1 ) 

where Ax is the axial extent of the boundary condition and r its radial 
coordinate. Using this concept one can see that only the initial 
leading edge region of the “ducted quadrupole” reported by H. W. 
Carlson in “Laboratory Sonic Boom Research and Prediction Tech- 
niques” in this volume would be expected to excite significant quad- 
rupole effect. Indeed the initial peaks of the Ap/p curves given for 
that configuration show a significant difference between the side and 
below and both of these initial peaks decay by the same factor of 
about 0.54 between the two values of h/l shown. The relative 
inefficiency of the tetrahedron model shown by Carlson is also easily 
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understood. Here the effect of the small portion of the leading edge 
which satisfies Az<</3 r is obliterated by the large volume of the 
model. A related test of a thin wing, having a much smaller volume 
(ref. 2), does show quadrupole effects clearly. 

For any given airplane configuration the equivalent multipole 
distributions are obtained by simply Fourier analyzing in 6 the far- 
field effective area or F-curve distributions. An aircraft generally has 
a source distribution corresponding to its volume and engine stream- 
tube area changes, a dipole distribution corresponding to its lift, and 
quadrupole and higher order multipole distributions corresponding to 
the off-axis positioning of various aircraft elements such as wing 
volume, pods, etc. For example, figure 1 shows the approximate 
source, dipole, and quadrupole contributions to the Whitham function 
for a version of the NASA SCAT 15F. As is discussed in detail in 
reference 1, additional quadrupole distributions can be used to reduce 
boom on lifting configurations. Also, because the wave drag can be 
related to the multipole distributions, boom reductions are made with 
consideration of wave drag changes. At present the possibilities for 
wave drag reduction are being studied. Examples of applications to 
boom reduction are given in reference 1. A sample calculation gave 
a 9.5-percent reduction in boom on the SCAT 15F at the expense of a 
14-percent increase in wave drag. Figure 2 shows the modifications 
made to achieve this reduction. 

WAVEFORMS 

In steady level flight an aircraft’s weight is ultimately transferred 
by the atmosphere to the ground in the form of a pressure field there. 
When the aircraft’s speed is greater than the speed of sound at the 
ground, then the main area of the ground available to the aircraft for 
its support is limited to the region between the intersection of the 
front and rear shock surfaces with the ground; the larger region behind 
the rear intersection contains the tail pressure wave left behind by 
these shock waves. Unfortunately, the primary region contains 
relatively large positive and negative contributions which, when they 
are combined with the negative tail wave pressures, nearly cancel: 
the difference is the pressure field that is required to support the air- 
craft. For radial distances that are large compared to a characteristic 
aircraft dimension, this pressure difference becomes negligible com- 
pared to the pressures that constitute the overall signature. Indeed, 
the term that must be present in order to provide a net force on the 
ground is proportional to r -5/4 , while the leading terms in the over- 
pressure signature are proportional to r“ 5/4 and r~ 7/8 . These terms 
provide no net force on the ground. The classical waveform formulas, 
which give a negative force on the ground, have been corrected so 
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that the leading terms yield, as they should, no net force on the ground 
(ref. 3). The further extension of this result to the order necessary 
to account for the net positive force on the ground is only of academic 
interest. Here we outline the derivation of the results; further discus- 
sion of their meaning and the difference between these results and the 
classical ones may be found in reference 3. It should be noted that 
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Figure 2. — Modification to the F-f unction for lower far-field boom. 

these results must be reinterpreted for an inhomogeneous atmosphere. 
The shape to which the signature for a homogeneous atmosphere has 
evolved in 7r/2 real atmospheric scale heights is the shape that will 
persist to any greater distance. In a real atmosphere the decay of the 
overpressure signature is not the r“ 3/4 law given here. For example, 
below the aircraft in an isothermal atmosphere the asymptotic decay 
law is (p/r) 1/2 . This point is discussed further in the submitted 
comments. 

For any finite body, the linearized pressure field at a given azimuthal 
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angle and radial distance can be integrated over the entire extent of 
the axial coordinate with the result 


J 'aeo 


Op (x, r, 6) dx=i2 2 sin Bjfir 


( 2 ) 


Here G v is the pressure coefficient and R\ sin $\ 1/2 is the ultimate base 
radius of the equivalent body of revolution, which is related to the 
lift by wR 2 =fi lift /p m U 2 . One way to arrive at this result is to con- 
sider the solution for the velocity potential, 


<p (x, r, 0)=— X) sin vB 

aTT n=0 


X 


■,-gr <r»(f ) cosh 1 (n cosh 


[(x—%) 2 —fi 2 r 2 ] 1/2 


d?, 


for various w-tuple distributions <x n (k), perform the indicated inte- 
gration and apply the closure conditions appropriate to the various 
distributions. If the result (2) is then integrated over the infinite 
lateral extent of the hyperbola formed by the intersection of the 
body’s fore Mach cone with a horizontal plane, then the lift is indeed 
equal to the force on the horizontal plane. After the attendant quasi- 
linear correction is applied to the linear theory, the resulting sonic 
boom waveform, when integrated over a horizontal plane, must again 
give a force that is equal to the lift : the quasi-linear correction simply 
relocates the pressure disturbance and introduces shock waves so 

/•CO 

that the value of the integral I C p (x, r, 0) dx is preserved. 

— QO 

Because we are interested in the leading terms in far field over- 
pressure signature, and as we have mentioned earlier, these dominate 
the right-hand side of equation (2) , we use the asymptotic form of the 
linear result 


/»« 

JflT 


G p (x, r, 0) dx-- 


-M 


F(y) d^=0 


( 3 ) 


Here we have introduced the ^-function of Whitham and the variable 


y = x —fir —kF(y) r 


1/2 


( 4 ) 


where k=2TM 4: /fi^j2fi; the parabolas y — constant are characteristics 
with their location corrected for quasi-linear effects. In those regions 
of the flow where characteristics of the same family intersect one 
another, the solution is now triple valued; the condition that determines 
the position of the shock wave that renders the solution unique is 
simply that 


J co 
-O 


G v das 
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be invariant. As a consequence the locations of the front and rear 
shock waves are given implicitly by 


f/F(t) &t=\kF\y t )r-™ 

F(t) df k[F 2 (y t ) — F 2 (y T ) ]r‘ /2 


where y t is the characteristic on which the rear shock terminates. 

To determine the asymptotic behavior of the overpressure signature, 
we note that as r— , y f and y T ^y 0 but that y t ~^ 00 ; thus we expand 
y r and y T in Taylor series about y 0i where F(y o )=0. It is easy to verify 
that 


and that 


'f.r=y 0 -(^r^^y r-^+Oir- 
^=(4^/4^) r 1/4 +0(r 1/8 ) 


. . dw 

The derivative ^ is obtained from equation (4) : 

dx R . RF(y) 
dy dr P ^~ 2 r 1/2 
dr - l-kF'(y)r 112 

Now at a shock wave we can write, to first order, 

( 7 ) 

where <r is the shock angle and A F is the jump in F{y) attributed to 
the jump in the flow deflection angle. Substituting the result into 
the expressions for dy/dr and carrying out the limiting process, we 
find: 

yr-Vo=I [^T7iJ7^-)]+0(r- I/2 ) 
y-Vo= -I \ympx ( y<) )+ 2/ 6 /V^ / (y 0 )]' 1 " 0(r_1/2) 

y t =i(r Ui - 7J ^)+od) 


{*J>) dtj* 


where 
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For simplicity we have tacitly assumed that F'(y 0 )^ 0 and that 
F(y) is continuous at y 0 . 

Finally, making use of the fact that 

f F(t) dt=0, 


we may reexpress (5) and (6) as 


and 


^cr-j>>*r 


The second integral in each formula is evaluated by expanding F{y) 
in a Taylor series about y 0 : 


The third integral in the second result is evaluated by using the 
asymptotic behavior of F(y) for large y. This gives 

kF(y,)=-I [l+jy.^ j r ,e,—psfs\ 

and 


my,)=~ 4 ^ 37 -s 

In an earlier derivation we expanded the square roots indicated in the 
first two formulas. In many instances the third term in the second 
formula is not sufficiently small for such an expansion to be valid. 
Consequently we shall leave the roots as they stand in these formulas. 
Furthermore, in computing the actual position x— &r of the shock 
waves, we combine y r with ~kF(y r )r 1/2 by first expanding the roots 
and then recasting the results back in to the form of a square root. 
This preserves the simple structure of the formulas and avoids using 
an expansion that is marginally valid. 
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The correct first-order results, which are given below, locate the 
rear shock wave so that the departure of the main portion of the 
signature from the final N-wave form compensates for the tail pressure 
wave to the extent that if l is a typical length for the configuration, 
then 


imi 

It is important to note that these corrections apply only to the 
asymptotic waveform formulas and follow directly from the basic 
theory as given by Whitham (ref. 4). 

The correct first-order results for the location of the front and 
rear shock waves, the pressure behind the front shock, the linear 
decrease with x to its value ahead of the rear shock wave and the 
pressure in the tail wave behind the rear shock follow directly from 
the results given above and are tabulated below. 

FAR-FIELD WAVEFORMS 

Front shock wave: 

ti n i/2 

1 ~W > (y.)rV*J r ' /4 

k VI Cp=I [ X + k J r " 3/4 

Between the shock waves: 


*-y/§ C p =—(x—Pr—y 0 ) £l- 


2 "I 1 / 

'kF'{y o y » J 




Rear shock wave: 

x—p T =y 0 +I^l 


kR 2 sin Of' 2 


kF'(y„)r U2 I 5/ V /8 


T 


- 1/4 


[i c /fi 

, 1 

kR 2 sin q 1 

/2 y L 

1 kF'(y„)r 1/2 

j5/2 r l/8 J 


-- 3/4 


Tail wave: 


/ fP o ~kB 2 sin o 
2 r™(x-pr-y 0 )™ 


where 


/= [ 2 *r >( ° d *i 1/2 


and k=2TM i /frl2j3 
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The NASA Langley computer program (ref. 5) has been used to 
determine the sonic boom waveforms in the azimuth planes below 
and above the B-70 aircraft (0=^ + x/2; 0== — 7r/2). Figures 3 and 
4 compare the computed pressure signatures given by the formulas 
above, and the uncorrected results. Excellent agreement between 
the far-field formulas and the actual waveform is obtained. 



Figure 3. — A comparison of the waveforms for the B-70 aircraft at a Mach 
number of M— 2.6, an altitude of 63 360 feet, a gross weight of 372 000 pounds, 
and (?~tt/2. 



Figure 4. — A comparison of the waveforms for the conditions used in figure 3, 

but with 6 = —71-/2. 
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ATMOSPHERIC EFFECTS 

Upon initiating an investigation of the turbulent dispersion of 
sonic boom signatures, it became apparent that the simple problem 
of propagation in a layered atmosphere was not universally under- 
stood. Results of various investigators were in disagreement. To 
clarify this situation a short study was made by Plotkin (related in a 
private communication) of the simplest possible case: the effect of a 
steady, horizontally stratified atmosphere on the overpressure directly 
below the airplane. Our analysis agrees with the results of Randall 
(ref. 6), The resulting correction factors were computed to better 
than 1 percent accuracy and compared with those of other investi- 
gators. We agree with Randall when we use his model atmosphere 
and essentially agree with the results of Kawamura and Makino 
(ref. 7) for a standard atmosphere. There are significant differences 
from some of the earlier published results (refs. 8 and 9). A com- 
parison of various results is shown in figure 5. 



Figure 5. — A comparison of the atmospheric correction factors reported by 

various authors. 
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Uniform Ray Theory Applied to Sonic Boom Problems 

M. B. Friedman and M. K. Myers 
Columbia University 


The subject of the present research is the propagation of shock 
waves from a slender body of revolution moving at supersonic speed 
through a stratified medium in which the ambient sound speed is a 
monotonically decreasing function of altitude. The problem is 
treated using a first-order Whitham correction to the linear field. 
The major difficulty arising in such a procedure is the development of 
a uniformly valid asymptotic expansion of the linear field which is 
appropriate as a basis for the Whitham technique. Such an expan- 
sion must take account of the interaction between diffraction effects 
attributed to the medium and the motion of the body and the primary 
disturbances. It must also account for the subsequent occurrence 
of focusing in the field attributed to refraction of the signal by the 
medium. 

The uniform linear field may be obtained in a parametric form 
using a ray technique which is a generalization of classical geometric 
acoustics. In recent years the subject of asymptotic solutions using 
rays has undergone wide development. However, the major portion 
of the existing theory of uniform expansions is concerned with prob- 
lems involving time-harmonic wave propagation. Thus, a major 
effort in the present work has been to derive corresponding gener- 
alizations of geometric acoustics to describe the propagation of 
arbitrary pulses. These uniform asymptotic representations serve 
to correct certain anomalies occurring in the simplest theory of 
geometric acoustics. In particular, the classical theory predicts in- 
finities in the field at points of focusing where the ray family forms an 
envelope, and also it predicts an infinite value for the field at points 
on the boundary between the region of the main signal and the 
diffracted field. From a study of certain exact solutions to typical 
problems of this type, it can be seen that such points in the field are 
located in regions where the linear solution varies most rapidly, and 
the anomalies above are a reflection of the failure of the geometric 
approximation to describe such regions. 
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These notions suggest that the problem of focusing is best analyzed 
by considering it as a stage in the transient development of the dis- 
turbance field. From this point of view, the process of developing 
an asymptotic theory for the treatment of a general problem involves 
several stages employing different asymptotic expansions. The first 
significant stage is governed by diffraction effects. A second is 
associated with the formation of ray envelopes in the field as a result 
of these interactions. Under certain circumstances, i.e., the presence 
of boundaries, the first stage of propagation is of major interest in 
itself. In any event, it provides the knowledge necessary to study 
the onset of focusing. 

It is characteristic of modem ray techniques associated with asymp- 
totic solutions to linear partial differential equations that certain 
functions appearing in the solution are left undetermined in the general 
formalism. A basic element of these theories is the hypothesis that 
these functions are determined by comparison of the asymptotic 
expansions with exact solutions of certain canonical problems. A 
canonical problem is one which is sufficiently simple to be solved 
exactly, and which is locally similar to the more general problem 
being considered. Consequently, for each stage in the development 
of asymptotic expansions for a general problem, it is necessary to 
identify and solve uniformly the appropriate canonical problem. 

UNIFORM RAY EXPANSION FOR STAGE ONE 

A canonical problem for the asymptotic development of the first 
stage is a two-dimensional point source set impulsively into uniform 
supersonic motion in a medium in which the sound speed is a constant; 
the source strength may be an arbitrary function of time. The 
geometry of the disturbance front is apparent and consists of a Mach 
envelope issuing from the source and an elementary initial disturbance 
front. The Mach envelope is tangent to the initial front. At the 
points of tangency, the potential <p(x, y ) t) of the point source has a 
directional singularity, and classical (nonuniform) ray theory predicts 
infinite values for the field at these points. 

A study of the exact solution has established that a uniformly 
valid ray expansion for <p exists of the form 

<p(x, /)=m [Cir) 1 '! s ”+2Xs b+i/2 (l) 

Here, s(X , t)= 0 represents the characteristic Mach envelope in space- 
time and s(X y t) — 0 represents the characteristic elementary cone of 
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the initial disturbance. 1 The quantities v n and v n represent the 
strength factors which propagate along the two sets of rays correspond- 
ing to the main rays (bicharacteristic of s=0) and the diffracted rays 
(bicharacteristics of s— 0). The canonical function F 0 is given by 



—sin*” 1 


W-&T 


the functions F„, n—l, 2, . . . are essentially successive integrals of 
Fo. 

The first terms of the series may be written explicitly : 

<P~v 0 H(sH { Vo s 1/2 -[% sin-{^) l/2 -(^=) 1/2 ] (2) 

H(s) is the unit step function. For ?<0, the leading term provides the 
usual first-order geometric acoustic approximation. As s— >0, the 
region of validity of the geometric approximation becomes vanish- 
ingly small. For a uniformly valid result, all three terms in equation 
(2) must be maintained. On the other hand, if only s is close to zero, 
an approximation can be obtained across the initial disturbance 
front, namely, 

<P~Vo(s) 1/2 

where v 0 — ^ QDSt „ : 0 O corresponds to the polar angle of the point 

of tangency. This is precisely the first-order geometric approximation 
across the initial disturbance front and becomes infinite along 0=0 O , 
the ray (shadow boundary) through the point of tangency. However, 
if all three terms are maintained to form the uniform expansion (2), 
then this expansion of <p is continuous across the shadow boundary and 
has the correct directional singularity at the point of tangency. 

This uniform expansion provides the information necessary to 
develop the asymptotic solution involving an arbitrary variation of 
sound speed. Consider a point source set into uniform supersonic 
motion in a medium in which the sound speed is described by a 
monotonically decreasing function of altitude c(y). The potential 
<p(Xjt) which describes the field attributed to the point source satisfies 
the equation 

<Pxx~\~<Pyy C 2 (y) ^ U ~ y) ( 3 ) 

where f(t)=0 for and <p(X ) 0)—<p t (X J 0)= 0. 

1 Because the canonical function F n are functions only of characteristic vari- 
ables, extension to similar problems involving more space dimensions is straight- 
forward. 
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The solution is postulated, to be of the form 


*=. S^ n [(^) 1/2 ]s”+ Si re s n+1/2 

where the functions s, s, v n , v n are to be determined. The wave 
equation (3) implies a set of differential equations to be satisfied by 
these unknown functions. The first few equations are 


Sz+Sy c 2 (?/) S * — ^ c 2 (y) S ^ ^ 


iPo)v^y c 2 (y) 1 ^ 


1 


(Z>o)sS*+ (Vo)ySy fS*— 0 


( 4 ) 

( 5 ) 


Equations (4) implies that 5 and s are characteristic surfaces of the 
wave equation (3) . In the usual manner, they can be solved by the 
method of characteristics which constructs the surfaces s and s out 
of the bicharacteristics (rays, in spacetime). Equations (5) are 
easily shown to be ordinary differential equations along the bichar- 
acteristics and are the transport equations which must be solved to 
determine the strength factors v 0 , v Q . 

NONLINEAR CORRECTION 

The availability of a uniform expansion permits the application of 
the Whitham method to many cases of interest. As an illustration, 
consider a two-dimensional, wedge-shaped airfoil set impulsively into 
uniform supersonic motion. In the case of arbitrary variations in 
sound speed, the corresponding nonlinear field is expressible in para- 
metric form based on computations which determine the ray system 
and the associated uniformly valid linear field. In the special case 
of constant sound speed explicit results can be obtained which illus- 
trate the importance of the uniform expansion as a basis for the 
nonlinear correction. 

By assuming a steady state in this problem, the Whitham method 
predicts a straight shock of constant strength standing ahead of the 
Mach lines and extending from the leading edge to infinity. If the 
full transient problem is considered, classical geometric acoustics 
can be applied and predicts constant strength along the Mach lines 
through the leading edge and infinite values for the pressure at the 
points of tangency of the Mach envelope and the initial disturbance 
front from the leading edge. Thus the Whitham approach will 
fail to predict correct shock pressures in the neighborhood of these 
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points. However, the uniform ray expansion predicts finite values 
for the linear solution at the points of tangency and the Whitham 
correction may be applied to this linear field. When this is done, 
it is found that the shock strength decreases continuously along the 
straight portion of the shock in the neighborhood of the diffraction 
boundary and along the curved portion of the shock associated with 
the diffraction front. 

The above results are valid until the onset of focusing. A canonical 
problem and the uniform ray theory for the next stage, which involves 
focusing, is presently under study. 
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The ARAP Sonic Boom Computer Program 

Wallace D. Hayes 
Princeton University 
and 

Rudolph C. Haefeli 

Aeronautical Research Associates of Princeton, Inc . 


Under the sponsorship of the NASA Langley Research Center, 
Aeronautical Research Associates of Princeton, Inc. (ARAP), has 
prepared a computer program for the calculation of sonic boom signa- 
tures on the ground. A NASA contractor's report entitled “Sonic 
Boom Propagation in a Stratified Atmosphere” (to be published) 
presents the program with instructions for its use, together with 
a review of the underlying theory. Here we briefly outline the com- 
puter program to indicate its theoretical basis and its special features. 

GENERAL DESCRIPTION 

The theory on which the computer program is based is an inte- 
grated “standard” theory (refs. 1 and 2), making only those simplify- 
ing assumptions which are made by all first-order theories on the 
subject. The atmosphere is assumed to be stratified with thermo- 
dynamic properties and horizontal winds functions of altitude alone. 
The signal leaving the aircraft is presumed known and must be speci- 
fied as input to the present program in terms of an i^-function. A 
discussion of the derivation of the ^-function may be found, for 
example, in reference 3, and it is dependent on Mach number, lift 
coefficient, and direction. This signal and the aircraft maneuver 
provide initial conditions for the wave propagation which is given 
by the theory of geometric acoustics. (See Hayes' paper, “Geometric 
Acoustics and Wave Theory,” in this volume.) The calculations 
provide ray trajectories for the signal, and also ray-tube areas to 
determine the strength of the signal. A modification of the signature 
is made with the help of an age variable to account for nonlinear 
distortion and the presence of shock waves. The atmospheric condi- 
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tions have an essential role in determining the age and degree of 
distortion. 1 

Particularly novel features of the computer program are that ray- 
tube areas are computed correctly according to geometric acoustics, 
with both arbitrary maneuvers of the aircraft and arbitrary stratifica- 
tion taken into account, and actual (midfield) signatures are computed 
without the common simplifying assumption of an N-wave (far-field 
signature). The calculations are made in the program for a set of 
selected rays originating at the aircraft at particular flight times and 
with specified initial azimuth angles (about the aircraft longitudinal 
axis). 

As a caustic is approached, the ray-tube area approaches zero. 
The theory appropriate near a caustic is more complicated than that 
underlying the computer program, and the program simply stops 
when a caustic is reached. If a ray traveling downward becomes 
horizontal, no attempt is made to follow the propagation further, 
and the propagation computation is stopped. 

PROGRAM INPUT DATA AND THEIR USE 
Stratified Atmosphere 

The temperature and pressure or density may be specified at a 
number of altitudes, together with the magnitude and direction of the 
horizontal wind. This input may be in English or metric units, and 
either geometric or geopotential altitude may be used. The program 
linearly interpolates these input tables. The input must be made 
consistent with hydrostatic equilibrium if such a condition is desired. 
The 1962 standard atmosphere is contained within the program and 
therefore may be specified as an alternate to the arbitrary atmosphere. 

Aircraft Maneuvers 

The aircraft Mach number, flight path direction, altitude, and 
bank angle must be specified for a succession of times during the 
flight of the aircraft. The aircraft maneuver equations then provide 
the initial location of the signal relative to the ground and also give 
values of the initial conditions and the time derivatives required by 
the ray-tube area equations. Two options are provided in the 
maneuver equations, depending on whether or not load factors are 
specified. In option 1, load factors are used, and time derivatives of 

1 For example, consider an aircraft in level, uniform flight without winds, com- 
paring signal propagation in a uniform atmosphere (constant density) with 
propagation in an exponential atmosphere (density varies exponentially with 
altitude). Then the signal shape traveling through the exponential atmosphere 
at an infinite distance below the aircraft is the same as the shape occurring in 
the uniform atmosphere only ir/2 scale heights (approximately 32 000 feet) below 
the aircraft (ref. 2) . 
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velocity, Mach angle, and two flight path angles, which are required 
by the ray-tube area equations, are calculated explicitly. In option 
2, these time derivatives are obtained using quadratic fits to tabular 
inputs, and the load factors are not needed. In the program, the 
differential equations of motion of the aircraft are integrated to give 
its location. 

F-Fufictions 

The ^-function for the aircraft is input for different values of the 
azimuth angle of the initial ray direction and as a function of distance 
along the aircraft (aircraft station). The F-f unction input for the 
azimuth angle closest to that under investigation is chosen for the 
then current calculation. 

Other Inputs 

The height of the ground above sea level and the choice of azimuth 
angles for which propagation is to be investigated need to be specified. 
The aircraft wing loading and length are also required. A ground 
reflection factor (which multiplies the pressure) may also be selected. 

PROGRAM OUTPUT INFORMATION 

1. For each chosen azimuth angle and aircraft time, the ray tra- 
jectory (including time) is given at selected altitude intervals, start- 
ing at the aircraft altitude and continuing to the ground intersection 
or to an intervening stop in the program. The values of ray-tube 
area and age variable are also listed at these altitudes. 

2. At each ground intersection, the signature in terms of pressure 
versus time is given without regard to shock waves. An auxiliary 
S-f unction versus time is also given which then provides the location 
of the shock waves in the actual pressure signature. 

3. If a stop occurs because of a horizontal ray, the information 
of item 2 is given for the signal at that point. 

4. Ground intersection and peak pressure information obtainable 
by interpolation from results from a number of rays is given. 

5. The input information is listed in the output for the purpose 
of a data check. 

SELECTED RESULTS 

Some results obtained with the sonic boom computer program are 
illustrated in this section only. These results were obtained using 
zero wind speed, although the program can give results for any strati- 
fied wind profiles. In figure 1 comparisons are shown for four ray 
traces from an aircraft flying at 20 000 feet above sea level to the 
ground level at 2300 feet (Edwards Air Force Base) . The location 
of the rays (wave front) at various times after leaving the aircraft 
is indicated on the traces. At Mach 1.4, the straight-line propaga- 
tion in a uniform atmosphere may be compared with the diffraction 
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Figure 1. — Ray traces for four flight conditions. 


in the 1962 standard atmosphere. The difference in ground inter- 
section locations is 1370 feet. At Mach 1.2 the rays travel con- 
siderably farther, of course. Ray traces are shown here for com- 
paring propagation from an aircraft in uniform level flight and in 
a shallow accelerating dive (n r =0.15 g, 7 — — 1 °) . The ground inter- 
section for the diving conditions occurs 1340 feet shorter and 1.1 
seconds earlier. 

Figure 2 shows an F-function for an F-104 aircraft as a function of 
aircraft station LfL A (L A —5 0 feet). This F-function has been used 
as an input to the sonic boom computer program with the above flight 
conditions. The complete pressure signature which results at the 
ground for propagation in a uniform atmosphere is shown in figure 3. 
The pressure increment p ambleat ) is plotted against a local 

time parameter | x which is proportional to the length of the signature. 
The shock structure has been drawn (heavy line) using information 



Figure 2. — An F-f unction for the F-104; Mach 1.4, <£=0, C^==0.06. 


2.0 
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Figure 3. — Computed pressure signature at ground and resulting shocks. 


on the location of the. shock provided by the S-curve of figure 4. 
This auxiliary curve is an integral of the signature, making use of 
area balancing (mentioned, for example, in ref. 3) to remove the 
inadmissible multivalued pressures. Specific £>-curve intersections, 
as indicated in figure 4, yield the shock locations. 

Figure 5 shows pressure signatures for three flight conditions using 
the 1962 standard atmosphere. The shock strengths and the lengths 
of the signatures are compared in table I. Here the results at Mach 
1.4 for uniform atmosphere and standard atmosphere are compared. 
The pressures in the uniform atmosphere (first column, values in 
parentheses) are multiplied by a factor of 1.50 commonly used to 
adjust for the standard atmosphere, and can be compared with the 
data for the standard atmosphere (second column). Although in 
this example the computed adjust pressure jumps would be the 
same in both cases, the lengths of the signatures are significantly 
shorter using the standard atmosphere solution (refer to previous 
footnote). 



Figure 4. — -S-curve for identifying shock locations. 




Table I . — Shock Strength and Length of Signature 


[Estimated precision of data: -2 ±'0.5 ft, Ap±0.03 psf] 




1.4 

2kT= 1.2 


Uniform 

Standard 


Standard atmosphere 


Parameter 

atmosphere 

atmosphere 







Uniform level flight 

Uniform 

Acceleration in level flight, 

Acceleration in dive, 




level flight 

nr— 0.15 g 

«t= 0.15 g, y— — 1° 

h, ft 

119. 9 

113. 8 

114. 2 

115. 4 

* 1. 05 % 

116.7 

‘ 2. 19 % 

l 2j ft 

131. 0 

128. 3 

128, 3 

128. 8 

0. 39% 

130. 2 

1.46% 

Api, psf 

1. 47, b 2. 20 

2. 17 

1. 99 

2. 27 

1. 41% 

2. 36 

1. 86% 

A pi, psf 

0. 60, b 0. 90 

0. 90 

0. S3 

0. 91 

0. 96% 

0. 94 

1. 33% 

AP3, psf 

I 60, b 2. 40 

2. 39 

1, 95 

2. 45 

2. 56 % 

2. 53 

2. 98% 




a Percent increase over value for uniform level flight. 

b Value of k A Ap where k A = 1.50 (commonly used factor to adjust to standard atmosphere). 
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Three sets of results for different maneuvers at Mach 1.2 are 
compared in table I: uniform level flight, accelerating level flight 
(0.15 g), and accelerating dive (0.15 g, 7 = — 1 °). These latter two 
examples represent very mild maneuvers. The resulting sonic boom 
strengths, as shown in the table, are about 1 to 2.5 percent larger with 
acceleration only, and about 1.3 to 3 percent larger in the shallow dive. 

These results have been presented to show some of the capabilities 
of the sonic boom program. It is well known that more severe over- 
pressures occur for other conditions: Further calculations are being 

planned to obtain results for various combinations of maneuvers, 
atmospheres, and wind conditions. 
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Geometric Acoustics and Wave Theory 

Wallace D. Hates 
Princeton University 


GEOMETRIC WAVE THEORY 

The purpose of this note is to present the theory of geometric 
acoustics as it emerges as a special case of the geometric theory of 
general linear wave propagation. The study of geometric theory 
of wave propagation of any type starts with a study of linear solutions 
in a uniform medium which are proportional to functions (generally 
sinusoidal) of a phase variable <t>=K-r—o)t, with r a distance variable 
in a suitable euclidean space and k a vector wave number. The study 
yields a relation w — fl(«) term a dispersion relation. If « is real 
when k is real, the waves are termed nondissipative. The solutions 
obtained are termed solutions for plane waves, the waves being planar 
in the r space. A typical example is that of gravity waves in a flat 
ocean of uniform depth, with the r space two-dimensional. 

In the general geometric theory for nondissipative waves, the 
strict conditions above are relaxed, and an asymptotic theory in a 
slowly varying nonuniform medium is sought for which the local 
solutions are very close to those obtained for plane waves, and « 
and k are considered large in some relative sense. The solutions are 
again proportional to functions (generally sinusoidal) of a phase 
variable 4>(r,t), and also to slowly varying amplitude functions. The 
frequency and wave number are defined by 


e 

II 

sd-i' 

(la) 

K=V<j} 

(lb) 

and are themselves functions of r and t, and satisfy 


1 7+^=0 

(2) 

plus the condition that Ak is symmetric. 
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The dispersion relation obtained for plane waves depends upon the 
medium, and when applied in the asymptotic theory to the quantities 
defined in equation (1) gives a dispersion relation 

o>=&(«c,r,6) (3) 

defined in an augmented space (#c,r,£) . The dispersion relation (3) 
applied to equation (2) gives a first-order partial differential equation 
for k alone. The method of characteristics gives then, in place of 
equation (2), the ordinary differential equation 

^+v r a=o (4) 

holding along characteristics defined by 


dr 
d t 




( 6 ) 


These characteristics are termed “rays.” The symmetry of A k is 
used in deriving equations (4) and (5). The quantity c is the group 
velocity (ref. 1), while the quantity k/u is the inverse phase velocity. 
The frequency « obeys 


do? 512 

d f~bt 


( 6 ) 


so that if the medium is steady, co is constant along rays. The phase 
variable obeys the relation 


K*C £0 = K< VK 12 — 0 


( 7 ) 


The rays may be parametrized by coordinates in a parameter space 
a of the same number of dimensions as the space r. In the obvious 
analog with fluid mechanics in which c becomes the particle velocity 
and the rays become particle paths, the parameter space a becomes 
a Lagrangian variable space. 

The slowly varying amplitude functions required for a quantitative 
solution may be obtained from separate equations. With the waves 
nondissipative this is best done through a conservation law, one in 
which volume integrals of an appropriate energy or wave action 
density are found to be conserved. In such a law the appropriate 
density times a measure V of an infinitesimal volume element is con- 
stant along rays. A convenient definition of V is as the determinant 
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The quantity V obeys the relation 

din V 


d* 




(9) 


along rays. 

NONDISPERSIVE WAVES 

A dispersion relation is nondispersive if it is of the form 

<o=KC*(n,r,0 (10) 

where #c— /cn and n is a unit vector. The group velocity is then 

c=n c n +V n c n (11) 


with V n a gradient (normal to n) in the unit sphere, while the inverse 
phase velocity is n/c„. Plots of c and n/c n are found to be dual in 
the sense that the procedure used to go from either plot to the other 
is the same. 

The most significant special property of nondispersive waves is that 
the functions of phase need not be sinusoidal but may be arbitrary. 
Anew phase variable which is a mono tonic function of the old may be 
introduced, if desired, and the variables c*> and k are far less significant 
than they are in the general case. If the medium is steady, it is 
convenient to define the phase so that co= — 1; the phase becomes 
then simply a time variable measured by a fixed observer, with a 
suitably defined zero point. 

Another special property is that the right-hand side of equation 
(7) is zero, so that the phase is constant along rays. This has as 
one consequence the result that the phase may be used as one 
component of the parameter space a, and a replaced by (a',<£). 

With a thus replaced, equation (8) may be rewritten 


where 


V — C n A n f 0) — ArJ K 



( 12 ) 

(13) 


is a measure of the area of a ray tube formed by rays for a given 
value of <t> as cut by surfaces of constant <p. If the right-hand side 
of equation (9) is divided into normal and tangential parts, we can 
identify 
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d In k 

~d T 


nvcn, 


(14 a) 


d InA - , c 

d t V ‘ C 


(14b) 


Here V t is the tangential gradient, normal to n. 

GEOMETRIC ACOUSTICS 

A study of linear inviscid acoustic theory following the approach 
outlined above leads to the conclusion that its geometric theory is 
of the nondissipative, nondispersive type. The perturbation velocity 
q and the perturbation pressure p' are related by 


np'=pa q (15) 

while the dispersion relation (10) is 


c n (n,r,t) =a(r ,t) +n * li(r,e) (16) 

where a is the speed of sound and u is the undisturbed fluid velocity 
(wind). Equations (4) and (6) are replaced by equation (13) and 

^-(V.u) -n (17) 

dln« 1 /da , 5u , ltA 

-ar _ B~(5( + 5 t’V (18) 

From equation (11) the group velocity is given by 

c=an+u (19) 


The quantity constant along rays is pq 2 c n 2 A n /w 2 a. The ray-tube 
area A n may be obtained by a quadrature, essentially that of equation 
(14b), after a differential equation for the wave-front curvature 
V*n has been solved. 1 If the medium is steady, so that co is constant 
along rays, the result reduces to the classical result of Blokhintsev 
that pq 2 c n 2 A n /a is constant along rays. 

PROPAGATION IN A STRATIFIED MEDIUM 

A stratified medium is one in which the dependence of 0 on r and 
t reduces to that in one Cartesian variable, here chosen to be z* We 


1 W. D. Hayes: The Energy Invariant for Geometric Acoustics in a Moving 

Medium. Phys. Fluids, vol. 11, 1968, to be published. 
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replace r by r'+sk and k by #c'+K z k. The general dispersion relation 
(3) takes the form 

^ <a=Sl(ic' ,k Z) z) (20) 


Equations (4) and (6) give 




( 21 ) 


A separate equation for k z may be given, but we may consider it 
given by equation (20) in terms of k! and w. 

We use the term Snell's law for a refraction law, once integrated, 
for wave propagation in a stratified medium. The result, equation 
(21) gives directly the general Snell's law that k and co are constant 
along rays. In particular, also, the horizontal component «'/« of the 
inverse phase velocity is constant along rays. 

It is convenient to replace t by z as a basic independent variable. 
The ray equations (5) then take the form 


d z c z 


(23) 


and can be integrated by quadratures to give the rays r' (z ) , t(z). 

When the wave propagation is nondispersive the dependence of 
K in equation (22) upon k ' and co reduces to dependence upon the 
horizontal component 


^r=N(a',«) (24) 

CO Cfi 

of the inverse phase velocity. The volume element V may be shown 
to be given by 

F=c 2 A/co (25) 

where A is a measure of the area of a ray tube formed by rays for a 
given value of 4> as cut by planes of constant z. The area A may 
conveniently be defined 

(26) 

To evaluate A we apply the operator Aa' to equation (22). The 
total derivative d/d z in equation (22) is a partial derivative in a 
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(a ' ,z) space, and commutes with the operator V a <. Thus we obtain 


^^=V a ,K=V a ,N-V N K (27) 

In this equation Va'N is constant, while with N given V N K is a func- 
tion of 2 alone. The solution of (27) is then 

V.*'=(V B ,r') 2 0 +V a .N. f V n K d 2 (28) 


The determinant then gives A according to equation (26) . 

In the acoustic case we take u to be horizontal and dependent 
only upon z. The horizontal vector K may be evaluated from equation 
( 22 ), 

H- a 2 N+(l-ti-N)ii (9 ^ 

a^il-u-Ny—afN 2 

We define the horizontal unit vector i and j so that N= M and 
u t j=u— iiu. We introduce the angle 9 such that 7i 2 =sin 6, n'= 
i cos 6, and note that c z ~a sin 9. The derivative 


u\ cos 3 9 \ I . 
a 2 sin 3 9 J ** 


(30) 


is obtainable from equation (29), and permits A to be calculated 
through equation (28) . The conserved Blokhintsev quantity is 
pqfCnCtA/a—N^pq^ sin 9 cos 9. 

The theory given here is equivalent to that used 2 in a recently 
developed computer program for calculating sonic boom pressure 
signatures. In such a calculation, results are first obtained using 
geometric acoustics and then modified for nonlinear effects. 
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Similarity Rules for Nonlinear Acoustic Propagation 
Through a Caustic 

Wallace D. Hayes 
Princeton University 


BASIC EQUATIONS AND SIMILITUDE 

A caustic in geometric acoustics is an envelope of rays, and also a 
locus of wave front cusps. A point on a caustic is one where a ray- 
tube area is zero, but corresponds to a lowest order type singularity 
for points of zero ray-tube area. A caustic has a characteristic scale 
R , which may conveniently be visualized as proportional to the in- 
verse of the relative curvature of a ray with respect to the caustic 
(both rays and caustics may be curved). This scale will be defined 
in terms of wave front shape. We are interested in acoustic signals 
of characteristic length Z, which may be the total length of a sonic 
boom signal, or which might be the spatial period of a periodic signal. 

When the ratio LjR is sufficiently small, a boundary-layer approach 
is valid (ref. 1). Separate inner and outer analyses are carried out 
and are matched. The outer analysis is one for geometric acoustics, 
in which weak nonlinear effects may be taken into account. The 
inner analysis, in a thin region including the caustic, is a local analysis 
with stretched coordinates. The thin region is interpreted as the 
boundary layer. Our purpose here is to examine how weak nonlinear 
effects are to be taken into account in this local analysis. 

One difficulty arises from the fact that rays are really trajectories 
in four-dimensional spacetime, with the caustic a hypersurface in 
this space. We use a particular moving coordinate system, one 
moving with velocity a Q n with respect to the fluid at a reference point 
on the caustic. Here a 0 is the speed of sound at this point. In this 
coordinate system the caustic will appear temporarily frozen. This 
procedure is analogous to that used with a general moving shock wave 
to reduce it to a temporarily stationary normal shock. 

Another difficulty arises from the fact that the motion is, in general, 
not irrotational. This may be attacked by defining a reduced 
velocity potential, but is most easily resolved by taking advantage 
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of the fact that the analysis is local. We assume that the character- 
istic scale of the vorticity in the disturbed flow is at least of order R. 
Then the vorticity is properly neglected in our local analysis, and a 
velocity potential may be used. 

We aline our x axis in the direction of the flow (— u), and let the 
velocity potential be a 0 ($+#). The linearized equation in a uniform 
sonic flow is 

% v + $zz — 2a 0 ~ 1 a 0 “ 2 $ «= 0 (1) 

This equation must be modified is several ways. The term in ® tt will 
turn out to be negligible, and we drop it here. A term in ® xx times a 
linear combination of x, y, and z must be added to give the caustic 
behavior. The coordinate system is alined so that the coefficient of z 
is zero, and the term added is —R~ l (y+bx)$ XXf where b is a dimen- 
sionless constant. Nonlinear effects are taken into account to lowest 

order with the term — 2r$ x $ xx , with F= 2 (t+ 1 ) & a perfect gas. 

The resulting equation is 


% y +^-[R- 1 (y+bx)+2T^ xx ^2a 0 - l ^t^O (2) 


The spatial coordinate system is illustrated in figure 1. 

We assume next that the characteristic time scale of unsteadiness 
is <r*af*R f and that the characteristic distance scale of z variations 
is dr l R y with c and d constants of order no greater than one. The 
small parameter X, defined by 

\==L l/z R~ 1/3 (3) 

is introduced, and a transformation of the independent variables is 
made, 

(4a) 

y=^L\'~ lf q~R\ 2 r] (4b) 

z—d^R? (4c) 

t = 0~^(%q~ 1 Rt (4d) 

Quantities of relative order X 2 will consistently be neglected. Equation 
(2) is thereby transformed into 

4>„+rf 2 X 4 ^r-(^+&XS+2rL- 1 X- 2 ^)^^2cX^ r ==0 (5) 

The errors in this equation are of relative order X 2 , and the s-derivative 
term in is here dropped. The term in equation (2) would 
appear here as of order c^X*. 
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For reasons that will be apparent when we discuss boundary condi- 
tions, the transformation 

$=eZA- 1/s ty (6a) 

(6b) 

is made. The quantity e is a measure of the strength of the signal. 
The quantity K, to be identified as a similarity parameter, is defined 

K=Te\-*' 2 (7) 

The basic equation (5) then becomes 

tf'w— (v + &X| + 22S^) — 2c\ypzr~ 0 (8) 


with relative errors that are of the order of X 2 . If b and c are of the order 
of X, or if they are of the order of one and X is negligibly small, the 
two terms in b\ and c\ may be dropped. We assume that this is the 
case, and obtain the nonlinear Tricomi equation 

(v + 2 ^ |{=0 ( 9 ) 

The quantity R has in essence been defined in terms of its appearance 
in equation (2). At a given instant, the characteristics of equation 
(2) are approximately 

| (y-\-bx<>) 3 —R(x—x<,) 2 (10) 

The quantity R thus characterizes the cusps in the wave fronts. In 
general, then, the wave fronts in the linear problem should be found; 
their shape then determines the parameter R. Such a definition is 
preferable to one in terms of curvatures of rays and caustics, because 
the wave fronts are definite geometric entities and are Galilean 
invariant. Rays and caustics are kinematic entities, and their 
geometric projections are not Galilean invariant. 

The boundary conditions must be chosen to correspond to the 
physical situation. The main boundary condition is one describing 
an incoming signal on the hyperbolic side of the caustic. On the 
elliptic side of the caustic we have simply the condition that the 
solution must approach zero. In a detailed analysis, the incoming 
signal should be described in terms of the theory of quasi-linear 
geometric acoustics, because this theory is consistent with equation 
(9) for large y or ??. Here we take the simpler course of describing 
the incoming signal in terms of linear geometric acoustics, consist- 
ently with the Tricomi equation. 


SI 5—766 0—68 12 
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The ray-tube area may be readily shown to be proportional to 
y li2 , so that the pressure according to linear geometric theory is 
proportional to y~ m . Thus we describe the incoming signal by 



for x negative and y sufficiently large, with / known (dependent also 
upon z and t or £ and r). The transformation equation (6) was 
chosen to be consistent with equation (11), and ghtes 

— ^=r) H/4 i(f+| V™) (12) 


as the boundary condition to be used with equation (9) . Equations 
(9) and (12) give us a problem for which the solution is of the form 


i (13) 

The function \f/ is a functional of /. The parameter K is a basic 
similarity parameter. 

With nonlinear effects taken into account there are more details 
in the boundary condition specification and the matching of inner 
to outer solution. There is little difference in principle, however, 
and the similitude result equation (13) will still hold with a slightly 
altered definition of/. 



a 0 * 


INCOMING 



Figure 1 . — Coordinate system for equation (2). 
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Figure 1 helps us interpret the constant ft. If b is nonzero, there is 
a normal relative motion of the caustic with respect to the medium. 
This is characterized by a normal velocity Ma 0j where M is a Mach 
number for this normal velocity given by 

M==6(l+6 2 )~ 1/2 (14) 


LINEAR SOLUTIONS 

In studying a nonlinear problem which reduces to a linear one in 
some limit, it is essential to understand the linear problem. In 
this case the linear problem is that for the Tricomi equation 

^—#*1=0 (15) 

Note that this equation is invariant with respect to ? differentiation. 
The most familiar special solutions are those in terms of Airy func- 
tions, sinusoidal in £. Solutions must be chosen that are exponentially 
decreasing rather than increasing for large — 17. With this condition, 
solutions for large rj have the form 


[exp i ir)+ exp i (k£— | w)J (16) 


The terms in give the interpretation that the reflected wave 
(second term) is the same as the incoming wave (first term) but with 
a phase advance of %tt. By decomposing a general incoming signal 
in terms of a Fourier integral we see that the reflected signal may have 
a much different shape. This phase shift leaves the power spectrum 
unchanged, but not the shape of the signal. A jump (linearized shock) 
reflects as a positive logarithmic singularity (ref. 2), for example. 

For our purposes we are interested in details of the solution in the 
caustic region, not just the outer asymptotic expansions. Though 
one attack would be through Fourier integral representations, a 
more direct and (in the author’s opinion) easier method is through 
self-similar solutions (ref. 3). It is possible to find a solution which 
gives a jump in the incoming signal and a logarithmic singularity in 
the reflected signal. This solution is a suitable linear combination of 


fi=y 


wj_ a I m 

Vl2, 12; 2; 4,7 


(17) 


and 


iiijD 


(18) 
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with particular choices of the possible analytic continuations across 
the characteristics. The solution has a — % power singularity in £ 
on the caustic. It may be interpreted as a solution for or as a 
solution for \f/ with a delta function in the incoming signal in terms of 

h- 

The solution corresponding to an arbitrary incoming signal may 
directly be represented as a superposition of solutions of the type 
discussed. If the incoming signal has a singularity of power n , 
there will arise on the caustic a singularity in £ of the power n~%. 
If the incoming signal in terms of / has no singularity of power % or 
worse, the solution on the caustic will be bounded. In this case, 
with K defined in equation (7) negligibly small, the similitude of 
equation (13) holds with K set equal to zero. Pressures, obtained 
from equation (6b), are proportional to eA“ 1/2 . 

GUIRAUD’S SIMILITUDE 

We consider finally the case in which K is very small but the func- 
tion/ has. jumps representing incoming shock waves. The nonlinear 
similarity result does not give more than peripheral information, and 
the linear theory gives infinite pressures both on the caustic and in 
the reflected wave. No matter how small K is, nonlinearity is 
important within a narrow region near the caustic including each 
shock and its reflection. The scale of this region in the £ direction 
we term L'. 

We repeat the analysis above with L replaced by L', and determine 
L' by equation (7) with the condition that K'—l. We obtain thereby 


x , =« 2 /«r 2/6 

(19a) 

L'=e m T m R 

(19b) 

The parameter t is defined here so that the jump in / as defined in 
equation (11) is one. The new reduced variables (for a shock whose 
linear trace passes through the origin) are related to the old ones by 


(20a) 

v[=Kr m n 

(20b) 


(20c) 


The pressure in the narrow region may be expressed 




( 21 ) 
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This similitude result is in accord with that found by Guiraud (ref. 4). 
Because the pressure predicted by the linear similitude is of the 
order of eX“ 1/2 , the pressure of equation (21) can predominate only 
if K 1/5 is very small. This rather severe condition may possibly be 
circumvented through additional analysis, treating the pressure field 
of equation (21) as a pressure difference relative to one given by some 
linear analysis. 
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General Remarks on Sonic Boom 

A. R. Seebass 
Cornell University 


I would like to begin these remarks by reiterating a comment that 
I made previously. We must not overlook the fact that advances 
in technology that result in an improved overall efficiency for SST 
aircraft, such as improvements in lift-to-drag and engine thrust-to- 
weight ratios, as well as in structural weight and specific fuel consump- 
tion, will provide direct reductions in the sonic boom intensities. 
Furthermore, attempts at sonic boom reduction cannot afford any 
penalties that result from a compromise in the overall efficiency of 
the aircraft if the reduction is to be a realizable one. Sonic boom 
reduction need not be at the expense of overall efficiency; indeed, the 
most direct route to sonic boom reduction is through the aircraft’s 
overall efficiency. 

During this meeting we have heard a number of papers that dis- 
cussed the utilization of aerodynamic techniques to reduce the 
“impact” of the sonic boom. We are all aware, thanks to Professor 
Busemann, that the volume contribution to the sonic boom can be 
eliminated in all azimuthal planes by utilizing interference effects; 
it can also be eliminated in a single azimuthal plane by a configuration 
that exactly balances the lift and volume contributions to the equiv- 
alent body shape in that plane. Such a configuration has no net lift. 
While the lift contribution can be reduced by the addition of suitable 
volume elements, there is no way of escaping the boom attributed 
to lift alone, if the Bernoulli constant of the flow is unchanged. For a 
given lift and volume there is a lower bound to the far-field overpres- 
sure within the framework of the above proviso. Professor Ferri has 
attempted to approach this lower bound with realistic aircraft con- 
figurations. Personally, I feel his results are encouraging in that 
they point the way for modest improvements in far-field overpressure 
levels. Professor George has indicated how quadrupole effects can 
also be used to approach this theoretical lower bound. 
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Several of the configurations we have seen in the past 2 days have 
achieved a reduction in the overpressure by lengthening the equivalent 
body in the azimuthal plane under the aircraft, without lengthening 
the actual aircraft. It is important to keep in mind that the theo- 
retical lower bound is for a given length of the equivalent body, and 
is proportional to the inverse fourth power of this length, not the 
actual length. Consequently, such reductions may accrue because 
the effective length of the aircraft is increased, not because the actual 
shape is near the optimum one. Indeed, reductions are sometimes 
obtained even though the new equivalent body shape is further from 
that prescribed by the lower bound than the original shape. 

As I mentioned earlier, we cannot escape the boom attributed to 
lift without modifying the Bernoulli constant of the flow* The 
equivalent slender body of revolution attributed to lift never closes 
in the azimuthal plane under the aircraft; it has a finite base area far 
downstream that is proportional to the lift. Aircraft engines, how- 
ever, modify the Bernoulli constant of the flow that passes through 
them, and it is not entirely unreasonable to think of aircraft in which 
a reduction in engine stream-tube area is used to compensate for the 
growth of the equivalent body of revolution attributed to lift. Hypo- 
thetically, it is possible to totally eliminate the boom attributed to 
lift in this way. Actually the second law of thermodynamics requires 
that infinitely far downstream the exit stream-tube area be larger 
than the entering one. This pluming of the engine exhaust, which 
takes many aircraft lengths to occur, will give rise to a weak pressure 
disturbance at the ground. A portion of this disturbance accounts 
for the support of the aircraft by the pressure field there. Because 
the disturbance is spread out over such a large longitudinal distance, 
it is reasonable to assume that no significant steepening of the dis- 
turbance into shock waves occurs. Professor Hosier's careful analysis 
of the gains we may achieve by this route shows how very difficult it is 
to obtain a reduction in the overpressure with realistic restrictions on 
the turbine inlet temperature. Total elimination of the exhaust 
stream-tube area in the present configurations would only reduce 
the boom attributed to lift by about 5 percent (assuming this elimina- 
tion is accomplished with no additional wave drag — which it cannot 
be). Nevertheless, this is an important area for further research 
because it gives us one more adjustable parameter with which to 
attempt to achieve overpressure reductions. 

I would like to discuss in more detail another aerodynamic means 
of reducing the impact of the sonic boom : the possibility of designing 
aircraft with their overpressure signatures modified to reduce human 
annoyance. Here I tread dangerous ground. Our understanding 
of- which features of the overpressure signature are the most undesir- 
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able ones is far from complete. If shock waves are absent from the 
signature and the compression rise times are on the order of hun- 
dredths of a second, then there is little acoustical energy present in 
the audible range. The most annoying feature of the sonic boom, as 
experienced outdoors, lies in the shock waves themselves. On the 
other hand, the low frequencies, which are inaudible outdoors, contain 
the majority of the acoustical energy and cause structures to vibrate 
and windows to rattle. As a consequence, an overpressure signature 
without shock waves may not be much less annoying indoors than 
its fully steepened counterpart. Let me leave these unanswered 
questions with the ad hoc assumption that one of the most annoying 
features of the sonic boom, at least as it is experienced outdoors, is 
attributable to the presence of shock waves in the overpressure 
signature, and turn to the technical and answerable question of 
whether or not it is possible to design aircraft with overpressure 
signatures that do not contain shock waves. 

Several years ago, Ed McLean of the NASA Langley Research 
Center noted from his calculations that it takes several hundred air- 
craft lengths for the overpressure signature to reach its asymptotic 
form. This led McLean to examine the requirements on aircraft 
length and weight, for a fixed altitude and Mach number, to achieve 
finite rise time overpressure signatures. The results of his most 
recent study were given during this meeting. These results may 
seem to be fairly pessimistic ones. We note that to achieve rise 
times on the order of 10 -2 seconds, instead of front and rear shock 
waves, requires that a 600 000-pound aircraft be nearly 1000 feet 
long. 

In his study McLean has assumed the only effect of the atmosphere 
is that it alters the overpressure level and not the signature shape. 
But Professor Hayes pointed out previously that this was not the 
case, and he reminded us at this meeting that the signature shape 
that exists x/2 real atmospheric heights from the aircraft in a homoge- 
neous atmosphere is, roughly speaking, the shape that persists to 
infinity in a real atmosphere. Because I disputed this point, and 
incorrectly so, I have taken it upon myself to correct McLean's 
results for atmospheric effects by calculating the advance, a, that 
would occur in isothermal atmosphere relative to that which would 
occur in a homogeneous atmosphere, a h} for two scale heights: one 
characteristic of the stratosphere (J9== 20 000 feet) and the other 
characteristic of the troposphere (H=30 000 feet). These advance 
ratios have then been used to calculate the aircraft length, Z, required 
for a finite rise time signature in an isothermal atmosphere, relative 
to the length, Z ft , computed by McLean for a homogeneous atmos- 
phere; the results are shown in figure 1. 



178 


SONIC BOOM RESEARCH 


H « 20,000 FEET 



McLean's analysis is also overly pessimistic from another point 
of view, in that it does not take advantage of the isentropic tail 
pressure wave that exists naturally in the signature to help eliminate 
the rear shock wave and reduce the aircraft length requirement. 
I have made a rough correction for this effect by employing the 
signature shown in figure 2, where I have taken the length of the 
aircraft to be 1.7X and completed the signature using the asymptotic 
form of the tail wave for the given lift. Figure 3 shows the aircraft 
length required as a function of altitude, for two Mach numbers 
and two weights. As you can see, at 61 000 feet there is about a 
30 percent reduction in the length required for a given weight when 
we include atmospheric and tail wave effects. Each contributes 



Figure 2. — Finite rise-time signature. 
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♦ NASA LANGLEY-600,000 LB UNIFORM ATMOSPHERE 

• .. .. 300,000 LB 

Figure 3. — Aircraft length for finite rise times. 

about half of the total decrease in length. For typical SST altitudes 
and weights, the lengths that are required for a Mach number of 2.7 
are about twice the length of the Boeing 2707. However, if we are 
willing to fly slower, and if we anticipate the usual improvements in 
specific fuel consumption, lift-to-drag ratio, structural weight and 
structural designs, then so far as I can judge, such aircraft are only 
one, or perhaps two, aircraft generations beyond the present SST 
design. 

In conclusion, it seems to me that we can safely prognosticate a 
continual evolution of SST designs with improved sonic boom char- 
acteristics. Major gains may be expected from improvements in the 
overall efficiency of such aircraft, as well as through novel design 
features such as those just discussed. Whether or not these evolu- 
tionary gains will be sufficient to make a domestic SST an economi- 
cally viable concept is unknown. We cannot hope to provide a 
definitive answer until we know what features of the overpressure 
signature are the most annoying ones, and what integrated over- 
pressure loadings are likely to prove to be acceptable. My prognosis, 
based on our tentative answers to these questions, is that such gains 
will indeed be sufficient to allow commercial supersonic flight over 
populated areas. 
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State of the Art of Sonic Boom Theory 

Wallace D. Hayes 
Princeton University 


This paper could be divided roughly into two main questions. The 
first is, How far have we proceeded in developing methods of analysis 
and the requisite understanding of sonic boom phenomena? The 
second is, How far have we proceeded in devising methods of reducing 
sonic boom for feasible aircraft designs? 

ANALYSIS OF SONIC BOOM 

With respect to the straightforward analysis of the propagation 
of sonic boom from a slender aircraft, there appears to be no problem 
now. Linearized theory for the flow about the aircraft has been 
understood for a long time. Geometric acoustics with winds and 
variable sound speed is now quite well understood. With the assump- 
tion that nonlinearity does not affect the rays, the nonlinear modifica- 
tions of the signature can be calculated directly. Developed second- 
order theory serves to justify the theory. 

One problem of practical importance, though better understood, 
remains unsolved. This is the problem of predicting sonic boom 
signatures near a caustic (ray envelope) . This problem, together with 
other theoretical problems as yet unsolved, is at least in focus. We 
know how they are to be attacked, and there is no reason to believe 
that practicable solutions will not emerge with time and labor. The 
other problems include that of nonlinear effects on ray tubes, diffrac- 
tion into shadow zones, and nonlinear effects near the aircraft. 

The effects of turbulence, in particular the creation of spikes in 
sonic boom signatures, are beginning to be understood (see the remark 
entitled “Atmospheric Turbulence” by Garrick in this volume)* 
More work is clearly needed in this field. 

REDUCTION OF SONIC BOOM 

It is now generally recognized that the one inescapable parameter 
is the lift of the aircraft plus pU 2 fi~ l times the increase of engine jet 
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area over capture area (discussed in NASA SP-147, p. 6). It is also 
recognized that general overall excellence in aircraft design is most 
important in reducing sonic boom. The first fact is important in that 
it gives the designer criteria for evaluating his design with respect 
to sonic boom production, and permits formulating an ideal toward 
which to aim. The second fact is important in establishing that 
sonic boom is but another element that must fit into the procedure 
of compromise that is aircraft design, that a preoccupation in design 
with sonic boom at the cost of other design factors is self-defeating. 

To illustrate these points, with the lift plus engine area increase 
given, we can clearly reduce the sonic boom intensities from the 
volume distribution by suitable redesign. This procedure may, and 
generally does, entail some increase in wave or frictional drag. This 
drag increase would affect the overall design and require a larger gross 
weight and lift for a given mission. 

How should we proceed in design studies to reduce sonic boom? 
Here we are faced with two completely different paths. The first 
path is somewhat better charted. We accept a signature with shock 
waves and assume that the signature is close to an N-wave. With 
the lift given and the drag estimated, we can design engines with 
small area increase. With the parameter lift plus engine area in- 
crease specified, we can minimize the boom beneath the aircraft. We 
are led to a family of problems in the calculus of variations, some of 
which have been solved. We then seek aircraft configurations that 
will come close to realizing our optimum solutions. 

The second path is quite different. We adopt the psychoacoustic 
principle that it is the shocks in sonic boom signatures which are 
objectionable and are to be avoided. We do not assume an N-wave 
shape, and aim toward a finite rise time signature near the aircraft 
that will still have no shock waves in it at the ground. Here we 
take advantage of the “freezing” effect of atmospheric stratification 
on downward-propagating signals. The calculus of variations has 
not yet been applied in an attack along this path. The main effort, 
of course, is again in detail design, with very low wing loadings 
required. 

Which of these two paths toward sonic boom reduction will be 
best is not known now. They must both be explored. 



Reduction of Peak Overpressure by Configuration Design 


Antonio Ferri 
New York University 


In connection with the comments and suggestions for future work, 
I believe the following points can be made: It appears that by some 
redesign of present configurations substantial changes in the peak 
overpressure can be obtained. In addition, it appears that near- 
field signature can be obtained if the total length of the airplane is 
utilized even if we fly at 60 000 feet. On this basis, a substantial 
effort must be made to investigate these two possibilities. 

In connection with improvement of the existing theory related 
to sonic boom, my feeling is that the improvement must be obtained 
in two directions. First, we must obtain more realistic information 
of the initial distribution for a practical vehicle at Mach numbers 
on the order of 3 or higher. It is my impression from the presenta- 
tions given that many of the higher order theories are based on expan- 
sion in terms of Mach number and similar to the expansion used by 
Busemanri for the two-dimensional wing. Such expansions are not 
valid at high Mach numbers. If this is correct, the higher order 
theories are not very useful. In addition, the present theory does 
not permit analyzing wings of two-dimensional type or close to 
two-dimensional type. Such wings permit producing a lot of lift 
in a short length and therefore are useful for sonic boom considera- 
tions because it permits obtaining distribution close to optimum. 

If you neglect tip effects by cutting the tips along the Mach wave, 
then the L/D of two-dimensional wing is proportional to the inverse 
of the angle of attack because the lift is proportional to a and C D is 
proportional to a 2 . Therefore, if we increase the span and decrease 
the cord leaving the area the same, and angle of attack the same, 
the L/D does not change. 

In connection with propagation theories, the effect of the density 
variations as a function of the altitude 2 has been considered already 
by several people. Therefore, I do not see much necessity of in- 
creasing activities in this field. The part that has not been considered 
is the variation of density in a given plane as a function of x , y , and 
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not z. This we are planning to do for the next year; I believe that 
such effects will justify some of the variations obtained for similar 
flight conditions attributed to atmospheric conditions. 

In addition, the work that Professor Hayes recommended in con- 
nection with maneuvering vehicles is also important. In my opinion 
the work that we have initiated connected with the effects of terrain 
configurations is very important. I feel that the linear theory is 
more than satisfactory if we are interested in pressure variations on 
surfaces from the stress analysis point of view. In any case, the 
higher order theory will not improve the situation because it can 
be done only for a far field and the reflections at the surface of the 
body depend only on the near field. Again, we plan to continue 
this work next year. 



Notes on the Sonic Boom Minimization Problem 

Harry W. Carlson 
Langley Research Center , NASA 


In view of the widespread current interest in sonic boom suppres- 
sion, it is believed to be appropriate to review some observations on 
the prospects for sonic boom reductions based on analytic studies 
and the supporting experimental research program conducted at the 
Langley Research Center. An extensive wind-tunnel test program 
beginning about 10 years ago has treated a great variety of models 
including rather basic research shapes, airplane components, and 
complete airplanes. In addition, the tunnel has provided data for 
evaluation of a number of exotic or nonconventional airplane concepts 
and for evaluation of certain minimum boom and minimum drag 
shapes as well. Recurrent conclusions in all of these studies have 
emphasized the remarkable ability of the simple linearized theory 
(with the Whitham corrections) to provide accurate predictions of 
sonic boom signatures. It has been noted that the prediction ability 
increases with increasing model slenderness and with increasing dis- 
tance from the model. Thus the theory is recognized to be all the 
more valuable as a minimization tool because the conditions of 
slenderness and extreme altitude are precisely those which encourage 
minimization. 

An examination of the minimization problem by the author in 
1962 led to the determination of an approximate far-held lower bound 
and to suggestions for achievement of this goal through airplane design 
(NASA TN D-1494). The equivalent body principle formulated by 
Hayes and the work of Whitham formed the theoretical basis for the 
study. Because the sonic boom had been shown to depend on the 
shape and magnitude of an airplane effective area distribution formed 
by both volume and lift components (fig. 1), and, because the magni- 
tude of the effective area at the base of the airplane is fixed in large 
part by the airplane weight, Mach number, and altitude, a study was 
made of the effect of the shape of the curve or the bow shock pressure 
rise in the far field. 
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It was found that the overpressure could be expressed as 

Ap=p (£)%'»K T -K,y[l C L f 2 +ji ( 1 ) 

where k s is a shape factor that may be evaluated by employment of 
the machine computing program described in NASA Report R-213 
or by analytic methods for curves that can be expressed as simple 
mathematical functions. The results of the aforementioned NASA 
study of many effective area distribution shapes indicated an approxi- 
mate lower bound value of 2£ s =0.60. It was then discovered that 
L. B. Jones had earlier made a more rigorous mathematical study of 
the problem (for volume alone or lift alone, but not the combination) 
and had arrived at factors 10 percent lower. Jones’s lower bound 
shape, however, requires an area development proportional to the 
square root of the distance from the airplane nose and, as a conse- 
quence, has a prohibitively high wave drag for practical consideration. 
The high drag is the result of the relatively intense shock strength 
near the blunt nose of the body. Only at extreme distances does the 
blunt shape have lower bow shock pressures than more conventional 
shapes. The k s value of 0.60 thus may be accepted as a more practi- 
cal lower bound. The corresponding effective area distribution shown 
in figure 2 is in fact the area development for a Yon Kdrmdn minimum 
drag body. 

Equation (1) may be plotted as shown in figure 3. The ordinate 
represents an overpressure parameter which contains altitude, Mach 
number, and reflection factors. The abscissa is dependent primarily 
on the airplane lift. The relationship between the two parameters 
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Figure 2. — Effective area distribution corresponding to far-field practical lower 

bound overpressure. 
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Figure 3. — Far-field practical lower bound overpressure in parametric form. 


is an indication of the sonic boom efficiency of the design. When a 
shape factor of 0.60 is chosen, the curve as shown in figure 3 defines 
the practical lower bound of attainable overpressure levels. In other 
words, the area under the curve represents overpressure levels that 
are unattainable in the practical sense. In this case, the base area 
A b was chosen to represent the engine exhaust effect for an SST at 
M— 2.7 cruise altitude. There is one further restraint that must be 
applied. The area under the lower bound effective area development 
may not contain enough volume for a practical airplane at lower 


Q £1 

values of the lift parameter ~ C h - and thus the lower bound curve 

2 


must be modified in that region as was done in TN-1494. 

Figure 3 and equation (I) form a basis for a study of the means 
which may be utilized in attempts at minimization. Airplane shaping 
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considerations which allow an approach to lower bound overpressures 
and which lead to sonic boom minimization will be discussed sub- 
sequently; first, the dependence of the boom on the operational 
factors of altitude and Mach number and on the gross design features 
of airplane length and weight will be treated. 

Figure 4 illustrates the lower bound pressure dependence on the 
primary factors of altitude, Mach number, airplane length, and air- 
plane weight. The circular symbol indicates conditions for SST 
with a length of 300 feet and a weight of 600 000 pounds cruising at 
a Mach number of 2.7 and an altitude of 60 000 feet. The correspond- 
ing A p is 1.83 psf. 

As is well known, altitude is a powerful factor in the determination 
of overpressure not only because of the three-fourths power of the 
nondimensionalized distance h/l in the overpressure parameter, but 
also because the reference pressure p which must represent a mean 
pressure between the airplane and the ground, decreases rapidly with 
increasing altitude. As altitude increases, however, the lift parameter 
also increases, and the benefits of increased altitude become less. 
There is, in fact, little benefit to be derived from operation of the 
airplane at altitudes much above the altitude which gives maximum 
lift-drag ratio or maximum range. For the same range, an airplane 
flying at the higher altitude requires a larger fuel load which com- 
pensates for the increased altitude. 

Increased Mach number is seen to have a beneficial effect beyond 
about M—2. Furthermore, increased speed will provide benefits 
beyond the pure Mach number effect because the optimum altitude 





Figure 4.— Variation of lower bound overpressure with the various factors. 
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for airplane range is increased as speed is increased. The effects of 
the first two curves of figure 4 may thus be combined. However, 
current designs for hypersonic speeds tend to have large volumes for 
fuel storage, and good shape factors approaching 0.60 may be difficult 
to achieve. 

An increase in airplane effective length is also seen to have a bene- 
ficial effect. Because the effective length is determined by supersonic 
area rule concepts, the effective length can be increased by employ- 
ment of wing dihedral or other vertical displacement of airplane 
components from a flat surface. However, attempts to elongate the 
airplane effective length for sonic boom benefits will be limited by 
consideration of the effects on the aerodynamic efficiency and on the 
airplane weight. 

It is clear that airplane weight reductions will be of utmost im- 
portance in sonic boom minimization. Weight reductions have a 
direct effect on sonic boom reduction and do not necessitate considera- 
tion of opposing factors as is the case for other minimization approaches. 
Actually the benefits of a given weight saving will tend to be amplified 
as the airplane design cycle is retraced. Unfortunately, the converse 
is also true, and history has shown far more weight increases than 
weight reductions in the path from concept to production airplane. 

It can be concluded that maximization of the product of altitude 
and Mach number consistent with the maintenance of aerodynamic 
and propulsive efficiency will be instrumental in sonic boom reduction. 
It is also apparent that minimization of the ratio of weight to length 
(or lift loading per unit length) consistent with aerodynamic efficiency 
considerations will lead to sonic boom reductions. 

Airplane shaping considerations will be discussed with the aid of 
figure 5. Superimposed on the lower bound plot taken from figure 3 
are configurations which permit a good approach to the lower bound 
effective area distribution for two values of the lift parameter. Note 
that for lower values of the lift parameter, near the volume restraint 
curve, a relatively short chord wing, well aft on the airplane, can be 
combined with the volume to give a near optimum effective area 
distribution. In fact, if the lift is not concentrated near the rear of 
the configuration, it will be all the more difficult to provide the neces- 
sary airplane volume within the constraints of the lower bound shape. 
At higher values of the lift parameter, however, it is necessary that 
the lift be distributed over most of the airplane length in order to 
avoid a rather bulbous drag-producing fuselage. A highly swept 
arrow wing with considerable cutout helps in meeting the require- 
ments. Present SST designs at cruise conditions would fall some- 
where between the lift parameters for the two examples. 

Another consideration worthy of note is that although the volume 
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Figure 5. — -An illustration of configuration requirements for lower bound 

overpressure. 


contribution to the pressure field could be eliminated by employment 
of concepts related to the Busemann biplane, there is not necessarily 
an advantage in so doing. Volume effects can be combined with lift 
effects in a favorable fashion to produce lower overpressures than 
would be the case for the lift alone. In both of the examples shown 
in figure 5, the wing lift is generated in an expansion field created by 
the fuselage to bring about a favorable lift-volume interaction. 

In this discussion, no account has been taken of the potential for 
near-field signature benefits such as the flat-top or finite rise time 
signatures discussed by McLean in another part of this document. 
The minimization concepts discussed herein do, however, apply to 
the impulse of the signature for all flight conditions, and from present 
indications apply to Ap for SST cruising conditions. 

It is the author’s contention that practical minimization concepts 
are well understood and are supported by a large body of experimental 
data. Furthermore, it is believed that the surviving SST designs in 
the national program have utilized this technology to nearly the full 
practical extent. Cruise overpressure levels are only 10 to 15 percent 
higher than the approximate lower bound used herein and 20 to 30 
percent higher than the unattainable Jones lower bound for airplanes 
of that size and weight. 


Atmospheric Turbulence 

I. Edward Garrick 
Langley Research Center , NASA 


Professor Hayes raised the discussion question on likely further 
work on the interaction of atmospheric turbulence and the sonic 
boom. In response, Mr. I. E. Garrick stated that the question was 
obviously very broad and the problem areas offered much room for 
research, for example, in both topological and weather effects. 
Although the present sonic boom theories seem to predict very well 
the mean values of the statistics of sonic boom pressures, it should 
be feasible also to predict the variances of these statistics by including 
meteorological-turbulence parameters. Moreover, because it is the 
highest values, which occur least often, that may be of special concern, 
the prediction of the extreme value statistics deserves additional 
attention and effort. 



Brief Comments on Sonic Boom Reduction 


A. R. George 
Cornell University 


At the risk of stating the obvious, it appears that several boom 
reduction methods were reported that could reduce boom significantly, 
particularly cumulatively. Unfortunately the likely performance 
and/or structural weight penalties that would be associated with 
application of these methods would probably be intolerable for the 
present designs. However, in design studies for any future domestic 
SST, these various methods should certainly be considered and 
evaluated. 
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General Remarks on Sonic Boom 

Adolf Busemann 
Colorado University 

If the psychological aspect of the sonic boom indicates that a reason- 
able finite rise time of the positive pressure steps (rather than a step 
function) is crucial for the acceptability of booms by man and animals, 
then increasing the flying altitude beyond 25 000 feet does not help 
alleviate the sonic boom effects. At ground level the absolute pres- 
sure jump attributed to lift is not decreased by changing the altitude 
from 25 000 feet to 70 000 feet, but the width of the audible path is 
doubled. The self-propagation of a pressure equal to 1/1000 of the 
atmospheric pressure on the ground (1/2000 prior to the doubling by 
reflection) corresponds to about 140 feet at 25 000 feet altitude and 
400 feet at 70 000 feet altitude (for M— 3). This is in addition to the 
original distance of the positive and negative pressure peaks. Spread- 
ing the near-field pressure rise created by the nose and the lifting 
portions of the aircraft over such a length to avoid forming the steep 
pressure rises of the common N-wave in the far field severely increases 
the minimum length of the airplane with increasing altitudes. 

Therefore, because these arguments recommend a lower altitude 
than the planned flight altitude, I would like to learn more about the 
engineering reasons for going toward the higher altitudes. One 
crucial aspect could be the encounter of water droplets and ice crystals 
at these speeds. Another problem may be the larger heat input as 
compared with the radiation output from the skin. Of course there 
are other reasons. Whatever they may be, such engineering reasons 
are physical and open to cures according to the laws of science, while 
the disturbances to humans by sonic booms are emotional and psy- 
chological and, therefore, not so open to tradeoffs. Thus, with the 
limitations placed on sonic booms, the feasibility of lowering the flight 
altitude even by unconventional methods and the finding of a per- 
manent or state-of-the-art limit for economical and safe flying seems to 
be a necessary study. The minimum length requirements mentioned 
above for the higher altitudes are much beyond the planned or state- 
of-the-art lengths as governed by engineering constraints. 

The prerequisite for such considerations is, of course, the proof that 
the finite rise time is a great relief. I would prefer to have it demon- 
strated to my own ears rather than rely on statistical reports as to its 
effect. 
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